Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,094,489 G→A 5.3% intergenic (‑217/+234) rhaD ← / ← rhaA rhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,094,4890GA 5.3% 340.4 / 7.8 131intergenic (‑217/+234)rhaD/rhaArhamnulose‑1‑phosphate aldolase/L‑rhamnose isomerase
Reads supporting (aligned to +/- strand):  ref base G (58/66);  new base A (2/5);  total (60/71)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.85e-01

Only 100 of 131 total aligned reads displayed.

ATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCT  >  NC_000913/4094343‑4094636
                                                                                                                                                  |                                                                                                                                                   
aTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTg                                                                                                                                                  <  1:667331/150‑1 (MQ=32)
  ccGGCAATCACTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCATACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGcc                                                                                                                                                <  2:139837/150‑1 (MQ=17)
   cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                               <  1:37785/150‑1 (MQ=32)
   cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCg                                                                                                                                               >  2:859689/1‑150 (MQ=32)
    ggCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCgt                                                                                                                                              >  1:1547281/1‑149 (MQ=17)
       aagcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCt                                                                                                                                           <  2:1785784/147‑1 (MQ=17)
            atGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAgg                                                                                                                                      <  1:1454576/150‑1 (MQ=34)
             tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGc                                                                                                                                     >  2:1500194/1‑150 (MQ=34)
             tGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGc                                                                                                                                     >  2:596911/1‑150 (MQ=34)
                cttattcGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCgg                                                                                                                                  <  1:1871923/144‑1 (MQ=21)
                cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCgg                                                                                                                                  <  2:431507/150‑1 (MQ=34)
                  taTGCGACGCTGTCGCGTCTTATCAGGCCTACTACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGAt                                                                                                                                <  1:764235/149‑1 (MQ=21)
                  gATGCTACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGAt                                                                                                                                <  1:378082/150‑1 (MQ=21)
                   ttGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATa                                                                                                                               <  1:1809253/149‑1 (MQ=21)
                      cGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAgg                                                                                                                            <  2:1849845/150‑1 (MQ=34)
                      cGACGCTGGCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAgg                                                                                                                            <  1:1562755/150‑1 (MQ=21)
                             gTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtgt                                                                                                                     >  1:1640454/1‑150 (MQ=37)
                             gTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTtgt                                                                                                                     >  1:2025091/1‑150 (MQ=37)
                                gcgTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                                                                                                                  <  1:372331/150‑1 (MQ=37)
                                  gTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGc                                                                                                                <  2:1713155/150‑1 (MQ=37)
                                  gTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGc                                                                                                                <  2:1704632/150‑1 (MQ=37)
                                   tCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCa                                                                                                               <  2:874861/150‑1 (MQ=37)
                                    cTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGGCTACAACTGTTGCCGACCGGTAGGCCGGAGACGGCACGTGTGCCGCag                                                                                                              >  1:138559/1‑149 (MQ=11)
                                    cTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCAt                                                                                                              >  1:977204/1‑150 (MQ=37)
                                     ttATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCCTc                                                                                                             >  1:1040039/1‑150 (MQ=37)
                                       aTCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGACTGATGCGACGCTGACGCGTCTTATCAGGCCTACAACTGTTGCCGACCAGTAGGCCGGATAAGGCTCTTGTCCCGCATCCg                                                                                                           >  1:690464/1‑150 (MQ=11)
                                           ggCCAACAACTATTGCCGACCTGTAGGCATGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGCCGCTGTCGCGTCTTATCAGGCCGACAACTGTTGCCGACCAGTACGCACGATAAGGCACTTGTGACCCATCAGGcaa                                                                                                       >  1:469338/1‑150 (MQ=11)
                                            gCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaat                                                                                                      <  1:791702/150‑1 (MQ=37)
                                              cTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCGGTAGGCCGGATAAGGCACTTGGGCCGCATCCGGcaatca                                                                                                    >  1:1390876/1‑150 (MQ=255)
                                                aCAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACAGTTGCCGACCTGTAGGACGCATAAGGCACTTGTGCCGCATCCGGcaatcaat                                                                                                  >  2:1018430/1‑150 (MQ=255)
                                                 cAACTATTGACGACCTGTAGGCTTGATAAGGCACTGGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                                                                 <  1:1071062/150‑1 (MQ=255)
                                                 cAACTATAGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGACAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCGGAATCCGGCATTCAATg                                                                                                 >  2:1024523/1‑150 (MQ=11)
                                                 aaaCTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATg                                                                                                 <  2:932653/149‑1 (MQ=37)
                                                         gCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGc                                                                                         <  1:530053/150‑1 (MQ=34)
                                                         accGACCTGTAGGCCTGATAAGGCACTTCTGCCGCATCCGCCAATCAATGCCTGATGCCACGCAGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCAGTAGGCCGGATAAGGCACATGTGACGCAAACGGCAATCAATGCCAGATGc                                                                                         >  1:1011102/2‑150 (MQ=11)
                                                           cGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGa                                                                                       <  2:355947/150‑1 (MQ=34)
                                                             aCCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACg                                                                                     <  1:1472816/150‑1 (MQ=34)
                                                              ccTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGc                                                                                    >  2:2096961/1‑150 (MQ=34)
                                                                tGTAGGCCTGATAAGGCACTTGTGCCGCCTCAGGCAATCAATGCCTGATGAGACGCTGTCGCGTCTTATCAGGCCTACAACTGGTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCATTCATTGCCTGATGCGACGCTg                                                                                  >  2:322841/1‑150 (MQ=11)
                                                                tGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTg                                                                                  >  1:551242/1‑150 (MQ=34)
                                                                tGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTg                                                                                  >  2:475154/1‑150 (MQ=34)
                                                                tGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGAAGCTGTCGCGTCTTATCAGGCATACAACAGTTGCCTACATGTAGGCAGGATAAGGACCTTGTGCCGCATCCGGCAATCAATGCCTGAGGCGACGCTg                                                                                  >  2:322844/1‑150 (MQ=11)
                                                                 gTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGt                                                                                 >  1:1169989/1‑150 (MQ=34)
                                                                    ggCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTcgc                                                                              >  2:1298515/1‑150 (MQ=34)
                                                                       cTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                           <  1:396591/150‑1 (MQ=33)
                                                                       cTGAGATGGCTCTTGGGCCGCATCCGGCAATCACTGCCTGATGCGTGGCGGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTc                                                                           <  1:396585/150‑1 (MQ=11)
                                                                        tGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCt                                                                          <  2:1482060/150‑1 (MQ=32)
                                                                          aTAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTa                                                                        >  2:784575/1‑150 (MQ=32)
                                                                           tAAGGCACTTGTGCCGCATCCGGGAATCAATGCCTGATGCGTCGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTAt                                                                       <  2:512121/150‑1 (MQ=17)
                                                                            aaGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGc                                                                                                                           >  2:1248920/1‑97 (MQ=33)
                                                                            aaGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGc                                                                                                                           <  1:1248920/97‑1 (MQ=33)
                                                                            aaGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATc                                                                      <  2:1412865/150‑1 (MQ=32)
                                                                                 ccTTGGGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGcc                                                                 <  1:447313/149‑1 (MQ=17)
                                                                                   ttGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTa                                                               <  2:2004628/150‑1 (MQ=32)
                                                                                    tgtgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc                                                                                >  2:142247/1‑132 (MQ=32)
                                                                                    tgtgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTc                                                                                <  1:142247/132‑1 (MQ=32)
                                                                                     gtgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACa                                                             >  2:542508/1‑150 (MQ=32)
                                                                                     gcgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACa                                                             <  2:1569811/148‑1 (MQ=17)
                                                                                      tgCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACaa                                                            >  2:104298/1‑150 (MQ=32)
                                                                                       gCCGCATCCGGCAATCTCTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACTACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                           <  1:1008728/150‑1 (MQ=11)
                                                                                       gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                           >  2:442634/1‑150 (MQ=32)
                                                                                       gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                           <  2:230253/150‑1 (MQ=32)
                                                                                       gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                           >  2:2037571/1‑150 (MQ=32)
                                                                                       gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                           >  1:1569267/1‑150 (MQ=32)
                                                                                       gCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTAGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAAc                                                           >  2:545658/1‑150 (MQ=17)
                                                                                            aTCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTg                                                      <  2:582700/150‑1 (MQ=34)
                                                                                              ccGGCAATCATTGCCTGATGCGACGCTGTCGCGTCTTATCAGGCGTACTACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGcc                                                    <  2:1665497/150‑1 (MQ=255)
                                                                                               cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACACCTATTGCCg                                                   >  2:1970023/1‑150 (MQ=21)
                                                                                               cGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGAGTCTTATCAGGCCTACAACTATTGCCg                                                   >  2:561819/1‑150 (MQ=21)
                                                                                                 gCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAc                                                 >  2:1841169/1‑150 (MQ=34)
                                                                                                  caatcaatGCCTGATGCGACGCTGTCGCGTCTTGTCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAcc                                                <  2:1808677/150‑1 (MQ=21)
                                                                                                  caatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACGACTATTGCCGAcc                                                <  2:466592/150‑1 (MQ=21)
                                                                                                  caatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGAcc                                                >  2:888126/1‑150 (MQ=34)
                                                                                                  caatcaatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGG‑‑AAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTg                                              <  2:954129/150‑1 (MQ=255)
                                                                                                      caatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTa                                            >  2:1432043/1‑150 (MQ=37)
                                                                                                       aatGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAg                                           <  2:721063/150‑1 (MQ=37)
                                                                                                       aatGCCTGATGCGACGCTGTCGAGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCTTCTTATCAGGCATACAACTATTGCCGACCTGTAg                                           <  1:1381429/150‑1 (MQ=255)
                                                                                                           ccTGATGCGACGCTGGCGCGTCTTTTCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCg                                       <  2:417350/150‑1 (MQ=255)
                                                                                                            cTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCTCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCgg                                      <  2:1003651/150‑1 (MQ=17)
                                                                                                                tGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCCGTAGGCCGGATaa                                  >  1:1814780/1‑150 (MQ=17)
                                                                                                                 gCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAg                                 <  2:1233478/150‑1 (MQ=17)
                                                                                                                 gCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAg                                 <  1:787715/150‑1 (MQ=17)
                                                                                                                  cGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAgg                                <  1:2037571/150‑1 (MQ=37)
                                                                                                                   gACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGc                               <  1:924392/150‑1 (MQ=37)
                                                                                                                          tCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTtgtg                        <  1:764266/150‑1 (MQ=37)
                                                                                                                           cgcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCCGGCCTACAACTATTGCCGACCGGTAGGCCGGATAAGGCACTAGTGc                       >  1:240607/1‑150 (MQ=255)
                                                                                                                            gcgTCTTATCAGGCCTAGAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                      <  1:564290/150‑1 (MQ=37)
                                                                                                                            gcgTCTTATCAGGCCTACTACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTTTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                      <  2:477433/150‑1 (MQ=255)
                                                                                                                            gcgTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGcc                      <  2:1970815/150‑1 (MQ=37)
                                                                                                                                  tATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACGATTGCCGACCTGTAGGCCGGATAAGGCCCTTGTGCCGCATcc                >  2:2031131/1‑150 (MQ=255)
                                                                                                                                   aTCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCg               <  2:600228/150‑1 (MQ=37)
                                                                                                                                     cAGGCCTACAACTGTTGCTGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTCCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGc             >  2:1639423/1‑150 (MQ=255)
                                                                                                                                     cAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGc             >  2:1663497/1‑150 (MQ=37)
                                                                                                                                      aGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCCATGCCTGATGCGACGCTGTCGCGACTTATCACGCCTACAACTATTGCCGACCAGTAGGCCCGATTATGCACTAGAGCCGCAACCGCca            >  1:66803/1‑150 (MQ=11)
                                                                                                                                          cTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGcaatca        >  2:648169/1‑150 (MQ=37)
                                                                                                                                          cTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCAGTCGCGTCTTATCAGGCCTACAACAATTGCCGACCGGTAGGCCGGATAAGGCCCTAGAGCCGCATCCGGcaatca        >  1:1218001/1‑150 (MQ=11)
                                                                                                                                           tACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTATTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCCGcaatcaa       >  1:871638/1‑150 (MQ=255)
                                                                                                                                               aCTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGcc   >  1:913757/1‑150 (MQ=37)
                                                                                                                                                cTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCt  >  2:835862/1‑150 (MQ=37)
                                                                                                                                                cTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCAGTAGGCCGGATAAGGCGCTAGTGCCGCATCCGGCAAACAATGcca  >  1:72253/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCTGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTGTTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTCGCGTCTTATCAGGCCTACAACTATTGCCGACCTGTAGGCCGGATAAGGCACTTGTGCCGCATCCGGCAATCAATGCCT  >  NC_000913/4094343‑4094636

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