Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,603 T→G intergenic (+73/+80) 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,6030TG98.3% 173.0 / NA 61intergenic (+73/+80)yaaX/yaaADUF2502 domain‑containing protein YaaX/DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (30/29);  total (31/30)
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.34e-01

ACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGC  >  NC_000913/5367‑5846
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
aCTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAatat                                                                                                                                                                                                                                        >  1:117925/1‑250 (MQ=255)
    ggATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTCTGAATGGCGAGGCAATCGCTGGCACACACACGGACCGCCGCCTCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGGCCAGGCACAAATACCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTAAGCGCCTGATGCGACGCTGGCGCGTCCTAACCGGACCACGTGTATTCTGCAATATTTTg                                                                                                                                                                                                                                    >  1:386087/1‑250 (MQ=255)
            gTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGCTCATCACGGCGGACATGGTCCAGGCAAACCTCAACGATAAATGACCAATACCGGGGAACAATCCCGCATTCAGCGCACGGTGCGAAGCTGGCGCGTCTTCTCAGGCCTACGCGAATTCTGCAATATATTTTATCTGCAt                                                                                                                                                                                                                           >  1:375152/1‑251 (MQ=255)
                                   tggAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCCTTCAACGCCCGATGCGACGCTGGGGCGTCTTATCAGGCCTACGTGAATTCTGCCATATATTTAATCTGCATGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                                                                    >  1:138324/1‑251 (MQ=255)
                                   gggAAACAACATTATGAAGGGCGAGGCAACCGCTGCCACCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTGATCAGGCTTACGGGAAGTCGGCAATATATTGAATCTGCATGCTTGTGTAGGCCGGAGAAGGCg                                                                                                                                                                                                    <  2:93906/250‑1 (MQ=255)
                                      aaaCACCATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                 <  1:179172/251‑1 (MQ=255)
                                                 tGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGGCCAAGCAAAAAACACAGCTAAATGACAAAAGTCGGGCAACAATCCGGCCATTAGCCCCTGATGCGACGCTGGCGCGGCTTATCCAGGCTAAGGGCAATTCTGAATATATTGAATCTGCATGGTTTTGTAGGACGGGTAAGGAGTTCACGCCGCacac                                                                                                                                                                                      >  2:176286/1‑248 (MQ=255)
                                                  gAATGGCCAGGCAACCGCTGCCACCCACCCGGCCCGCCGCCCCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                     <  2:550629/251‑1 (MQ=39)
                                                     tGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTTAATCTGCATGCTTTTTTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGc                                                                                                                                                                                   >  1:149492/1‑250 (MQ=255)
                                                         gAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCCGGCAAACATCACCGCTAAATGACAAAAGCCGGGGAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGGATTTGCATGCTTTTTTAGGCCGGATAAGGCGGTCCCGCCGCCTCCGGCAttt                                                                                                                                                                               >  1:182571/1‑249 (MQ=255)
                                                          aGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGAc                                                                                                                                                                             <  1:286175/251‑1 (MQ=255)
                                                             cAATGGCTGGCCCCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCCCGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc                                                                                                                                                                          <  1:478938/251‑1 (MQ=38)
                                                                  gCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGCTCATCACGGCGGACCTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGGAACAATCCGGCATTAAGCGCCTGCTGCCACGCTGGCGCGTCTTATCAAGCCTCCGTGAATTCTTCCAAATATTGAATCTGCATGCTTTTGTTGGCCGGATAAGGCGTTTACGCCCCCTCCGGCCTTGACTGCCAACt                                                                                                                                                                     >  2:283693/1‑251 (MQ=255)
                                                                     ggCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTa                                                                                                                                                                   >  1:43710/1‑250 (MQ=255)
                                                                     gcgtccccaccggcccgccgCCCCCGCCCCGCCACCATAAGAAAGCTCCTCATGAGCAACACGTCGGTCATGGTCCAGGCAAACATCACCGATAAACGACAAATGCCGGGTAACAATCCGGCATTCAGCGTCTGAGGCTACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCGTCAATATATGGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGCAAACTTaa                                                                                                                                                                  <  2:219553/237‑1 (MQ=25)
                                                                          ccacacGGACCGCCGCCACCCCCTCGCCACCATAAGAACGCTCCTCATGATCATCACGGCGGCCATGGTCCAGTCAAACATCGCCGATAAATGCCAAATGCCGGGTAAAAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTAATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctg                                                                                                                                                             <  2:346485/249‑1 (MQ=38)
                                                                                 gacgcccgccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGTTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAACGCCTGATTCGACGCTGGCGCGTCTGATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGTTTTTGTAGGCCGTATAAGGCGTTCACGCCTCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAg                                                                                                                                                      <  2:89584/246‑1 (MQ=255)
                                                                                 gACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAAAGCCGGGGAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGGGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGGTCACGACGCATCCGGCATTTACTGCAAAATTTAAGCTGCTCgg                                                                                                                                                        >  1:604917/1‑248 (MQ=255)
                                                                                  aCCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATTCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCACTATATTGAATCTGCATGCCTTTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACACTGCTCGTAGc                                                                                                                                                     >  1:219502/1‑251 (MQ=255)
                                                                                      ccgccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGttt                                                                                                                                                 <  1:411276/251‑1 (MQ=255)
                                                                                        gccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCCAACATCAACGCTAAATGAAAAAAGACGGGAAACAATCCGGCAATCAGCGCCACATGCGACGCTGTCGCCGCTTTTAAGGCCTACCTGAATTCTGCAATATTTTGAATCTGCCTGCTTTTTTAGGCCTGATAAGGCGTTCCCGCCACATCCGGCATTTACTACaaaactcacgcttctcctaacgtttcc                                                                                                                                               >  1:356292/1‑227 (MQ=255)
                                                                                           aCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAaca                                                                                                                                            <  1:283693/251‑1 (MQ=255)
                                                                                               ccgcGCAACCATAAGAAAGCTCCACATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAg                                                                                                                                        <  2:411235/251‑1 (MQ=255)
                                                                                                       ccaTAAAAAAGCCCCCCATGATCATCACGGCGGTCATGGCCCAGGCAAACATCACCGCTAAATGACAAATGCCGGTTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCAt                                                                                                                                <  2:130876/251‑1 (MQ=39)
                                                                                                         aTAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTg                                                                                                                              <  2:67416/251‑1 (MQ=255)
                                                                                                          tAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCATGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAACGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTTTAGGCCGGGTAAGGCGGTCACGCCGCATACGGCCTTGACTGCAAAATTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGc                                                                                                                             >  2:546843/1‑251 (MQ=255)
                                                                                                            aGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTg                                                                                                                           <  1:529777/251‑1 (MQ=255)
                                                                                                               aaGCTCCTCATGATCATCACGGCGGTCAGGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAAGCCGGCATCCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTggag                                                                                                                        <  1:162828/251‑1 (MQ=255)
                                                                                                                    cctcATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATGCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGaag                                                                                                                    <  1:147303/250‑2 (MQ=255)
                                                                                                                          gagcatcaCGGCGGTCATGGTCCAGGCAAACATCCCCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCTACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTtcatc                                                                                                             <  2:227347/248‑1 (MQ=255)
                                                                                                                                      ggTCATGGTCCAGGCAAACATCACCGCTAAATTACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATGGCTGGAGTAAGCTTCATCAAAGAAGTAAcc                                                                                                 <  1:109328/251‑1 (MQ=255)
                                                                                                                                       ggcATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCt                                                                                                <  1:525338/249‑1 (MQ=255)
                                                                                                                                          aTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCg                                                                                             <  2:138324/251‑1 (MQ=255)
                                                                                                                                               ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTAtt                                                                                        >  2:343330/1‑251 (MQ=255)
                                                                                                                                               ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTAtt                                                                                        <  1:366763/251‑1 (MQ=255)
                                                                                                                                                 aGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTaa                                                                                      <  2:181074/251‑1 (MQ=255)
                                                                                                                                                     aaaCATCACCGCTAAATAACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                  >  1:443451/1‑251 (MQ=255)
                                                                                                                                                      aaCATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                  >  1:513330/1‑250 (MQ=255)
                                                                                                                                                         aTCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAACGAAGTAACCTTCGCTATTAAAACcagtc                                                                              >  2:128946/1‑251 (MQ=255)
                                                                                                                                                            cccGCTAAAGGACAAATGCCGGGTAACATTCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCTTACGTGAATTCTGCAATATATTGAATCTGCATGCTGTTGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCTTTTAAACTCCAGTTCGCCATTGCTGGATGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACctgtcagt                                                                           <  2:601037/250‑1 (MQ=255)
                                                                                                                                                                  aaaTGACAAATGCCGGGTAACAATCCGACATGCGGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATGGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCTGTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATGGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGctc                                                                      <  2:559210/250‑1 (MQ=255)
                                                                                                                                                                         aaaTGCCGGGTAACAATCCGGCATTCAGCGCATGATGCGAAGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAATCCAGTCAGTTGCTCTGGTTTgg                                                              >  2:77020/1‑251 (MQ=39)
                                                                                                                                                                          aaTGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGCATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAACCCAGTCAGTTGCTCTGGTTTGGt                                                             >  2:101231/1‑251 (MQ=255)
                                                                                                                                                                           aTGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGTCGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATGGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTGCGCTATTAAAACCAGTCAGTTGCGCTGTTTTGGTc                                                            <  2:568926/251‑1 (MQ=255)
                                                                                                                                                                            ggCCGGTTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATTCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTTTTAAAACCAGTCAGTTGCTCTGGGGTGGTCa                                                           <  2:386087/250‑1 (MQ=255)
                                                                                                                                                                             gCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCa                                                           <  2:182571/250‑1 (MQ=255)
                                                                                                                                                                                    aaCAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGAttt                                                   >  2:439937/1‑251 (MQ=255)
                                                                                                                                                                                           cGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCaata                                             >  1:588982/1‑250 (MQ=255)
                                                                                                                                                                                            ggCATTCAGAGCCTGATGCGTCGCTGGGGGGTCTTATCAGGCCTCGGTGCATTCTGCAATATATTGAATCTGCATGCTTTTGTAGCCCGGATAAGGCGTTAACGCAGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCCGTTCGCCATTGCTGGAGGAAGCTTAAACAAAGTAGTAACCTTAGCTATTAACACGAGTCTGTTGCTCTGGTTTGGTCATCCGATCTTCtataa                                            >  2:30322/1‑250 (MQ=255)
                                                                                                                                                                                                   cAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACg                                     >  1:444429/1‑251 (MQ=255)
                                                                                                                                                                                                        ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCATTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtca                                <  1:499542/251‑1 (MQ=255)
                                                                                                                                                                                                           gaggCGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCTTTAACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGGTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtcatca                             <  2:581382/248‑1 (MQ=255)
                                                                                                                                                                                                               cGACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGAcc                         >  2:563882/1‑251 (MQ=255)
                                                                                                                                                                                                                gACGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACcg                        >  1:121295/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGCTGGCGCGTCTTATCCGGCCTACGTGAATTACGCCATATATTGCATCTGCATGCTTTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGTATTGACTGCAAACTTAACGCTGCTCGTCGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGACGTCACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGCCTCATCAGACagcg                      >  2:30890/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACcgcg                      <  1:439937/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCTGGCGCGTCTTATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAACACGACTCATCAGACcgcat                     >  1:371235/1‑249 (MQ=255)
                                                                                                                                                                                                                              ttATCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGc          >  1:23864/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTCAGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATACTAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt        >  2:400483/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCCGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt        >  1:343904/1‑248 (MQ=255)
                                                                                                                                                                                                                                       ccTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGc  >  2:325508/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGCATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGC  >  NC_000913/5367‑5846

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