Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,587 A→C 11.2% intergenic (‑58/+35) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5870AC11.2% 481.6 / 56.4 161intergenic (‑58/+35)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (66/77);  new base C (10/8);  total (76/85)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

Only 100 of 161 total aligned reads displayed.

TGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATA  >  NC_000913/83441‑83734
                                                                                                                                                  |                                                                                                                                                    
tGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAt                                                                                                                                                    >  3:808699/1‑149 (MQ=255)
 gCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                                   >  1:311779/1‑149 (MQ=255)
   aGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgc                                                                                                                                                 >  3:1265014/1‑149 (MQ=255)
     tttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGGCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCCCAATGgcgc                                                                                                                                               >  4:964059/1‑149 (MQ=255)
      ttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGcgcg                                                                                                                                              <  4:510071/149‑1 (MQ=255)
         tctttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑t                                                                                                                                          <  2:1471892/149‑1 (MQ=255)
           tttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑ttt                                                                                                                                        >  4:9748/1‑149 (MQ=255)
           tttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑ttt                                                                                                                                        <  3:123728/149‑1 (MQ=255)
            ttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tttt                                                                                                                                       >  3:606870/1‑149 (MQ=255)
             tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttt                                                                                                                                      <  1:1289127/149‑1 (MQ=255)
             tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttt                                                                                                                                      <  2:651606/149‑1 (MQ=255)
             tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑ttttt                                                                                                                                      >  2:509175/1‑149 (MQ=255)
              cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tttttt                                                                                                                                     <  2:205151/149‑1 (MQ=255)
                cacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttttta                                                                                                                                   <  3:1460906/149‑3 (MQ=255)
                 acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtt                                                                                                                                  <  4:415035/149‑1 (MQ=255)
                  cTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCACTGTGGTATCGAAAATAATAACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgttt                                                                                                                                 >  3:162149/1‑149 (MQ=255)
                   tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtttg                                                                                                                                >  4:424717/1‑149 (MQ=255)
                   tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtttg                                                                                                                                >  1:616291/1‑149 (MQ=255)
                   tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgtttg                                                                                                                                >  1:821782/1‑149 (MQ=255)
                    cAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGa                                                                                                                               >  3:402155/1‑149 (MQ=255)
                         ttGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCg                                                                                                                          >  1:722466/1‑149 (MQ=255)
                          tGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGt                                                                                                                         <  3:849194/149‑1 (MQ=255)
                            ccTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGc                                                                                                                       <  2:799702/149‑1 (MQ=255)
                             cTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCt                                                                                                                      >  3:518389/1‑149 (MQ=255)
                               gTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑ttttt                                                                                                                                      >  2:778737/1‑131 (MQ=255)
                                tAACGCCTGTTCACCGTCGCGCAATGTGGTATCTAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGc                                                                                                                   <  3:742432/149‑1 (MQ=255)
                                    gCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttt                                                                                                                                        <  4:265219/124‑1 (MQ=255)
                                    gCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTaa                                                                                                               >  2:1409736/1‑149 (MQ=255)
                                     ccTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAAt                                                                                                              <  1:80304/149‑1 (MQ=255)
                                      cTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                                                                             <  4:606870/149‑1 (MQ=255)
                                        gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                                                                           <  1:813294/149‑1 (MQ=255)
                                        gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                                                                           <  3:1358684/149‑1 (MQ=255)
                                        gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                                                                           >  4:32769/1‑149 (MQ=255)
                                        gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                                                                           <  1:1477576/149‑1 (MQ=255)
                                        gTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                                                                           <  4:676396/149‑1 (MQ=255)
                                         ttCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                                                                          <  3:381984/149‑1 (MQ=255)
                                         ttCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                                                                          <  2:11653/149‑1 (MQ=255)
                                                     cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                                                                              <  4:493585/149‑1 (MQ=255)
                                                      aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACt                                                                                             >  2:100442/1‑149 (MQ=255)
                                                      aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACt                                                                                             >  1:153461/1‑149 (MQ=255)
                                                                 gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                                   <  3:1019486/85‑1 (MQ=255)
                                                                 gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                                   >  4:1019486/1‑85 (MQ=255)
                                                                 gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgt                                                                                                                                   <  2:599703/100‑1 (MQ=255)
                                                                     ataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaaa                                                                              <  3:424717/149‑1 (MQ=255)
                                                                       aatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaaga                                                                            <  4:574146/149‑1 (MQ=255)
                                                                        atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGGAaagaa                                                                           <  4:162149/149‑1 (MQ=255)
                                                                          gACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAagaaga                                                                         <  2:153461/149‑1 (MQ=255)
                                                                             ttGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAGGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                                                                      >  1:1001492/1‑149 (MQ=255)
                                                                                cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                                   <  2:863510/149‑1 (MQ=255)
                                                                                cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                                   <  1:1136107/149‑1 (MQ=255)
                                                                                 tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                  >  1:459426/1‑149 (MQ=255)
                                                                                  ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                                                                 >  4:366263/1‑149 (MQ=255)
                                                                                    cTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt                                                               <  2:388322/149‑1 (MQ=255)
                                                                                     tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                              >  2:1388834/1‑149 (MQ=255)
                                                                                     tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                              >  1:1143143/1‑149 (MQ=255)
                                                                                     tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                              <  3:1358887/149‑1 (MQ=255)
                                                                                        tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                           <  1:903976/149‑1 (MQ=255)
                                                                                        tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                           <  1:618340/149‑1 (MQ=255)
                                                                                           tttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagta                                                        >  1:1499349/1‑149 (MQ=255)
                                                                                              gggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                     >  1:245723/1‑149 (MQ=255)
                                                                                              gggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                     <  3:122846/149‑1 (MQ=255)
                                                                                               ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                                                                >  1:509181/1‑137 (MQ=255)
                                                                                                   cTTGTCTCTTTTAGAGCGCCTCGCTTCGGACATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATa                                                >  2:1345905/1‑149 (MQ=255)
                                                                                                    ttGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                                                                           <  2:944215/121‑1 (MQ=255)
                                                                                                    ttGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                                                                           >  1:944215/1‑121 (MQ=255)
                                                                                                      gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                             <  2:931433/149‑1 (MQ=255)
                                                                                                      gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                             <  4:1265014/149‑1 (MQ=255)
                                                                                                         tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGa                                          >  3:589714/1‑149 (MQ=255)
                                                                                                         tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGa                                          >  4:457935/1‑149 (MQ=255)
                                                                                                               agagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                    >  2:192180/1‑149 (MQ=255)
                                                                                                               agagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                    >  4:1408007/1‑149 (MQ=255)
                                                                                                                      cTCGCTTCGGGCGTAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                             >  1:1407363/1‑149 (MQ=255)
                                                                                                                      cTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                             <  2:1143143/149‑1 (MQ=255)
                                                                                                                       tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAt                            >  4:1440819/1‑149 (MQ=255)
                                                                                                                       tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAt                            <  2:1417728/149‑1 (MQ=255)
                                                                                                                        cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAtt                           >  1:964987/1‑149 (MQ=255)
                                                                                                                               ggCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                    >  1:680709/1‑149 (MQ=255)
                                                                                                                                 aaacaaaaaaCCCGCGCCATTGCGCGGGTTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                        >  2:407757/5‑128 (MQ=38)
                                                                                                                                 aaacaaaaaaCCCGCGCCATTGCGCGGGTTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                        <  1:407757/124‑1 (MQ=38)
                                                                                                                                   acaaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTCAGTAGTAGTAGACCGATAAAGCGAACGAGGTGAGTCATTAAATCAGCTcc                >  1:512123/3‑149 (MQ=255)
                                                                                                                                    aaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                   >  4:1062210/1‑129 (MQ=255)
                                                                                                                                    aaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                   <  3:1062210/129‑1 (MQ=255)
                                                                                                                                     aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTCAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg              >  2:645044/1‑149 (MQ=255)
                                                                                                                                      aaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGa             <  3:451832/149‑1 (MQ=255)
                                                                                                                                        aaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa               >  4:1453410/1‑145 (MQ=255)
                                                                                                                                        aaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa               <  3:1453410/145‑1 (MQ=255)
                                                                                                                                        aaaCCCGCGCAATGGCGCGGG‑TGTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg           >  1:200678/1‑149 (MQ=255)
                                                                                                                                         aaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                 >  2:153020/1‑126 (MQ=255)
                                                                                                                                          aCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa         >  2:1039000/1‑149 (MQ=255)
                                                                                                                                          aCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACCGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa         >  4:84651/1‑149 (MQ=255)
                                                                                                                                          aCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACCGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTATACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa         >  3:1031484/1‑149 (MQ=255)
                                                                                                                                            ccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATg       >  1:72357/1‑149 (MQ=255)
                                                                                                                                            ccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATg       >  1:519628/1‑149 (MQ=255)
                                                                                                                                              gcgcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCg     >  1:1023167/1‑149 (MQ=255)
                                                                                                                                              gcgcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCg     >  2:502451/1‑149 (MQ=255)
                                                                                                                                                gcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGat   <  1:192180/149‑1 (MQ=255)
                                                                                                                                                gcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGa    <  1:1202829/148‑1 (MQ=255)
                                                                                                                                                gcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGa    >  2:1202829/1‑148 (MQ=255)
                                                                                                                                                 cAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                   <  1:367106/116‑1 (MQ=255)
                                                                                                                                                 cAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgt                                   >  2:367106/1‑116 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATA  >  NC_000913/83441‑83734

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑