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

ACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGC  >  NC_000913/5367‑5846
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
aCTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAatat                                                                                                                                                                                                                                        >  1:117925/1‑250 (MQ=255)
    ggATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTCTGAATGGCGAGGCAATCGCTGGCACACACACGGACCGCCGCCTCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGGCCAGGCACAAATACCCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTAAGCGCCTGATGCGACGCTGGCGCGTCCTAAC‑CGGACCACGTGTATTCTGCAATATTTTg                                                                                                                                                                                                                                    >  1:386087/1‑250 (MQ=255)
            gTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGCTCATCACGGCGGACATGGTCCAGGCAAACCTCAACGATAAATGACCAATACCGGGGAACAATCCCGC‑ATTCAGCGCACGGTGCGAAGCTGGCGCGTCTTCTC‑AGGCCTACGCGAATTCTGCAATATATTTTATCTGCAt                                                                                                                                                                                                                           >  1:375152/1‑251 (MQ=255)
                                   tggAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑CTTCAACGCCCGATGCGACGCTGGGGCGTCTTATC‑AGGCCTACGTGAATTCTGCCATATATTTAATCTGCATGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                                                                    >  1:138324/1‑251 (MQ=255)
                                   gggAAACAACATTATGAAGGGCGAGGCAACCGCTGCCACCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTGATC‑AGGCTTACGGGAAGTCGGCAATATATTGAATCTGCATGCTTGTGTAGGCCGGAGAAGGCg                                                                                                                                                                                                    <  2:93906/250‑1 (MQ=255)
                                      aaaCACCATTATGAAGGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                 <  1:179172/251‑1 (MQ=255)
                                                 tGAATGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGGCATGGGCCAAGCAAAAAACACAGCTAAATGACAAAAGTCGGGCAACAATCCGGCCATT‑AGCCCCTGATGCGACGCTGGCGCGGCTTATCCAGGCTAAGGGCAATTCTG‑AATATATTGAATCTGCATGGTTTTGTAGGACGGGTAAGGAGTTCACGCCGCacac                                                                                                                                                                                      >  2:176286/1‑248 (MQ=255)
                                                  gAATGGCCAGGCAACCGCTGCCACCCACCCGGCCCGCCGCCCCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                     <  2:550629/251‑1 (MQ=39)
                                                     tGGCGAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTTAATCTGCATGCTTTTTTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGc                                                                                                                                                                                   >  1:149492/1‑250 (MQ=255)
                                                         gAGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCCGGCAAACATCACCGCTAAATGACAAAAGCCGGGGAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGGATTTGCATGCTTTTTTAGGCCGGATAAGGCGGTCCCGCCGCCTCCGGCAttt                                                                                                                                                                               >  1:182571/1‑249 (MQ=255)
                                                          aGGCAATCGCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGAc                                                                                                                                                                             <  1:286175/251‑1 (MQ=255)
                                                             cAATGGCTGGCCCCCACCCGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCCCGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGc                                                                                                                                                                          <  1:478938/251‑1 (MQ=38)
                                                                  gCTGGCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGCTCATCACGGCGGACCTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGGAACAATCCGGC‑ATTAAGCGCCTGCTGCCACGCTGGCGCGTCTTATC‑AAGCCTCCGTGAATTCTTCCAAATATTGAATCTGCATGCTTTTGTTGGCCGGATAAGGCGTTTACGCCCCCTCCGGCCTTGACTGCCAACt                                                                                                                                                                     >  2:283693/1‑251 (MQ=255)
                                                                     ggCACCCACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTa                                                                                                                                                                   >  1:43710/1‑250 (MQ=255)
                                                                     gcgtccccaccggcccgccgCCCCCGCCCCGCCACCATAAGAAAGCTCCTCATGAGCAACACGTCGGTCATGGTCCAGGCAAACATCACCGATAAACGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGTCTGAGGCTACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCGTCAATATATGGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGCAAACTTaa                                                                                                                                                                  <  2:219553/237‑1 (MQ=25)
                                                                          ccacacGGACCGCCGCCACCCCCTCGCCACCATAAGAACGCTCCTCATGATCATCACGGCGGCCATGGTCCAGTCAAACATCGCCGATAAATGCCAAATGCCGGGTAAAAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTAATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACgctg                                                                                                                                                             <  2:346485/249‑1 (MQ=38)
                                                                                 gacgcccgccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGTTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAACGCCTGATTCGACGCTGGCGCGTCTGATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGTTTTTGTAGGCCGTATAAGGCGTTCACGCCTCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAg                                                                                                                                                      <  2:89584/246‑1 (MQ=255)
                                                                                 gACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAAAGCCGGGGAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGGGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGGTCACGACGCATCCGGCATTTACTGCAAAATTTAAGCTGCTCgg                                                                                                                                                        >  1:604917/1‑248 (MQ=255)
                                                                                  aCCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATTCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCACTATATTGAATCTGCATGCCTTTTTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACACTGCTCGTAGc                                                                                                                                                     >  1:219502/1‑251 (MQ=255)
                                                                                      ccgccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGttt                                                                                                                                                 <  1:411276/251‑1 (MQ=255)
                                                                                        gccACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCCAACATCAACGCTAAATGAAAAAAGACGGGAAACAATCCGGC‑AATCAGCGCCACATGCGACGCTGTCGCCGCTTTTA‑AGGCCTACCTGAATTCTGCAATATTTTGAATCTGCCTGCTTTTTTAGGCCTGATAAGGCGTTCCCGCCACATCCGGCATTTACTACaaaactcacgcttctcctaacgtttcc                                                                                                                                               >  1:356292/1‑227 (MQ=255)
                                                                                           aCCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAaca                                                                                                                                            <  1:283693/251‑1 (MQ=255)
                                                                                               ccgcGCAACCATAAGAAAGCTCCACATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAg                                                                                                                                        <  2:411235/251‑1 (MQ=255)
                                                                                                       ccaTAAAAAAGCCCCCCATGATCATCACGGCGGTCATGGCCCAGGCAAACATCACCGCTAAATGACAAATGCCGGTTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCAt                                                                                                                                <  2:130876/251‑1 (MQ=39)
                                                                                                         aTAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTg                                                                                                                              <  2:67416/251‑1 (MQ=255)
                                                                                                          tAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCATGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAACGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTTTAGGCCGGGTAAGGCGGTCACGCCGCATACGGCCTTGACTGCAAAATTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGc                                                                                                                             >  2:546843/1‑251 (MQ=255)
                                                                                                            aGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTg                                                                                                                           <  1:529777/251‑1 (MQ=255)
                                                                                                               aaGCTCCTCATGATCATCACGGCGGTCAGGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAAGCCGGC‑ATCCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTggag                                                                                                                        <  1:162828/251‑1 (MQ=255)
                                                                                                                    cctcATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATGCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGaag                                                                                                                    <  1:147303/250‑2 (MQ=255)
                                                                                                                          gagcatcaCGGCGGTCATGGTCCAGGCAAACATCCCCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCTACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTtcatc                                                                                                             <  2:227347/248‑1 (MQ=255)
                                                                                                                                      ggTCATGGTCCAGGCAAACATCACCGCTAAATTACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATGGCTGGAGTAAGCTTCATCAAAGAAGTAAcc                                                                                                 <  1:109328/251‑1 (MQ=255)
                                                                                                                                       ggcATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCt                                                                                                <  1:525338/249‑1 (MQ=255)
                                                                                                                                          aTGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCg                                                                                             <  2:138324/251‑1 (MQ=255)
                                                                                                                                               ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTAtt                                                                                        >  2:343330/1‑251 (MQ=255)
                                                                                                                                               ccAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTAtt                                                                                        <  1:366763/251‑1 (MQ=255)
                                                                                                                                                 aGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTaa                                                                                      <  2:181074/251‑1 (MQ=255)
                                                                                                                                                     aaaCATCACCGCTAAATAACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                  >  1:443451/1‑251 (MQ=255)
                                                                                                                                                      aaCATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAAcc                                                                                  >  1:513330/1‑250 (MQ=255)
                                                                                                                                                         aTCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAACGAAGTAACCTTCGCTATTAAAACcagtc                                                                              >  2:128946/1‑251 (MQ=255)
                                                                                                                                                            cccGCTAAAGGACAAATGCCGGGTAACATTCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCTTACGTGAATTCTGCAATATATTGAATCTGCATGCTGTTGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCTTTTAAACTCCAGTTCGCCATTGCTGGATGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACctgtcagt                                                                           <  2:601037/250‑1 (MQ=255)
                                                                                                                                                                  aaaTGACAAATGCCGGGTAACAATCCGAC‑ATGCGGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATGGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCTGTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATGGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGctc                                                                      <  2:559210/250‑1 (MQ=255)
                                                                                                                                                                         aaaTGCCGGGTAACAATCCGGC‑ATTCAGCGCATGATGCGAAGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAATCCAGTCAGTTGCTCTGGTTTgg                                                              >  2:77020/1‑251 (MQ=39)
                                                                                                                                                                          aaTGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGCATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAACCCAGTCAGTTGCTCTGGTTTGGt                                                             >  2:101231/1‑251 (MQ=255)
                                                                                                                                                                           aTGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGTCGCTGGCGCGTCTTATC‑AGGCCTACGTTAATTCTGCAATATATGGAAGCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTGCGCTATTAAAACCAGTCAGTTGCGCTGTTTTGGTc                                                            <  2:568926/251‑1 (MQ=255)
                                                                                                                                                                            ggCCGGTTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATTCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTTTTAAAACCAGTCAGTTGCTCTGGGGTGGTCa                                                           <  2:386087/250‑1 (MQ=255)
                                                                                                                                                                             gCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCa                                                           <  2:182571/250‑1 (MQ=255)
                                                                                                                                                                                    aaCAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGAttt                                                   >  2:439937/1‑251 (MQ=255)
                                                                                                                                                                                           cGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCaata                                             >  1:588982/1‑250 (MQ=255)
                                                                                                                                                                                            ggC‑ATTCAGAGCCTGATGCGTCGCTGGGGGGTCTTATC‑AGGCCTCGGTGCATTCTGCAATATATTGAATCTGCATGCTTTTGTAGCCCGGATAAGGCGTTAACGCAGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCCGTTCGCCATTGCTGGAGGAAGCTTAAACAAAGTAGTAACCTTAGCTATTAACACGAGTCTGTTGCTCTGGTTTGGTCATCCGATCTTCtataa                                            >  2:30322/1‑250 (MQ=255)
                                                                                                                                                                                                   cAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACg                                     >  1:444429/1‑251 (MQ=255)
                                                                                                                                                                                                        ccTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCATTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtca                                <  1:499542/251‑1 (MQ=255)
                                                                                                                                                                                                           gaggCGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCTTTAACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGGTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACtcatca                             <  2:581382/248‑1 (MQ=255)
                                                                                                                                                                                                               cGACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGAcc                         >  2:563882/1‑251 (MQ=255)
                                                                                                                                                                                                                gACGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACcg                        >  1:121295/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGCTGGCGCGTCTTATC‑CGGCCTACGTGAATTACGCCATATATTGCATCTGCATGCTTTTGTAGGCCGGATACGGCGTTCACGCCGCATCCGGTATTGACTGCAAACTTAACGCTGCTCGTCGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGACGTCACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGCCTCATCAGACagcg                      >  2:30890/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACcgcg                      <  1:439937/251‑1 (MQ=255)
                                                                                                                                                                                                                   gCTGGCGCGTCTTATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAACACGACTCATCAGACcgcat                     >  1:371235/1‑249 (MQ=255)
                                                                                                                                                                                                                              ttATC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGc          >  1:23864/1‑251 (MQ=255)
                                                                                                                                                                                                                                aTC‑AGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATACTAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt        >  2:400483/1‑251 (MQ=255)
                                                                                                                                                                                                                                    aGGCCTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCCGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGt        >  1:343904/1‑248 (MQ=255)
                                                                                                                                                                                                                                       ccTACGTGAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAAGCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATAAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGc  >  2:325508/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
ACTGGGATGGAGGTCACTGGCGCGACCACGGCTGGTGGAAACAACATTATGAATGGCGAGGCAATCGCTGGCACCTACACGGACCGCCGCCACCGCCGCGCCACCATAAGAAAGCTCCTCATGATCATCACGGCGGTCATGGTCCAGGCAAACATCACCGCTAAATGACAAATGCCGGGTAACAATCCGGC‑ATTCAGCGCCTGATGCGACGCTGGCGCGTCTTATC‑AGGCCTACGTTAATTCTGCAATATATTGAATCTGCATGCTTTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCATTGACTGCAAACTTAACGCTGCTCGTAGCGTTTAAACACCAGTTCGCCATTGCTGGAGGAATCTTCATCAAAGAAGTAACCTTCGCTATTAAAACCAGTCAGTTGCTCTGGTTTGGTCAGCCGATTTTCAATAATGAAACGACTCATCAGACCGCGTGCTTTCTTAGCGTAGAAGC  >  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: ‑