Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 70,659 G→A W91* (TGG→TGA araC → DNA‑binding transcriptional dual regulator AraC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091370,6590GA98.2% 195.0 / NA 55W91* (TGG→TGAaraCDNA‑binding transcriptional dual regulator AraC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (31/23);  minor base C (1/0);  total (32/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCC  >  NC_000913/70416‑70904
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
gccggtttattcgtttacactccctCTGGTTGCGGTTTAAACGCCGATTGAGTCCAACGGTTGTCTCGTTTTTTTTAGCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAAGCTCACCATCCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGATAGTTTTTCTGCCGCCGGGGTTATAGTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGATGCTCGCGAATGATATCAc                                                                                                                                                                                                                                                  <  1:308925/226‑1 (MQ=255)
 cctggATTCTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGTTGATATTTTGCTGGTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATAGCACCa                                                                                                                                                                                                                                                <  2:421370/248‑1 (MQ=255)
                             tggCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGCCCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGAGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCCCTACGGTCGTCATCCGGCGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTAc                                                                                                                                                                                                                    >  1:244600/1‑251 (MQ=255)
                                 gggTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGc                                                                                                                                                                                                                <  1:215085/251‑1 (MQ=255)
                                                 gAGGCCAGCGGTTATCTCGATTTTTTTATCGACCGACCGGTGGAAATGAAAGGTTATATTCTCAATCCCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGATAATTTGTCTGCCTACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCTTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTGCGTCCGCGCGCCTACTGGCATGAATGGCTTAACt                                                                                                                                                                                                 <  1:445395/250‑1 (MQ=255)
                                                 gAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTg                                                                                                                                                                                                <  1:397852/251‑1 (MQ=255)
                                                 gAGGCAAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACt                                                                                                                                                                                                 <  1:15157/250‑1 (MQ=255)
                                                           gttttctcGATTTTTTTATCGACGGGCCGCTGGGAATGAAAGGTTATATTCACAATCTCACCATGCTCGGGCAGGGGGTGTTGAAAAATCAGGGACTAGAATTTGCCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTGCGGCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAata                                                                                                                                                                                      <  2:33278/247‑1 (MQ=255)
                                                              aTCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGCATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCAGTCAATAttt                                                                                                                                                                                   >  1:238438/1‑251 (MQ=255)
                                                                   gATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCa                                                                                                                                                                                                                                                   >  2:387730/1‑182 (MQ=255)
                                                                   gATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCa                                                                                                                                                                                                                                                   <  1:387730/182‑1 (MQ=255)
                                                                        ttttATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCTCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACg                                                                                                                                                                          <  2:240009/250‑1 (MQ=255)
                                                                        tttgtacgaccgacCGCTGGGAATGAAAGGTTTTATTCTCAATCTCACCATTCGCGGTCAGGGGGTGTTGAAAAATCAGGGACGAGATTTTGTCTGCCGCCCGGGTGATATTTTGCGTTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCATACTGGCATGAATGGCTTAACTGGCCGTCAATCTTTGCCAATACg                                                                                                                                                                          <  2:231683/244‑1 (MQ=255)
                                                                          ttATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGt                                                                                                                                                                       >  1:68585/1‑251 (MQ=255)
                                                                                 ccgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGtttctttc                                                                                                                                                                <  1:411550/251‑1 (MQ=255)
                                                                                  cgaccgCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGATGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACAAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCg                                                                                                                                                               >  1:296015/1‑251 (MQ=255)
                                                                                       gCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCgg                                                                                                                                                          >  1:267365/1‑251 (MQ=255)
                                                                                        cTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCgg                                                                                                                                                          >  1:182368/1‑250 (MQ=255)
                                                                                                 aaaGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATAAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGAGTAAgcgc                                                                                                                                                 >  2:457070/1‑250 (MQ=255)
                                                                                                     gTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCAcc                                                                                                                                              >  2:163841/1‑249 (MQ=255)
                                                                                                        atatTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCg                                                                                                                                         >  1:164199/1‑251 (MQ=255)
                                                                                                               cAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTc                                                                                                                                   >  1:468529/1‑250 (MQ=255)
                                                                                                                aaTCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTc                                                                                                                                   >  2:130229/1‑249 (MQ=255)
                                                                                                                      atcATTCGCGGGAAGGGGGTGGAGAAATATCGGGGCAGAGAATTTGTCTGCAGAGCGGGTGATGTTTTGCTGTTCCCGCCAGGAGAGATTCATCACAACGGTCGTCACCCGGAGGCGCGCGAATGATATCACCAGTGGGTTTACGTTCGTCCGCGCGCCTGCTGGCATAAATGGCTTAACTGGCCGTCACTATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTGCAGCGAccg                                                                                                                          <  2:381331/249‑2 (MQ=255)
                                                                                                                      aCCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCt                                                                                                                          >  1:61623/1‑251 (MQ=255)
                                                                                                                         aTTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGtt                                                                                                                       >  1:406134/1‑251 (MQ=255)
                                                                                                                            cgctgccAGGGGGTGGTGAAAAAGCAGGGACGAAAATTTGTCTGCCGACCGGGGGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTgg                                                                                                                    <  2:497291/245‑1 (MQ=255)
                                                                                                                                    gggtgagatgcaaagtCAGTGACGAGACTTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCCATACGGTCGTCATTCGGAGGCTCTCGCATGATATCACCAGTGGGGTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTAGCCGTCAATATTTTCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCCCATTTCTAGGCCCTGTTTGGGCAAATcc                                                                                                            >  2:323645/16‑250 (MQ=255)
                                                                                                                                     ggggTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCAt                                                                                                           >  2:255976/1‑251 (MQ=255)
                                                                                                                                                    cAGGGACGAGAATTTGTCTGCCGACCGGGTGGTATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGCTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTCACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCCGGCAAgg                                                                                            >  1:253596/1‑251 (MQ=255)
                                                                                                                                                     aGGGCCGAGATTTTGTCTGCCGCCCGGGTGTTATTTTGCGGTTCCCGCCAGGAGAGATTCATCACTACGGGCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAgtg                                                                                           <  2:193448/251‑3 (MQ=255)
                                                                                                                                                     aGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGttt                                                                                                                      >  1:419159/1‑224 (MQ=255)
                                                                                                                                                      gggccgagtttttGTGTGCCGCCCGGGTGATTTTTTTCTGGTCCCGCCAGGAGAGATTCATCACTACGTTCGTTATCCGGAGGCTCGCGAATGATATCACCAGTGTGTTTACGTTCGTCCGCGCTCCTACTGGCATTAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGACGCGCACCAGCCGCATTTCAGCGACCTGTTTGTGCAAATCTTTAACGCCGGGCAAgggg                                                                                          <  2:195818/241‑1 (MQ=255)
                                                                                                                                                               aTTTTTCTGCCGACCGGGTTATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGACTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTATTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGTCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCt                                                                                 <  1:403620/251‑1 (MQ=255)
                                                                                                                                                                          gACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGCAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGACTGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCCGCGACCTGTTTGGGCAAATCATTAACGCCGGGCACGGGGAAGGGCGCTATTCGGAgctg                                                                      >  1:56655/1‑251 (MQ=255)
                                                                                                                                                                            cagggTGATCTTTTGCTGTTCACGCCAGGAGAGATTCATCACTACGGTCGTAATCCGGAGGCTCGCGAATGATATCACAAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGACGTAAATATTTGCCAATACGGGTTTCTTTCGCCCGGTTGAAGCGCACCAGCCGCATTTCCGCGACCTGTTTTGGCATATCATTAACTCCGGGCAAGGGGAAGGGCGCTATTCGGAgctgct                                                                    >  2:331754/3‑251 (MQ=255)
                                                                                                                                                                             cGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGctgctg                                                                   <  1:255976/251‑1 (MQ=255)
                                                                                                                                                                              gggTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCg                                                                                                                              >  1:490982/1‑191 (MQ=255)
                                                                                                                                                                               ggTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTTCTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGc                                                                 <  1:114955/251‑1 (MQ=255)
                                                                                                                                                                               ggTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTAttag                                                                            >  2:352806/1‑238 (MQ=255)
                                                                                                                                                                               ggTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCg                                                                            <  1:352806/240‑1 (MQ=255)
                                                                                                                                                                                               tCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTg                                                   >  2:32216/1‑249 (MQ=255)
                                                                                                                                                                                                          gagaTTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACAAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACt                                      >  1:299059/1‑251 (MQ=255)
                                                                                                                                                                                                                 atcaCTACGGTCGTCATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGgcgc                               >  1:195328/1‑251 (MQ=255)
                                                                                                                                                                                                                          gtcgtcATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCa                                                                                                                                             <  1:145033/133‑1 (MQ=255)
                                                                                                                                                                                                                          gtcgtcATCCGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCa                                                                                                                                             >  2:145033/1‑133 (MQ=255)
                                                                                                                                                                                                                                  ccGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACg               >  1:269732/1‑250 (MQ=255)
                                                                                                                                                                                                                                  ccGGAGGCTCGCGAATGATATCACCAGTGGGTTACCTTTCGTCCGCGCGCATACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGGTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCACGGGGAAGGGCGCTATTCGGCGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACg               >  2:449974/1‑250 (MQ=255)
                                                                                                                                                                                                                                   cGGAGGCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAg             <  2:61623/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gCTCGCGAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCt        >  2:139916/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cgcgAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTa                                        <  1:469624/216‑1 (MQ=255)
                                                                                                                                                                                                                                           cgcgAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTa                                        >  2:469624/1‑216 (MQ=255)
                                                                                                                                                                                                                                              gAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCtccatcc  >  2:193993/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gAATGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCtccatcc  <  1:29110/251‑1 (MQ=255)
                                                                                                                                                                                                                                               aaTGATATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCtccatcc  >  1:359585/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GCCGGGATACTCGTTTAACGCCCATCTGGTGGCGGGTTTAACGCCGATTGAGGCCAACGGTTATCTCGATTTTTTTATCGACCGACCGCTGGGAATGAAAGGTTATATTCTCAATCTCACCATTCGCGGTCAGGGGGTGGTGAAAAATCAGGGACGAGAATTTGTCTGCCGACCGGGTGATATTTTGCTGTTCCCGCCAGGAGAGATTCATCACTACGGTCGTCATCCGGAGGCTCGCGAATGGTATCACCAGTGGGTTTACTTTCGTCCGCGCGCCTACTGGCATGAATGGCTTAACTGGCCGTCAATATTTGCCAATACGGGTTTCTTTCGCCCGGATGAAGCGCACCAGCCGCATTTCAGCGACCTGTTTGGGCAAATCATTAACGCCGGGCAAGGGGAAGGGCGCTATTCGGAGCTGCTGGCGATAAATCTGCTTGAGCAATTGTTACTGCGGCGCATGGAAGCGATTAACGAGTCGCTCCATCC  >  NC_000913/70416‑70904

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