Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,951,559:1 +A 28.8% pseudogene (18/663 nt) ilvG → pseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,5591.A28.8% 309.5 / 241.5 184pseudogene (18/663 nt)ilvGpseudogene, acetolactate synthase 2 large subunit, valine‑insensitive; acetolactate synthase II, large subunit, cryptic, interrupted
Reads supporting (aligned to +/- strand):  ref base . (76/55);  new base A (26/27);  total (102/82)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 197 total aligned reads displayed.

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCG  >  NC_000913/3951460‑3951658
                                                                                                    |                                                                                                   
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGTAACACTgc‑g                                                                                                    >  1:826266/1‑100 (MQ=255)
  cAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgc‑gc                                                                                                   <  1:1703990/100‑1 (MQ=255)
     gCTGCGCCAGGCACATGTGGCATTACAAGGTGATTTACATGCTCTGTTACCAGCATTACAGCAGCAGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCa                                                                                                >  1:1369125/1‑100 (MQ=255)
       tGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGc                                                                                              <  2:447576/100‑1 (MQ=255)
        gCGTCGGGGACATGTGGCAGTACAATGTGATGTAAATGCTCTGGGACCAGCATTACAGCAGCCGGTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCt                                                                                             >  1:3680410/1‑100 (MQ=255)
        gCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCt                                                                                             >  1:50578/1‑100 (MQ=255)
          gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTg                                                                                            <  2:2141009/100‑1 (MQ=255)
           tCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGc                                                                                           >  1:4257349/1‑100 (MQ=255)
            cAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTCGCAGCAACACTGC‑GCGCAGCTGCGt                                                                                         >  1:1450878/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                        <  2:3330253/100‑1 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtga                                                                                       >  2:3278255/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtga                                                                                       >  2:3562017/1‑100 (MQ=255)
                cacaTTTTGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtgatg                                                                                     <  1:4133457/100‑1 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtgatg                                                                                     >  2:2082036/1‑100 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtgatg                                                                                     >  2:3799347/1‑100 (MQ=255)
                 acaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtgatg                                                                                     <  2:4115549/100‑1 (MQ=255)
                 acaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGtgatga                                                                                    <  2:2514773/100‑1 (MQ=255)
                  caTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGaa                                                                                   <  1:2453167/100‑1 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGaaga                                                                                 >  1:872217/1‑98 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACa                                                                                 <  1:3278237/100‑1 (MQ=255)
                      tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAt                                                                                >  1:1749144/1‑100 (MQ=255)
                      tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACAtt                                                                               >  2:1543375/1‑100 (MQ=255)
                       ggCATTACAAGGTTATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAtt                                                                               <  2:2188207/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCt                                                                            >  1:19527/1‑100 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTg                                                                           <  2:2233131/100‑1 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTg                                                                           <  2:3262934/100‑1 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTgg                                                                          >  2:320382/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTccagat                                                                         >  2:1713182/1‑96 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGc                                                                         >  1:1614203/1‑100 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGc                                                                         <  2:434312/100‑1 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        <  1:2762628/100‑1 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        <  2:1808415/100‑1 (MQ=255)
                                aGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                     <  1:718107/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                    >  2:3938103/1‑100 (MQ=255)
                                  gTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACg                                                                   <  1:3084134/100‑1 (MQ=255)
                                   tGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                  >  1:2028894/1‑100 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                                >  1:3963935/1‑100 (MQ=255)
                                            tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCACAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                         <  2:3355240/100‑1 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GGGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        >  2:1761569/1‑100 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        <  2:427191/100‑1 (MQ=255)
                                               tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        <  1:1879275/99‑1 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     >  1:2060044/1‑100 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     >  2:4060813/1‑100 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     <  1:519840/100‑1 (MQ=255)
                                                     aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGCCGCGAt                                                >  1:1764571/1‑100 (MQ=255)
                                                       cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCTTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  2:2777955/1‑100 (MQ=255)
                                                          cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                            >  2:1151351/1‑100 (MQ=255)
                                                              aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                       >  1:3214551/1‑100 (MQ=255)
                                                               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                      >  1:925913/1‑100 (MQ=255)
                                                                agcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                      >  2:1013610/1‑100 (MQ=255)
                                                                agcagcCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  2:3220697/100‑1 (MQ=255)
                                                                 gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  2:3478796/100‑1 (MQ=255)
                                                                 gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  2:4257381/100‑1 (MQ=255)
                                                                   agcCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  >  2:2281969/1‑100 (MQ=255)
                                                                    tctgTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  <  2:1233405/97‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 <  2:501116/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 <  2:3424127/100‑1 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGCTGAACATTCCTGGCGATACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               >  2:526500/1‑100 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  2:1694690/100‑1 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  1:3799320/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               >  2:2803781/1‑100 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              <  2:2060051/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              >  2:884950/1‑100 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              <  2:4383125/100‑1 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             <  2:4362813/100‑1 (MQ=255)
                                                                         tAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             >  2:3992418/1‑100 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                            >  2:1924347/1‑100 (MQ=255)
                                                                          aaaTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaa                            >  2:64723/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaaa                          >  2:2036847/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaaa                          >  1:1100082/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACCCTATCTACGCGCCGTTGTTGTTaaaa                          >  2:2669377/1‑100 (MQ=255)
                                                                            aTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                           <  1:1029857/99‑1 (MQ=255)
                                                                            aTCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                         <  2:3913909/100‑1 (MQ=255)
                                                                            aTCAATGACTGGCAGCAACACTGC‑GCGCACCTGCGTGATGAACCATCCTGGCCTTACGACCATCCCGCTGACGCCATCTACGCGCCGTTGTTGTTAacac                         >  1:3039844/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                       >  2:4385959/1‑100 (MQ=255)
                                                                               aaTGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                      >  1:2931046/1‑100 (MQ=255)
                                                                                 tGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                    <  2:4375640/100‑1 (MQ=255)
                                                                                 tGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                    >  2:1458869/1‑100 (MQ=255)
                                                                                 tGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                    <  1:3246919/100‑1 (MQ=255)
                                                                                   aCTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                  >  1:2371449/1‑100 (MQ=255)
                                                                                   aCTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                  <  2:830735/100‑1 (MQ=255)
                                                                                    cTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 >  2:3198643/1‑100 (MQ=255)
                                                                                    cTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:3888355/100‑1 (MQ=255)
                                                                                    cTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:1785440/100‑1 (MQ=255)
                                                                                     tGGCAGAAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:1082584/100‑1 (MQ=255)
                                                                                      ggCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa               >  2:365070/1‑100 (MQ=255)
                                                                                        cagcaACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             >  2:2911741/1‑100 (MQ=255)
                                                                                        cagcaACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             >  2:1605524/1‑100 (MQ=255)
                                                                                         agcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             <  1:209385/100‑1 (MQ=255)
                                                                                          gcaACACTGC‑GCGCAGCTGCGTGATGAACACTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAACACAACTGTCGGATCGt           >  1:3481641/1‑100 (MQ=255)
                                                                                           caACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa          <  2:349525/100‑1 (MQ=255)
                                                                                           caACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa          >  2:2259126/1‑100 (MQ=255)
                                                                                            aaCACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa         <  1:2129879/100‑1 (MQ=255)
                                                                                             acacTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa         <  2:1701347/100‑1 (MQ=255)
                                                                                              cacTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGATGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa        >  2:409667/1‑100 (MQ=255)
                                                                                               acTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGTTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc      >  2:90917/1‑100 (MQ=255)
                                                                                                cTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt     >  2:3142315/1‑100 (MQ=255)
                                                                                                  gc‑gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGc   <  2:243014/100‑1 (MQ=255)
                                                                                                   c‑gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg  <  2:3713624/100‑1 (MQ=255)
                                                                                                   c‑gcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg  >  1:3564882/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGC‑GCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCG  >  NC_000913/3951460‑3951658

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑