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.

AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCG  >  NC_000913/3951460‑3951658
                                                                                                    |                                                                                                   
 aCAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGTAACACTgcg                                                                                                    >  1:826266/1‑100 (MQ=255)
  cAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTgcgc                                                                                                   <  1:1703990/100‑1 (MQ=255)
     gCTGCGCCAGGCACATGTGGCATTACAAGGTGATTTACATGCTCTGTTACCAGCATTACAGCAGCAGTTAAATCAATGACTGGCAGCAACACTGCGCGCa                                                                                                >  1:1369125/1‑100 (MQ=255)
       tGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGc                                                                                              <  2:447576/100‑1 (MQ=255)
        gCGTCGGGGACATGTGGCAGTACAATGTGATGTAAATGCTCTGGGACCAGCATTACAGCAGCCGGTAAATCAATGACTGGCAGCAACACTGCGCGCAGCt                                                                                             >  1:3680410/1‑100 (MQ=255)
        gCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCt                                                                                             >  1:50578/1‑100 (MQ=255)
          gTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTg                                                                                            <  2:2141009/100‑1 (MQ=255)
           tCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGc                                                                                           >  1:4257349/1‑100 (MQ=255)
            cAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTCGCAGCAACACTGCGCGCAGCTGCGt                                                                                         >  1:1450878/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtg                                                                                        <  2:3330253/100‑1 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtga                                                                                       >  2:3278255/1‑100 (MQ=255)
              ggCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtga                                                                                       >  2:3562017/1‑100 (MQ=255)
                cacaTTTTGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtgatg                                                                                     <  1:4133457/100‑1 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtgatg                                                                                     >  2:2082036/1‑100 (MQ=255)
                cacaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtgatg                                                                                     >  2:3799347/1‑100 (MQ=255)
                 acaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGtgatg                                                                                     <  2:4115549/100‑1 (MQ=255)
                 acaTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGtgatga                                                                                    <  2:2514773/100‑1 (MQ=255)
                  caTGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGaa                                                                                   <  1:2453167/100‑1 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGaaga                                                                                 >  1:872217/1‑98 (MQ=255)
                    tgtgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACa                                                                                 <  1:3278237/100‑1 (MQ=255)
                      tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAt                                                                                >  1:1749144/1‑100 (MQ=255)
                      tgGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACAtt                                                                               >  2:1543375/1‑100 (MQ=255)
                       ggCATTACAAGGTTATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACAtt                                                                               <  2:2188207/100‑1 (MQ=255)
                         cATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCt                                                                            >  1:19527/1‑100 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTg                                                                           <  2:2233131/100‑1 (MQ=255)
                          aTTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTg                                                                           <  2:3262934/100‑1 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTgg                                                                          >  2:320382/1‑100 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTccagat                                                                         >  2:1713182/1‑96 (MQ=255)
                            tACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGc                                                                         >  1:1614203/1‑100 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGc                                                                         <  2:434312/100‑1 (MQ=255)
                             aCAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        <  1:2762628/100‑1 (MQ=255)
                              cAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCg                                                                        <  2:1808415/100‑1 (MQ=255)
                                aGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTa                                                                     <  1:718107/100‑1 (MQ=255)
                                 ggTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTAc                                                                    >  2:3938103/1‑100 (MQ=255)
                                  gTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACg                                                                   <  1:3084134/100‑1 (MQ=255)
                                   tGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGa                                                                  >  1:2028894/1‑100 (MQ=255)
                                     aTTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGAcc                                                                >  1:3963935/1‑100 (MQ=255)
                                            tGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCACAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCgg                                                         <  2:3355240/100‑1 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGGGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        >  2:1761569/1‑100 (MQ=255)
                                             gCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        <  2:427191/100‑1 (MQ=255)
                                               tctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGt                                                        <  1:1879275/99‑1 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     >  1:2060044/1‑100 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     >  2:4060813/1‑100 (MQ=255)
                                                ctGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGAc                                                     <  1:519840/100‑1 (MQ=255)
                                                     aCCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGCCGCGAt                                                >  1:1764571/1‑100 (MQ=255)
                                                       cAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCTTGGCGTTACGACCATCCCGGTGACGCTATCt                                              >  2:2777955/1‑100 (MQ=255)
                                                          cATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTAc                                            >  2:1151351/1‑100 (MQ=255)
                                                              aCAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGcc                                       >  1:3214551/1‑100 (MQ=255)
                                                               cagcagCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                      >  1:925913/1‑100 (MQ=255)
                                                                agcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCg                                      >  2:1013610/1‑100 (MQ=255)
                                                                agcagcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  2:3220697/100‑1 (MQ=255)
                                                                 gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  2:3478796/100‑1 (MQ=255)
                                                                 gcagcCGTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgt                                     <  2:4257381/100‑1 (MQ=255)
                                                                   agcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  >  2:2281969/1‑100 (MQ=255)
                                                                    tctgTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgt                                  <  2:1233405/97‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 <  2:501116/100‑1 (MQ=255)
                                                                    gcCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgtt                                 <  2:3424127/100‑1 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGCTGAACATTCCTGGCGATACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               >  2:526500/1‑100 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  2:1694690/100‑1 (MQ=255)
                                                                      cGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               <  1:3799320/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgt                               >  2:2803781/1‑100 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              <  2:2060051/100‑1 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              >  2:884950/1‑100 (MQ=255)
                                                                       gTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCgttgttgtt                              <  2:4383125/100‑1 (MQ=255)
                                                                        ttAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTa                             <  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)
                                                                           aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaaa                          >  2:2036847/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaaa                          >  1:1100082/1‑100 (MQ=255)
                                                                           aaTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACCCTATCTACGCGCCGTTGTTGTTaaaa                          >  2:2669377/1‑100 (MQ=255)
                                                                            aTCAATGACTGGCAGCAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTaaa                           <  1:1029857/99‑1 (MQ=255)
                                                                            aTCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaac                         <  2:3913909/100‑1 (MQ=255)
                                                                            aTCAATGACTGGCAGCAACACTGCGCGCACCTGCGTGATGAACCATCCTGGCCTTACGACCATCCCGCTGACGCCATCTACGCGCCGTTGTTGTTAacac                         >  1:3039844/1‑100 (MQ=255)
                                                                              cAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaa                       >  2:4385959/1‑100 (MQ=255)
                                                                               aaTGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAaacaac                      >  1:2931046/1‑100 (MQ=255)
                                                                                 tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                    <  2:4375640/100‑1 (MQ=255)
                                                                                 tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                    >  2:1458869/1‑100 (MQ=255)
                                                                                 tGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTg                    <  1:3246919/100‑1 (MQ=255)
                                                                                   aCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                  >  1:2371449/1‑100 (MQ=255)
                                                                                   aCTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTc                  <  2:830735/100‑1 (MQ=255)
                                                                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 >  2:3198643/1‑100 (MQ=255)
                                                                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:3888355/100‑1 (MQ=255)
                                                                                    cTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:1785440/100‑1 (MQ=255)
                                                                                     tGGCAGAAACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCg                 <  1:1082584/100‑1 (MQ=255)
                                                                                      ggCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGa               >  2:365070/1‑100 (MQ=255)
                                                                                        cagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             >  2:2911741/1‑100 (MQ=255)
                                                                                        cagcaACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             >  2:1605524/1‑100 (MQ=255)
                                                                                         agcaACACTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATc             <  1:209385/100‑1 (MQ=255)
                                                                                          gcaACACTGCGCGCAGCTGCGTGATGAACACTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAACACAACTGTCGGATCGt           >  1:3481641/1‑100 (MQ=255)
                                                                                           caACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa          <  2:349525/100‑1 (MQ=255)
                                                                                           caACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTa          >  2:2259126/1‑100 (MQ=255)
                                                                                            aaCACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa         <  1:2129879/100‑1 (MQ=255)
                                                                                             acacTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaa         <  2:1701347/100‑1 (MQ=255)
                                                                                              cacTGCAGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGATGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTaaa        >  2:409667/1‑100 (MQ=255)
                                                                                               acTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGTTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAAcc      >  2:90917/1‑100 (MQ=255)
                                                                                                cTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCt     >  2:3142315/1‑100 (MQ=255)
                                                                                                  gcgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGc   <  2:243014/100‑1 (MQ=255)
                                                                                                   cgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg  <  2:3713624/100‑1 (MQ=255)
                                                                                                   cgcgcAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCg  >  1:3564882/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AACAAGCTGCGTCAGGCACATGTGGCATTACAAGGTGATTTAAATGCTCTGTTACCAGCATTACAGCAGCCGTTAAATCAATGACTGGCAGCAACACTGCGCGCAGCTGCGTGATGAACATTCCTGGCGTTACGACCATCCCGGTGACGCTATCTACGCGCCGTTGTTGTTAAAACAACTGTCGGATCGTAAACCTGCG  >  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: