Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,451,569 T→G 9.3% M1290L (ATG→CTG)  rpoB ← RNA polymerase subunit beta

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0077793,451,5690T→G 9.3% 69.7 ‑0.7 0.0246M1290L (ATG→CTG) rpoBRNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  new base (9/14):  ref base (105/117):  total (115/131)
Fisher's exact test strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 258 total aligned reads displayed.

AAGGATTCTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGC  >  NC_007779/3451469‑3451669
                                                                                                    |                                                                                                    
aaGGATTCTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCAt                                                                                                      >  1:5144632/1‑101 (MQ=255)
  ggATTCTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCAttt                                                                                                    <  1:597914/101‑1 (MQ=255)
   gATTCTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTc                                                                                                   <  1:8496539/101‑1 (MQ=255)
    aTTCTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTcc                                                                                                  <  1:6707404/101‑1 (MQ=255)
       cTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGc                                                                                               <  1:8862591/101‑1 (MQ=255)
        tGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCa                                                                                              >  1:10011395/1‑101 (MQ=255)
         ggCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAg                                                                                             <  1:11679794/101‑1 (MQ=255)
          gCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAgg                                                                                            >  1:5274541/1‑101 (MQ=255)
            aTGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGgt                                                                                          >  1:13970143/1‑101 (MQ=255)
             tGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGgtg                                                                                         >  1:3024430/1‑101 (MQ=255)
              gCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGgtgt                                                                                        <  1:12318293/101‑1 (MQ=255)
               cccGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTa                                                                                       <  1:1849569/101‑1 (MQ=255)
                 cGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtg                                                                                     <  1:11901768/101‑1 (MQ=255)
                  ggCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgc                                                                                    >  1:1637904/1‑101 (MQ=255)
                  ggCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgc                                                                                    <  1:1319054/101‑1 (MQ=255)
                  ggCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgc                                                                                    >  1:3519297/1‑101 (MQ=255)
                    cTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgctg                                                                                  <  1:310465/101‑1 (MQ=255)
                     tCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgctgc                                                                                 >  1:1906966/1‑101 (MQ=255)
                     tCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgctgc                                                                                 <  1:13418838/101‑1 (MQ=255)
                     tCCATCTGATGGCTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTAtgctgc                                                                                 >  1:1364755/1‑101 (MQ=255)
                        aTCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGcc                                                                              <  1:9382220/101‑1 (MQ=255)
                         tCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCg                                                                             <  1:12118955/101‑1 (MQ=255)
                          cTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGt                                                                            <  1:4736292/101‑1 (MQ=255)
                           tGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTa                                                                           >  1:1766577/1‑101 (MQ=255)
                           tGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTa                                                                           >  1:6237050/1‑101 (MQ=255)
                             aTGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATg                                                                         <  1:12634520/101‑1 (MQ=255)
                             aTGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATg                                                                         <  1:5010807/101‑1 (MQ=255)
                                gTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCtt                                                                      >  1:12993739/1‑101 (MQ=255)
                                 gtGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTc                                                                     <  1:10182652/100‑1 (MQ=255)
                                  tGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTcc                                                                    >  1:10811598/1‑101 (MQ=255)
                                   gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                   <  1:3579643/101‑1 (MQ=255)
                                   gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                   <  1:9694574/101‑1 (MQ=255)
                                   gCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCa                                                                   <  1:12400068/101‑1 (MQ=255)
                                    ccGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAg                                                                  <  1:4799761/101‑1 (MQ=255)
                                     cGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAgc                                                                 >  1:4480438/1‑101 (MQ=255)
                                         cACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGccc                                                             >  1:6956173/1‑101 (MQ=255)
                                         cACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGccc                                                             >  1:3850691/1‑101 (MQ=255)
                                          aCGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCca                                                            <  1:8209153/101‑1 (MQ=255)
                                          aCGATCTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCca                                                            <  1:7967346/101‑1 (MQ=255)
                                           cGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCcac                                                           >  1:12811254/1‑101 (MQ=255)
                                            gATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCcaca                                                          <  1:9291161/101‑1 (MQ=255)
                                               gTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACtt                                                       <  1:11221180/101‑1 (MQ=255)
                                                   ttATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCAt                                                   <  1:9178193/101‑1 (MQ=255)
                                                    tataCATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCCCACTTCCAtc                                                  >  1:405542/1‑101 (MQ=255)
                                                    tataCATCTTAGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCAtc                                                  >  1:9197732/1‑101 (MQ=255)
                                                      taCATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCAtctc                                                <  1:663187/101‑1 (MQ=255)
                                                         aTCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcccc                                             <  1:6255800/101‑1 (MQ=255)
                                                         aTCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTcccc                                             <  1:2344784/101‑1 (MQ=255)
                                                           cTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGa                                           >  1:4169222/1‑101 (MQ=255)
                                                            ttGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGaa                                          >  1:6303915/1‑101 (MQ=255)
                                                              ggTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACATTTCCATCTCCCCGAAAc                                        <  1:11317152/101‑1 (MQ=255)
                                                              ggTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAAc                                        <  1:11482843/101‑1 (MQ=255)
                                                                 aCGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCt                                     >  1:14087124/1‑101 (MQ=255)
                                                                  cGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTg                                    <  1:860525/101‑1 (MQ=255)
                                                                     ccGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacc                                 >  1:490180/1‑101 (MQ=255)
                                                                     ccGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacc                                 <  1:8406974/101‑1 (MQ=255)
                                                                      cGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacca                                >  1:2819029/1‑101 (MQ=255)
                                                                      cGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacca                                <  1:7057871/101‑1 (MQ=255)
                                                                      cGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGacca                                >  1:7022047/1‑101 (MQ=255)
                                                                        ttCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccacc                              >  1:5218806/1‑101 (MQ=255)
                                                                        ttCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGaccacc                              <  1:2802301/101‑1 (MQ=255)
                                                                          cACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGa                            >  1:7796507/1‑101 (MQ=255)
                                                                           aCGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGaa                           <  1:11468572/101‑1 (MQ=255)
                                                                            cGTCATCGGACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAAc                          <  1:6570421/101‑1 (MQ=255)
                                                                            cGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAAc                          >  1:4203476/1‑101 (MQ=255)
                                                                             gTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACt                         >  1:10561546/1‑101 (MQ=255)
                                                                             gTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACt                         >  1:10431583/1‑101 (MQ=255)
                                                                              tcatcgGACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACtg                        <  1:12918595/101‑1 (MQ=255)
                                                                                atcgGACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACtgag                      <  1:717515/101‑3 (MQ=255)
                                                                                 tcggACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGc                     <  1:2659638/98‑1 (MQ=255)
                                                                                 tcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGc                     >  1:9453790/1‑101 (MQ=255)
                                                                                 tcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGc                     <  1:9091286/101‑1 (MQ=255)
                                                                                 tcaGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGc                     <  1:7512857/101‑1 (MQ=255)
                                                                                   ggACTTGATGGTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCt                   <  1:3310051/100‑1 (MQ=255)
                                                                                   aGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCt                   >  1:10309089/1‑101 (MQ=255)
                                                                                    gACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCtt                  <  1:13003023/101‑1 (MQ=255)
                                                                                    gACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCtt                  <  1:7998329/101‑1 (MQ=255)
                                                                                     aCTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTa                 >  1:1966938/1‑101 (MQ=255)
                                                                                       ttAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTacc               >  1:4548567/1‑101 (MQ=255)
                                                                                        tAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTacca              <  1:9479115/101‑1 (MQ=255)
                                                                                        tAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTacca              <  1:10993415/101‑1 (MQ=255)
                                                                                         aaCGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTaccac             <  1:13369897/101‑1 (MQ=255)
                                                                                         aaCGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTaccac             >  1:8702131/1‑101 (MQ=255)
                                                                                           cGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCAccc           >  1:4399922/1‑101 (MQ=255)
                                                                                            ggTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCa          <  1:5660813/101‑1 (MQ=255)
                                                                                            ggTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCa          <  1:8655425/101‑1 (MQ=255)
                                                                                             gTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAg         >  1:2193434/1‑101 (MQ=255)
                                                                                             gTGAGCAGCTCCTGCAGGGTGTAGGCAGCGCCGTATGCCTGCATTGCCCACACCTCCATCTCACCGAAGCGCTGACCACCGAACTGAGCCTTACCACCCAg         <  1:12693554/101‑1 (MQ=255)
                                                                                              tGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGc        <  1:12575999/101‑1 (MQ=255)
                                                                                               gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg       <  1:13807084/101‑1 (MQ=255)
                                                                                               gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg       >  1:2565329/1‑101 (MQ=255)
                                                                                               gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg       <  1:10578869/101‑1 (MQ=255)
                                                                                               gAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCg       <  1:10922588/101‑1 (MQ=255)
                                                                                                aGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCgg      >  1:4256667/1‑101 (MQ=255)
                                                                                                 gCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgc     >  1:587552/1‑101 (MQ=255)
                                                                                                  cATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgct    <  1:7134941/101‑1 (MQ=255)
                                                                                                  cATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgct    <  1:14262544/101‑1 (MQ=255)
                                                                                                   aTTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgctg   <  1:7347439/101‑1 (MQ=255)
                                                                                                   aTTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgctg   <  1:6511984/101‑1 (MQ=255)
                                                                                                    tttCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGgctgc  <  1:6087654/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
AAGGATTCTGGCATGCCCGGCTCCATCTGATGGTTGCCGTCCACGATGTTTTTATACATCTTGGTACGACCGTTCACGTCATCAGACTTAACGGTGAGCATTTCCTGCAGGGTGTATGCTGCGCCGTATGCTTCCAGCGCCCACACTTCCATCTCCCCGAAACGCTGACCACCGAACTGTGCCTTACCACCCAGCGGCTGC  >  NC_007779/3451469‑3451669

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.