Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,642 A→C intergenic (+112/+41) yaaX → / ← yaaA DUF2502 domain‑containing protein YaaX/DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,6420AC93.6% 128.6 / ‑1.0 48intergenic (+112/+41)yaaX/yaaADUF2502 domain‑containing protein YaaX/DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (16/28);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-01

GGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCG  >  NC_000913/5396‑5878
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
ggcgggtgaAAACAACATTATGAATGGCGGGGCAACCGTGGGCCCCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGAt                                                                                                                                                                                                                                           <  2:421047/247‑1 (MQ=255)
    ggttgAAACAACATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATaa                                                                                                                                                                                                                                         <  1:125621/245‑1 (MQ=255)
             aacaTTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGACCTGGGCCAGGCAAACATCACCGCCAAATGACAAATGCCGGGGAACAATCCCGCATTCAGCGCCCGATGCGACGCTGGCGCGGCTTTTCAGGCCTACGTGAATTTTTCCATATATTTAATCTGCATGCTTTTTTTAGACGGATAAGGCGTTCAcgc                                                                                                                                                                                                                              >  1:215961/1‑251 (MQ=255)
              acaTTATGAATGGCGAGGCAATCGCTGGCCCCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                                                                                                                                             <  1:193845/251‑1 (MQ=255)
                     gAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCCGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTGCGTGAATTCTGCAATATATTTAATCTGCATGCTTTTGTAGGCCGGATAAGGCGGTCACGCCGCATCCg                                                                                                                                                                                                                      >  1:460848/1‑251 (MQ=39)
                                cAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCCTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGGAAAAATCCGGCATTCAGCGCCCGATGAGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGAATGGTTTTTTAGTACGGATAAGGGGTTCCCGCCGCATCCGGCCATGACTGc                                                                                                                                                                                                           >  2:336895/1‑251 (MQ=255)
                                                 acggcccccccccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCCCGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCg                                                                                                                                                                                          <  2:402508/240‑1 (MQ=255)
                                                       ccagccccccgccgCGCCACCATAAGAAAGCTCCTCAGGATCATCCCGGCTGTCATGGTCCAGGCAAACATCACCGCGAAATGCCAAAGGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACGTGAATTCGGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGt                                                                                                                                                                                    <  2:76186/243‑1 (MQ=37)
                                                                               aGAAAGCTCCCCAAGATCATCACGGCGTTCATGGTCCAGCCAAACATCACCGCTAAATGACAAAGGCCGGGAAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTg                                                                                                                                                            <  2:92059/251‑1 (MQ=38)
                                                                                  aaGCTCCTCAGGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCTCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGAGTCGACGCTGACGCGTCTTATCAGGCCTACGTGAATTCTGCAATAAATTGAATCTGCATGCTTTTGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGAAGTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTggag                                                                                                                                                         <  1:304209/251‑1 (MQ=255)
                                                                                              atcatcACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATTACAAAAGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGGGTCTTATCAAGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAAGCCGGATAAGGCGTTCACGCCGCATCCGGCCTTGACTGCCAAATTAACGCTGCTCGTAACGTTTAAACCCCAGTTCGCCCTTGCTTGGGGAAGCTtcatc                                                                                                                                              >  2:149252/1‑250 (MQ=255)
                                                                                                 atcaCGGCGGTCCTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCCTGTGCGACGCTGGCGCGTCTTATCAGGCCTACGTGGATTCTGCAATATATTGGAACTGCATGCTTTTGTTGGCCGGGTAAAGCGTTCCCGCCCCATCCCGCATTTTTTGCACACTTAACGCCGCTCGTATCGTTTTAACACCCGTTCGCCACTTCTGGAGGGAGCTTCATCaaa                                                                                                                                           >  2:95002/1‑250 (MQ=255)
                                                                                                 aacaCGGCGGTCATGGTCCAGGCAAACATCCCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAGGCTTCATCAaag                                                                                                                                          <  2:141466/249‑1 (MQ=255)
                                                                                                              tGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGc                                                                                                                             >  1:14356/1‑251 (MQ=255)
                                                                                                               ggTCCAGGCAAACATCACCGCTAAATGTCAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGTTGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGTCATTGACTGCAAACTGAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCt                                                                                                                            <  2:69473/251‑1 (MQ=38)
                                                                                                                      gCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTaaaa                                                                                                                    <  1:85385/251‑1 (MQ=255)
                                                                                                                       cAAACATCACCGCTAAATGACAAAGGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAc                                                                                                                   <  2:98296/251‑1 (MQ=255)
                                                                                                                         aaCATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGTTAAGGCGTTCACGCCGCATCCGGCATTGACTGCAACCTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                                                  >  2:296565/1‑250 (MQ=255)
                                                                                                                              cACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTTCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGTCTGCAAACTTAACGCTGCTCGTAGCGTTTAAACCCCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACcagtcag                                                                                                            <  1:439835/251‑1 (MQ=39)
                                                                                                                               ccccatAAATTACAAATGCCGGGTAACAATCCGGCATTCAGCTCCTTATGCGACGCTGGCTCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCTGCATTGACTGCAAACTTAACTCTGCTCGTAGCGTTTAAACGCCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACcagtcagt                                                                                                           <  1:367707/246‑1 (MQ=37)
                                                                                                                                 cGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTg                                                                                                         >  1:325564/1‑251 (MQ=255)
                                                                                                                                       aTGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCGGAGGCGACGCGGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCGGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTgg                                                                                                   <  2:124790/251‑1 (MQ=39)
                                                                                                                                            acaTGCCGGGTAACAAACAGGCATTCAGCGCCTGATGCGAAGCTGGCGCGTCTTATAAGGCCTACGAGAATTCTGCAATATATTGAATCTGCATGCGTTTGTAGGCAGGATAAGGCGTGCACGCCGCATCCGGCATTGACTGCAAACTTAAAGCTGCTCGTAGCGTTTAAACAACAGGTAGCCATTGCTGGAGGAAGCTTAATAAAAGAAGAAACCATAGAATGTAAAGCCAGTCAGTGGATCTGGTTTgg                                                                                              >  2:321392/3‑251 (MQ=255)
                                                                                                                                             aaTGCCGGGTAAAAATCCGTCATTCGGCGCCTGATGCGACTCTGGCGCGTTTTATCAGTCGTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCAGCAAAGAAGGAACCTTCTCTATTAAAACCAGTCAGGTGCTCTTGTTTGGt                                                                                             <  2:228191/251‑1 (MQ=37)
                                                                                                                                                     gTAACAATCCGGCATTCAGCGCCTGATTCGCCGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGAt                                                                                     <  1:96595/251‑1 (MQ=255)
                                                                                                                                                     gTAAAAATCCGGCTTTCAGCGCCTGAGGCGACGCTGGCGGGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGAt                                                                                     <  1:8848/251‑1 (MQ=39)
                                                                                                                                                        aCAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGATCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGAtttt                                                                                  <  1:361870/251‑1 (MQ=255)
                                                                                                                                                         cAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTc                                                                                 >  1:256757/1‑251 (MQ=255)
                                                                                                                                                               ggCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGATGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCaataat                                                                           >  1:413063/1‑251 (MQ=255)
                                                                                                                                                                     cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAg                                                                                          >  1:195469/1‑230 (MQ=255)
                                                                                                                                                                     cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAg                                                                                          <  2:195469/230‑1 (MQ=255)
                                                                                                                                                                        cgcCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACt                                                                  <  1:336895/251‑1 (MQ=255)
                                                                                                                                                                           cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGATCAGCCGATTTTCAATAATAAAACGACtcat                                                               >  1:367893/1‑251 (MQ=255)
                                                                                                                                                                              gATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACTCCAGTTCGCCATTTCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtcat                                                               <  1:53679/249‑1 (MQ=255)
                                                                                                                                                                                 gCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATGGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGGCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGAc                                                          <  2:215961/251‑1 (MQ=255)
                                                                                                                                                                                  cGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGAcc                                                         >  1:165730/1‑251 (MQ=255)
                                                                                                                                                                                   gccGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGAc                                                                   >  1:318833/3‑240 (MQ=255)
                                                                                                                                                                                   gagtcTGGCTCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGAc                                                                   <  2:318833/236‑1 (MQ=38)
                                                                                                                                                                                      gCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATCTTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGt                                                     >  1:181019/1‑251 (MQ=255)
                                                                                                                                                                                          gcgcGTCTTATCAGGCTTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCtt                                                 <  1:149252/251‑1 (MQ=255)
                                                                                                                                                                                             cgTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGACGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCt                                              >  1:165775/1‑251 (MQ=255)
                                                                                                                                                                                               tCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCtt                                             <  1:155971/250‑1 (MQ=255)
                                                                                                                                                                                                   aTCAGGCCTGCGTGAATTCTGCAATATATGGAATCTGCATGCTTTTGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCg                                         <  2:460848/250‑1 (MQ=255)
                                                                                                                                                                                                     ccggCCGACGTGAATTCTGCAATATTTTGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCGCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCAGCAGACCGCGTGCTTTCTTAGCGTAg                                      <  2:256757/249‑1 (MQ=255)
                                                                                                                                                                                                         ccTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGc                                  >  1:18648/1‑251 (MQ=255)
                                                                                                                                                                                                                              atTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCg              <  1:95002/250‑1 (MQ=255)
                                                                                                                                                                                                                                   atcggCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGttcttc        <  2:325564/247‑1 (MQ=255)
                                                                                                                                                                                                                                         aTGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATTATCTTAAATTTGCCGTTCTTCTCATCg  >  1:409818/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             
GGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGCTGATGATCTTAAATTTGCCGTTCTTCTCATCG  >  NC_000913/5396‑5878

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: