Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,563,778 (G)6→5 coding (148/258 nt) ytiA ← uncharacterized protein YtiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,563,7780G.90.9% 158.5 / 13.0 44coding (148/258 nt)ytiAuncharacterized protein YtiA
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base . (20/20);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

ACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAG  >  NC_000913/4563532‑4564019
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ccccacTTCCACAATTGTAAAAAGAACCACAAAAATGGAGATGCTTTTACTGAACGAAACACTTTTTGTCACCGTAATCTATATAAGGCAAATCAATAAGATATCGAGAAGTATACGGAGCGACCTCAATGCTTCCTTCACGTGTAGAAACCCCGTAATGTGCCATGCCGCGCCAGACTCCATTCTGCTGGCAATAAATTGGCTCAGTCAGCCCAGGGATAACTCATTGTATTCCCTCGCCGGGGGACAgg                                                                                                                                                                                                                                              <  2:378616/248‑1 (MQ=255)
    accgCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCTCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAg                                                                                                                                                                                                                                          <  2:38078/248‑1 (MQ=255)
     cTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGc                                                                                                                                                                                                                                         >  1:319607/1‑251 (MQ=255)
            aaTGGCAAAAAGAATCACAATAATGGCGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGCTATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAGATTGACTCAGTCAGCCAAGGGATCACCCATTGTCTTTCATCGCCGGGGGACAGGTCAGCCAAATGc                                                                                                                                                                                                                                  >  2:488378/1‑251 (MQ=255)
             aTGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGc                                                                                                                                                                                                                                  >  2:328635/1‑251 (MQ=255)
                       gAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGCGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGc                                                                                                                                                                                                                         >  2:425075/1‑249 (MQ=255)
                                taatGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGt                                                                                                                                                                                                              >  1:20096/1‑251 (MQ=255)
                                 aatGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTa                                                                                                                                                                                                             <  1:145308/251‑1 (MQ=255)
                                         tGCGTTTCCTTAACGAAATACTTTTTGTCACGGTACTCCCTTTAAGGCAAATCCCTAAGATATCGAGCAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCCCCAGACACCATCCTGCTGGCACTAAATTGACTCAGTAAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACCGGTCAGCCAAATGCATACCCCGCAGCCGCCGCGTAAATCCCgg                                                                                                                                                                                                     >  2:93885/1‑251 (MQ=255)
                                               tACTTAACGAAATACTTTTTGTCACGGTCATCCCTTTAAGGCAAATCAATACGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAat                                                                                                                                                                                                >  1:270393/1‑251 (MQ=255)
                                                aCTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAata                                                                                                                                                                                               <  1:297203/251‑1 (MQ=255)
                                                 attAACGAAATTCTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                                                                                                                                                                             <  1:47881/250‑1 (MQ=255)
                                                    aaCGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGCCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATAGCt                                                                                                                                                                                          >  2:655109/1‑251 (MQ=255)
                                                              tttttGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaa                                                                                                                                                                                 >  1:146307/1‑250 (MQ=255)
                                                               ttttGTCACGGTAATCCCTTTAAGGCAAATCAATAAGTTATCGTGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGCCAGGTCAGCCACATGCATAACCCGAAGCCGCCGCGTAAATCACGGGGAATCTCGCTGACGCGaaa                                                                                                                                                                                 >  2:528132/1‑249 (MQ=255)
                                                                   gTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                                                                                                                                                                            <  1:328635/251‑1 (MQ=255)
                                                                       cGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGccc                                                                                                                                                                       >  1:522931/1‑251 (MQ=255)
                                                                           aaTCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctg                                                                                                                                                                   <  2:146307/251‑1 (MQ=255)
                                                                           aaTCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGct                                                                                                                                                                    <  1:2715/250‑1 (MQ=255)
                                                                                  ttAAGGCAAAACAATAAGATATCTAGAAGTATATGGAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCCGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATGTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCa                                                                                                                                                             <  2:11149/250‑1 (MQ=255)
                                                                                    aaGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAACTCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCGGCCTGCaccg                                                                                                                                                          >  2:624245/1‑248 (MQ=255)
                                                                                                  aGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACAATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTAAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCCTGCCCGCTGCCTGCAACGCTTTGCGGGCTTGt                                                                                                                                            >  1:520262/1‑251 (MQ=255)
                                                                                                   gATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTa                                                                                                                                           <  2:319607/251‑1 (MQ=255)
                                                                                                                ataCGGAGCGACCTTAATAAAGCCTTAACGCGGAGAAACCCCGTAATGTTCAATGCAGCGCCAGACACTATCCTGCTTGCAATAAAGTGGCTCAGTCAGCCCTGGGATCTCCCATTGTATTTCATCGCCGGGGGACAGGTCAGCCAAATGCATAACCCTCAGCCGCTGCGTAAGTCACGGGGAATATCGCTGACTCGAAAAGTCATGCCCGCTGCCTGCAATGCTTTGCGGGTTTGTATGActcagatcg                                                                                                                               <  2:44597/250‑9 (MQ=255)
                                                                                                                        cGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGa                                                                                                                       >  1:476171/1‑250 (MQ=255)
                                                                                                                         gACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaa                                                                                                                     <  1:482307/251‑1 (MQ=255)
                                                                                                                           ccTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaatc                                                                                                                   >  1:417435/1‑250 (MQ=255)
                                                                                                                                 ataaTCCTTAACGCGGAGAAACCCCGTAATGTGCAATCCCGCCCCAGTCACCTCCCTCCGTGCAATAAATTGACTAAGTCAGCCCAGGGAGCACCCATTTTAGTTCCTCGCCGGTGGACAGGTCAGCCAAATTCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAAGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGaaataaaaata                                                                                                             <  2:244246/249‑1 (MQ=255)
                                                                                                                                  gAATCCTTAACGCGGCGAAACCCCGTAATGAGCACTGCAGCGACAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCAAAATGCATAACCCGCAGCCGCCGAGAAAATCACGGGGAATATCGCTGACGCGAAAAGTAATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCACCCGTGGAGTGaagtacaacta                                                                                                             >  2:626653/1‑241 (MQ=255)
                                                                                                                                   aaTCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCGCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAaataaaaataaa                                                                                                           <  2:119091/251‑1 (MQ=255)
                                                                                                                                    aTCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAaataaaaataaa                                                                                                           >  1:217454/1‑250 (MQ=255)
                                                                                                                                          accggggAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTTCTGGCAAGAAATTGACTCAGCCAGCCCAGTGATCACCCATTGGATTGCCGCGCCGGTGTACAGGTCATCCAAATGCAGAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGTCCGTTGCCTGCAACGCATTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCtt                                                                                                    <  2:677108/246‑1 (MQ=255)
                                                                                                                                                  gaAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGc                                                                                            <  1:270446/251‑1 (MQ=255)
                                                                                                                                                       cccGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc                                                                                                                                                           <  2:524228/183‑1 (MQ=255)
                                                                                                                                                       cccGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc                                                                                                                                                           >  1:524228/1‑183 (MQ=255)
                                                                                                                                                          gTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCGCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTc                                                                                    <  2:164497/251‑1 (MQ=255)
                                                                                                                                                            aaTGTGCAATGCAGCGCCAGACCCCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCAGGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGGGACTTAAAAATAAAAACTCTTTCATCGGCGTGTCCTccc                                                                                  >  2:414866/1‑249 (MQ=255)
                                                                                                                                                                            ccAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAcc                                                                  <  1:691843/251‑1 (MQ=255)
                                                                                                                                                                            ccAGACACCATCCCTCGGGCAATAAATTGACTCAGGCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTAcc                                                                  <  2:192087/251‑1 (MQ=255)
                                                                                                                                                                                                        gACTCAGTCAGCCCAGGGATCCCCCATTTTTTTCCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTgcg                                      <  1:488378/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccAGGGATCCCCAATTGTATTCCCTCGCCGGGGGACAGTTCCGCCAACTGCATACCCCGCAGCCGCCGCGTAAATCCCGGGGAATATCGCGGACTGGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGTTTTTTATGACGCCGACCGTGGAGTGAAATAAAAAAAAAAACGCTTCCATCGGCGTGTCCTCACTCTTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATggg                          <  2:20096/251‑1 (MQ=255)
                                                                                                                                                                                                                         gATCACCCATTGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTTCGCCAGCAAATTTCGCATCAGGGTGGGTAttt                     >  1:157816/1‑251 (MQ=255)
                                                                                                                                                                                                                              cccATTGTATTTACTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGACGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGCCATGCCCGCTGCCTGCAACGATTTGCGGGTTTGTATGACGACGACCGTGGAGTGAAATAACAATAACAAATCTTTCATCGGCGTGTCCTCCATCGTTGACCAATTAACCCCGCGCCAAATGAGCCCGAAAATTGCGCATCAGGCTGGGGATTTAACgg                >  1:361034/1‑251 (MQ=255)
                                                                                                                                                                                                                                  ttGTATTTCCTCGCCGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGAGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGGTGACCAATTTCCGCCGCCCCCATTGCGCCAGCAAATTGCGCCTCAGGATGGGGTTTTACCGGCAtt            >  2:149/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ccgccggggaACAGGACAGCCAAAGCCAAAACCCGCAGCCGCCGCGGAAAGCACGGGGAATATCGCTGACGATAAAAGTCATGCCCGCTGCCTACAACGCTTGGGGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACCCTCTCACCGGCGTGTCCTCACTCGTTGACAAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTgag  <  2:217454/246‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                 
ACACACTGCCACAATGGCAAAAAGAATCACAATAATGGAGATGCGTTTACTTAACGAAATACTTTTTGTCACGGTAATCCCTTTAAGGCAAATCAATAAGATATCGAGAAGTATACGTAGCGACCTTAATGAATCCTTAACGCGGAGAAACCCCGTAATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTTGCGGGTTTGTATGACGCCGACCGTGGAGTGAAATAAAAATAAAAACTCTTTCATCGGCGTGTCCTCACTCGTTGACCAATTACCGCCGCGCCAATTGCGCCAGCAAATTGCGCATCAGGATGGGTATTTACCGGCATTCGCAGGTGAG  >  NC_000913/4563532‑4564019

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: