Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913949,1180CT27.3% 108.9 / 18.2 66E151K (GAA→AAA) ycaNputative DNA‑binding transcriptional regulator YcaN
Reads supporting (aligned to +/- strand):  ref base C (15/33);  new base T (18/0);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 3.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACGGCCTCCAGTTCAGCATCAACATCATCCAGAATGATATTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCA  >  NC_000913/948879‑949368
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             
aaCGGCCTCCAGTGCAGCATCAAAATCAGCCTGAATGATATTTCCTGCGACAGCAATTGCCAGCTCCCTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCTGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACAGAATTTCGCTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAAc                                                                                                                                                                                                                                                    <  2:617097/251‑1 (MQ=255)
          aGTTCAGCATCAACATCATCCAGAATGATATTTCCTGCGACAGCAATTTCCAGCGCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGAAGATGCAACTAAGTCGAAc                                                                                                                                                                                                                                          <  2:596754/251‑1 (MQ=255)
                   tCAACATCATCCAGAATGATATTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCCCCCTGATCATATCTTTTTTGACGGTGCAAATAAGTCGAAAGCCCGCGTc                                                                                                                                                                                                                                 >  1:302182/1‑251 (MQ=255)
                    cAACATCATCCCGCATTATAGTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGAGAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTTCTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCCAACGCCAGCGTc                                                                                                                                                                                                                                 <  2:474514/250‑1 (MQ=255)
                      aCATCATCCAGAATGATATTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCTACGATGCAACTAAGTCGAACGCCAGCGTC‑aa                                                                                                                                                                                                                              >  1:171541/1‑250 (MQ=255)
                             ccAGAATGATATTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCCCGCTGGTCATATCTTTTTTGACGATGCAACTAAGACGAAAGCCAGCGTC‑AAAACCCTGt                                                                                                                                                                                                                      >  1:322000/1‑251 (MQ=255)
                               aGAATGATATTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTg                                                                                                                                                                                                                    <  1:438541/251‑1 (MQ=255)
                                            cTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGGGCG‑CCACGCTGAGCATATCTTTTTCGACGATGCCAATAAATCGAACGACCGCGGC‑CAAGACCTGTTTGAAGAAATCGAc                                                                                                                                                                                                        >  1:444415/1‑250 (MQ=255)
                                            cTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTTGACGATGCAACTAAAACGAACGCCAGCGTCCAAACCCCT‑TTTGACGAAATCGAc                                                                                                                                                                                                        >  1:140039/1‑250 (MQ=255)
                                                   aGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTTGCCGATGCAACTAAGTCGAACGCCAGCGTC‑AAAACCCCGTTGGGCGAAATCGACAagagagg                                                                                                                                                                                                >  1:249447/1‑250 (MQ=255)
                                                      aaTTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAgtcgt                                                                                                                                                                                             <  2:15192/251‑1 (MQ=255)
                                                      aaTTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCA‑CCCCGCTGGTTATATCTTTTTTGACGATGGAACTAAATCGAAAGCCACCGTC‑CAAGACCTGTTGGGCGAAATCGACAagaggatggt                                                                                                                                                                                             >  1:203463/1‑245 (MQ=255)
                                                          tCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGGTCG‑CCACGCCGATCATATCTTTTTTGCCGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTgggcggt                                                                                                                                                                                                              >  1:350320/1‑225 (MQ=255)
                                                          tCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGAGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTAATCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCGGCCCCG‑TGGTCATATCTTTTTTGTCGACACAACTAAGGCGAACACCAGCGTC‑CAAGCCCTGTTGGACGCTATCGACAagagggtgggc                                                                                                                                                                                            >  1:194816/1‑241 (MQ=255)
                                                          tCCAGCTCCTNTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGat                                                                                                                                                                                                              <  2:350320/230‑1 (MQ=255)
                                                                        aaaTTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGGTCATATCTTTTTTGACGATGCAACTAAGTCGAACGCCAGCGGC‑AAAGCCCCGTTGGACGATATCGACAAGAGAGTGGGCAGTTGTGGGCTCGACtt                                                                                                                                                                           >  1:287796/1‑251 (MQ=255)
                                                                        aaaTTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGAACG‑CCACGCTGATCCTATCTTTTTTGACGAGGCCAATAAGACGACAGCCCGCGGC‑AAAGCCCTGTTGGACGATCTCNACAAGAGAGTCGTCCGTTGTGAGCTCGACtt                                                                                                                                                                           >  2:361735/1‑251 (MQ=255)
                                                                              ccAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGCACG‑CCACGCTGATCATATCTTTTTTTAAGATTCAACTAAATCGAACGCCAGCGTC‑AAAGCACTCTTGGACGATATCGACAAGAGAAGCGTCCGGTGTGAGc                                                                                                                                                                                  >  1:593680/1‑238 (MQ=255)
                                                                              ccAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTGCGACGATGCAACGAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGc                                                                                                                                                                                  <  2:593680/238‑1 (MQ=255)
                                                                               cAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCCCGCTGATAATATCTTTTTTGACGATGCAAATAAGCCGAACGCCCACGGC‑AAAGCCCCGTTGGACGAAATTGAAAAAGGAGGCGGCAGTTGTGAGGTCGGCTTTAatatt                                                                                                                                                                    >  1:679853/1‑250 (MQ=255)
                                                                                     aaaGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGAGCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACGCTGGCGTTGCTGCTACACATAATTTCACGGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGGCGTCAGTTGTGAGCTCGACTTTAATATCAGGata                                                                                                                                                              <  2:32679/251‑1 (MQ=255)
                                                                                       aGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATAtt                                                                                                                                                            <  1:600265/251‑1 (MQ=255)
                                                                                              ccGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAAGCGAACGCCAGCGTC‑ACAGACCTGTTGGACGATATCGACAAAAGAGGCGTCAGTATTGAGCTCGACTTTAATATCAGGAAATTAACgtgt                                                                                                                                                     >  1:670581/1‑251 (MQ=255)
                                                                                                gCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGa                                                                                                                                                   <  1:141807/251‑1 (MQ=255)
                                                                                                    ggATAAAGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACAGAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAaacc                                                                                                                                               <  2:708592/251‑1 (MQ=255)
                                                                                                      aTAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGACCGATCC‑CCACGCTGATCCTATCTTTTTCGACGATGCCACTAAGTCGAACGCCAGCGGC‑AAAAACCCGTTGGACGATATCGACAAGAGAGTCGCCAGATGTGAGCTCGACTTTAATATCAAGATATTCACGtgtggaaaca                                                                                                                                              >  1:32505/1‑244 (MQ=255)
                                                                                                       tAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGATTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGGACG‑CCACGCCGATCATATCTTTTTTGACGATGCAACAAAGTCGAACGCCAGCGTC‑AAAGACCAGTTGGACGAAATAGACAAGAGAAGAGTCCGTTGGGCGGTCGACTTTAATATTAGGATATTTACGAGTGAaaccaa‑c                                                                                                                                           >  1:396199/1‑251 (MQ=255)
                                                                                                        aaCGAAAAACAACACATTGATGATTAAGCAGNTCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAaaccaa‑cc                                                                                                                                          <  2:331038/251‑1 (MQ=255)
                                                                                                                    cacaTTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGCTCC‑CCACGCCGATCATATCTTTTTTGACGATGCAACTAAATCGAACGCCAGCGGC‑AAAGACCCGTTTGACGCTATCGACCAAAGAGTAGTCAGTTGTGAGCTCGAATTTAATATAAGGATATTTACGGGTGAAACCAA‑ACAGTAAAGAAAtt                                                                                                                              >  1:263506/1‑251 (MQ=255)
                                                                                                                         tgatgaTTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTaaaa                                                                                                                          <  2:679853/250‑1 (MQ=255)
                                                                                                                          gatgatTAAGCAGATCGTGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACGGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACTCCAGCGTC‑AAAGCCCTGTTGGACGATTTCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTaaaaa                                                                                                                         <  2:656563/250‑1 (MQ=255)
                                                                                                                                tAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACCCTGATCATATCTTTTTTGACGATGCAACTAAGTCGAACGCCAGCGTC‑CAAGCCCCGTTGGACGATATCGACAAGAGAGGCGTCAGTTGTGAGCTCGGCTTTTAAATCAGGATATTCACGTGTGAAAACAA‑CCCATAAAAACAATAAAAATATACGt                                                                                                                  >  1:372330/1‑251 (MQ=255)
                                                                                                                                 aaGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGGC‑AAAACCCTGTTGGACGATATCGACAAGAAGGGCGTCAGGTGTGGGCTCGACTTTAATATCAGGATATTTACGGGTGAAACCAA‑CCAGTAAAAACATTAAAAATATAAGTg                                                                                                                 >  1:217136/1‑251 (MQ=255)
                                                                                                                                      gATCATGAGGATGGCGGGGTTTGCCATATCTTTCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGTCGGGCCGATCG‑CCACGCTGAACATATCTTTTTCGACGAt                                                                                                                                                                                                                                                         <  2:622079/112‑1 (MQ=255)
                                                                                                                                      gATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGAACC‑CCACGCCGAATATATCTTTTTTGACGAt                                                                                                                                                                                                                                                         >  1:622079/1‑112 (MQ=255)
                                                                                                                                               gATGGCGGGGTTTGCCATATCTTGCAAAGTACTCGGGCGTTGCTGCTTCACATAATTTCCCTGGCGGGCCGATCG‑CCACGCTGATCATATCTTGTGCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGggcg                                                                                                  <  2:670581/251‑1 (MQ=255)
                                                                                                                                                   gCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCAt                                                                                              <  2:249447/251‑1 (MQ=255)
                                                                                                                                                        gTTTGCCATATCTTGCCAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGAGCC‑CCACGCTGTTCCTATCTTTTTTGAAGATGCAACTAAGACGGACCCCAACGGC‑AAAACCCTGTTTGACGGTAACGACAACAGAGTCGTCAGTTGGGAGACCGACTTTAATAACAGGATATTAACATGTGAAAACAA‑ACAGTAAAGAAATTAAAAAAATACGTGCAGC‑TCCCCGGGCGGCATTGAtt                                                                                         >  2:619947/1‑251 (MQ=255)
                                                                                                                                                           tGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGAtttt                                                                                       <  2:287796/250‑1 (MQ=255)
                                                                                                                                                              cATATCTTGCAAAGCACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCC‑CCACGCTGGGCAAATCTTTTTCGACCATGCAACTAAGTCGGACCCCAGCGGC‑AAAGCCACGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTTGGCCCTCCTTTTATATCAGGAAATTTCCGTGTGAAAAAAA‑CCCGTAAAGACATTAACAAAATACGAGCAGC‑TCACCGGGCGGCATTGATTTTGatg                                                                                   >  2:109102/1‑249 (MQ=255)
                                                                                                                                                               atatCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCCCTCG‑CCACGCTGCTCATAACTTTTTTGACGCTGCAAATAAGGAGAACGCCAGCGGC‑AAAGACCTGTTGGGCGATATCGCCAAAAGAGACGTCAGTTGTGGGGCCGTCTTTTATATCAGGATATTCCCGTGTGAAACCAAACCCG‑AAAAAAATTAAAAAAATACTNGCAGC‑TACCCCGGCGGGATTGAtttttggcg                                                                                  >  2:333369/1‑246 (MQ=255)
                                                                                                                                                                     tGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCTATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGGCGAACGCCATCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTTCAGC‑TACACGGGCGGCATTGATTTTGAGCGtt                                                                                <  2:508629/247‑1 (MQ=255)
                                                                                                                                                                     tGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCCATCC‑CCACGCTGATCCTATCTTTTTTGACCGTGCCACCAAGTCGAACGCCAGCGTCCAAAACCGTGT‑GGGCGATATTGCCAAGAGAGGCGTACGTTTTGAGCTCGGCTTTTATATATGGATATTCACGTGTGAAACCCA‑CCCATAAAAACATTAAAAAAATACGTGCCGCCTACC‑TGGCGGCATTGCTTTTTAGCGTTCCtt                                                                            >  2:240464/1‑250 (MQ=255)
                                                                                                                                                                      gcaaatgaCTCTGGCGTTGCTGCTACACATAATTCCACTGGCGGGACGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACTCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGATAGTCGTCAGTTGTGAGCGCGACTGTAATATCAGTATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTg                                                                            <  2:669831/243‑1 (MQ=255)
                                                                                                                                                                      gCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGt                                                                           <  1:228438/251‑1 (MQ=255)
                                                                                                                                                                          aGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTTACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGTGCGTTCCTGTCggg                                                                       >  1:379382/1‑250 (MQ=255)
                                                                                                                                                                             aCTCTGGCGTTGCTGCTACACATAATTTCACTGTCGGGCCGATCG‑CCACGCTGAACATATCGTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGgtgt                                                                     <  1:382565/250‑1 (MQ=255)
                                                                                                                                                                             aCTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTc                                                                    <  1:386965/251‑1 (MQ=255)
                                                                                                                                                                                  gtcGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCTCGCTGATCAGATCTTTTTCGACGATGCAACTAAGTCGAACGCCGGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGGCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCATTAAAGACAGTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCTTTCCTGTCGGTGTCAGACg                                                               <  2:23213/249‑1 (MQ=255)
                                                                                                                                                                                     gTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCCTATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTTGACGATATAGACAAAAGAGTCGTCAGTTGTGAGCTCGACCTTAATATCAGGGTATTCTCGTGTGAAAACAA‑CCAGTAAAAAAATTAAAAAAAAACGTGCAGCCAACCGGGGGGGCAT‑GATTTTTAGCGGTCCTTTCGGTGTCAgcagaa                                                             >  1:333677/1‑245 (MQ=255)
                                                                                                                                                                                     gTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGaaa                                                            <  1:136145/251‑1 (MQ=255)
                                                                                                                                                                                         ctgctACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTaaa                                                                                                                                    <  2:716741/176‑1 (MQ=255)
                                                                                                                                                                                           tctACTCATAATTGCACTGGGGGGCCGATCG‑GCACGCGGATCATACCTTTTCCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AACGTCCTGTTGGACAAGAGCGTCAAGAGAGTCGTCAGTTGTTAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATAAAAAATATACGTGCAGC‑TACCCGTGCGGCATTGATGTTGAGCGTTTCTGTCGGTGTCAGACGAAAAtcat                                                       <  2:549727/249‑1 (MQ=255)
                                                                                                                                                                                                    aatttcccgtgCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCACCTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATACCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACAGTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGACAATCATtcatctcat                                              <  2:395345/240‑1 (MQ=255)
                                                                                                                                                                                                    aaTTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATtcatctcat                                              <  1:32320/250‑1 (MQ=255)
                                                                                                                                                                                                             ggCGGGCCGATCG‑GCACGCTGATCATATCTTTTTCGACGATCGAACTAAGTCGGACGTCCGCGTC‑AAAGCCCTATTGGACGATATAGACAAGAGAGTCGCCAGTAGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGCAACCAA‑CCAATCCAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTACGTCGAAAATCATTCATCTCATCTacaatgat                                    >  2:101177/1‑244 (MQ=255)
                                                                                                                                                                                                                     gATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAAttt                            >  1:171914/1‑251 (MQ=255)
                                                                                                                                                                                                                        cG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTtcgt                         >  2:648407/1‑251 (MQ=255)
                                                                                                                                                                                                                           ccACGCTGATCATAGCTTTGGCGAAGATGCAACTAAGGCGAACGCCAGCGAT‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATGCGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCGGAATTtcgtcg                       <  2:444415/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCCCATCTAACATGATCTGAATTTCGTCGAAAGCTg                >  2:265871/1‑251 (MQ=255)
                                                                                                                                                                                                                                      atatCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGa                                                                                     >  2:324635/1‑178 (MQ=255)
                                                                                                                                                                                                                                      atatCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGa                                                                                     <  1:324635/178‑1 (MQ=255)
                                                                                                                                                                                                                                        atCTTTTGCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGgcgc           <  2:154715/250‑1 (MQ=255)
                                                                                                                                                                                                                                            ttttCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATc       <  2:171541/250‑1 (MQ=255)
                                                                                                                                                                                                                                             tttCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGt     >  2:126071/1‑251 (MQ=255)
                                                                                                                                                                                                                                              ttCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTAGTCGACAGCTGGGCGCAATCGtt    >  1:417162/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cgacgaTGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTc   >  2:508520/1‑250 (MQ=255)
                                                                                                                                                                                                                                                cgacgaTGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCCCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCa  <  1:265871/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cgacgaTGCAACTAAGGCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCa  <  2:263506/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                             
AACGGCCTCCAGTTCAGCATCAACATCATCCAGAATGATATTTCCTGCGACAGCAATTTCCAGCTCCTTTGCAAATTGCCAGTGGAAAGGTTTGCCGCTGGGATAACGAAAAACAACACATTGATGATTAAGCAGATCATGAGGATGGCGGGGTTTGCCATATCTTGCAAAGTACTCTGGCGTTGCTGCTACACATAATTTCACTGGCGGGCCGATCG‑CCACGCTGATCATATCTTTTTCGACGATGCAACTAAGTCGAACGCCAGCGTC‑AAAGCCCTGTTGGACGATATCGACAAGAGAGTCGTCAGTTGTGAGCTCGACTTTAATATCAGGATATTCACGTGTGAAACCAA‑CCAGTAAAGACATTAAAAATATACGTGCAGC‑TACCCGGGCGGCATTGATTTTGAGCGTTCCTGTCGGTGTCAGACGAAAATCATTCATCTCATCTAACATGATCTGAATTTCGTCGAAAGCTGGGCGCAATCGTTCA  >  NC_000913/948879‑949368

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