Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,678 T→C 9.5% intergenic (+287/+428) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,6780TC 9.5% 157.1 / 3.4 61intergenic (+287/+428)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (3/3);  total (32/29)
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.96e-01

CGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCA‑AGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGC  >  NC_000913/2304397‑2304952
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                   
cGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCGTCACCCGGCATCTAACcccc                                                                                                                                                                                                                                                                              >  1:262204/1‑289 (MQ=255)
         cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTACGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGCCGGATAAGGCGGTTACGCCGCCTCCGACATCTAACGCCCAAGCCCGtt                                                                                                                                                                                                                                                                     >  1:543744/1‑289 (MQ=255)
              ccccGAAGGGCAGATAAGGGGCGCAGGCCGCATCGGACGTCCGCGGCCGTAGCCGGTTGCCTGATGCGACGCTGGCGCTTCGTTCAAGGGCTAAACCGCTGTGAATGGCGGCACCCCGTAGGTCGGACAAGGCGGTCAAGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCGGATGCGACGCTGGGGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCt                                                                                                                                                                                                                                                                 <  2:313274/288‑1 (MQ=11)
                 cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTCCGCCGCATCCGACATCTCACGCCCAAGCCGGTTGCCTGAt                                                                                                                                                                                                                                                              >  1:631835/1‑288 (MQ=255)
                 cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGAt                                                                                                                                                                                                                                                              >  1:619200/1‑288 (MQ=255)
                                gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTACCCCGCTGTGAAGCGCTCCACCCCGTAGATCAGACAAGGCGGTCACGCCGCATCCGCCATCCAACGCCCGCGCGAGTTGCCTGGTCACACGCTggaggc                                                                                                                                                                                                                                                >  2:224076/1‑283 (MQ=11)
                                     gacggggcatccgaacgctagcgcgcaggCCGGTTGCCGGATGAGGAGCTGGGGCGTGTTGTCAGGTCTAATCCGGTGTGAAGGGCGGAACCCCGTAGGTCAGACAGGGCGGGCACGCCGAGTCCGCCATCCAGCGCCCGAGCC‑GGTGGCCTGATGCGACGCTGGCGCGGCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  2:699504/260‑1 (MQ=11)
                                     gaagacgcaaccgagtcccaccgcgctggCCGGTTGCCTGATGGGACGCTGGCGGTTCGTAGCAGGGCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGGCAAGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  2:457013/260‑1 (MQ=11)
                                     cacggctcgaccgcagcgccacgccacggccggttgccggttGCGCCGCGGGCGCGTCTTTTCAGGCTTAGCCCGCGGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTTCGGCGGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  2:543744/260‑1 (MQ=11)
                                     cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  1:111810/288‑1 (MQ=40)
                                     cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  1:192430/288‑1 (MQ=40)
                                     cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  1:562736/288‑1 (MQ=40)
                                       cgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                                                                                                                                                                       <  1:276022/289‑1 (MQ=40)
                                        gcgaagtccgtgatacgacccccGAGGCGGGTGCGGGAGGCGAAGCTGCCGAGGCTGATCAGGCCTACACAGCTGGGAAGTGCGGCACCCCGGAGGGCGGACAAGGCGGTCACGCCGCATCCGACATCCGACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                                                                                                                                                                       <  2:545226/268‑1 (MQ=11)
                                              tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCCACATCTAACGCCAAAGCCGGTTGCCTGATGCGAAGCTGGCGCGTCTTATCAGGCCTac                                                                                                                                                                                                                                 >  1:561088/1‑288 (MQ=255)
                                              ccggactaccAACGCCCGACCGGGTTGCCTGAGGCAACGATGTCCCGTCTTATAAGGCCCAAACGGCTGGGAAGTGCGCCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                <  2:312272/281‑1 (MQ=11)
                                                                                  acgctgacggggctgagagggCCTAAACCGCTGTGAAGTGCGGCACCCCGCAGGTCAGACAAGGCGGTCAGGCGGAGTCCGACATCCGACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                                                                                                                                                                                             <  2:146619/269‑1 (MQ=11)
                                                                                                  tCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCCACACCGCTGTGCAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCCACATCTAACGCCCAAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCCTCTCAGACCTACACCGCTGCGAAGCGCCCGACCCCGTAGGCTGGCTAAGCCGGCTCCGCCGCAt                                                                                                                                                                             >  2:178056/1‑288 (MQ=11)
                                                                                                      gCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCgcc                                                                                                                                                                        >  1:183111/1‑287 (MQ=38)
                                                                                                           cacCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCCCAACCGACATTCa                                                                                                                                                                   >  1:461371/1‑289 (MQ=11)
                                                                                                           accgcccgtagaggacGGTACCGCGTGGTTCGGACAAGGCGGCCACGCCGGATACGAGACACAAAGCCCGGGCC‑GGGTGCCTGATGCGACGCTGGGACGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACAGCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                                                                                                                                                                   <  2:212980/274‑1 (MQ=11)
                                                                                                                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCGGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                                                                                                                                               >  1:624709/1‑288 (MQ=255)
                                                                                                                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                                                                                                                                               >  1:145468/1‑288 (MQ=40)
                                                                                                                        tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCgg                                                                                                                                                      >  1:569715/1‑289 (MQ=255)
                                                                                                                        tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCAACCGACATCCAACGTCCGAGCCgg                                                                                                                                                      >  1:606910/1‑289 (MQ=255)
                                                                                                                           ggCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGtt                                                                                                                                                    >  2:174808/1‑288 (MQ=40)
                                                                                                                               ccccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTg                                                                                                                                               >  1:205148/1‑289 (MQ=39)
                                                                                                                                  cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCATCGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCTCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                >  1:553132/1‑205 (MQ=255)
                                                                                                                                  cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGAt                                                                                                                                             >  1:46373/1‑288 (MQ=37)
                                                                                                                                                 gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                                                                                                                                         >  1:70889/1‑181 (MQ=33)
                                                                                                                                                     tgaCGCGGCATCCGACAGCCCACGCCCGATCA‑GGTTGCCAGATGCGAGTCTGGGGGGTATTATCAGGCCTACACCGCGGTAAAGTGCTACACCCCGTAGGTCGGAAAAGGCGGTTACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                         <  2:619200/287‑1 (MQ=11)
                                                                                                                                                      taCGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          <  2:80743/174‑1 (MQ=14)
                                                                                                                                                      caggccgagtCCGACACCAAACGCCCGGGCG‑GGTTGCCGGATGGGACGCTGGGGAGGCTTATCAGGCCTACACTGCTGTGAAGGCCTCCCCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTGGCCTGAGGGGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                         <  2:721856/279‑1 (MQ=11)
                                                                                                                                                      cacgcgtcatcagagagtcAACGGCCGGGTC‑TGTGGCCTTATGCGAGGCTGGCGCGTCTGATCGGGGCTAAACCGCGGTGAGGTGCTCCACCCCGTAGGTCGGACAAGGCGGTTACGCCGCATCCGAGATCCGACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                         <  2:537944/270‑1 (MQ=11)
                                                                                                                                                      cACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                          >  1:80743/1‑175 (MQ=32)
                                                                                                                                                      cACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                                                        >  1:32912/1‑289 (MQ=34)
                                                                                                                                                             cATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                                                                                 <  1:174808/289‑1 (MQ=37)
                                                                                                                                                               tCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                               <  1:386668/289‑1 (MQ=37)
                                                                                                                                                                ccGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTg                                                                                                                                                                                                                   <  2:297488/188‑1 (MQ=31)
                                                                                                                                                                   acctccgacaccgaaGCC‑TGGTGCCTGAGGGGAGGCTGGAGCGTCTTATCAGGCCTAAACCGCGGTGAAGTGCTCCACCCCGTAGGGCGGACAAGGCGGTCACGCCGCATCCGACATCTAGCGCCCAAGCCGGTTGCCTGATGGGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCg                                                                                                            <  2:474778/274‑1 (MQ=11)
                                                                                                                                                                       ccAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                                                                                                                                                                                                     >  2:141667/1‑195 (MQ=31)
                                                                                                                                                                        cgaCGGCCGAGGCGGGTAGCCTGATGCGACCATGGCGGGTGGTATCAGTGCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGAGAAGGCGGGTACGCCGCATCCGAGAGCTAACGCCCAAGCCGGTTGCCTGATGCGACGCGGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtgtg                                                                                                       <  2:372030/287‑1 (MQ=14)
                                                                                                                                                                         aaCGCCCGTGCC‑GGTTGCCTGATGCGACGTTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGTTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                     <  1:469768/289‑1 (MQ=255)
                                                                                                                                                                         aaCGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                     <  1:75605/289‑1 (MQ=255)
                                                                                                                                                                          aCGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCg                                                                                                                                                                                                         >  1:297488/1‑188 (MQ=31)
                                                                                                                                                                                                    aCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCa                                                                                             >  1:540442/1‑271 (MQ=39)
                                                                                                                                                                                                    aCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCa                                                                                             <  2:540442/271‑1 (MQ=39)
                                                                                                                                                                                                                            tcaccgcgtcgAAGTGTGCCCCCCCGTAGTTCGGATGAGGCGTTTACGCCGCCTCCGACAGCTACCGCCGAAGCGTGTTGCCTGAGGCGACGCTGGCGCGTTTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTAACGCCGCATCCGACATCCAACGTCCGCGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCt                                                   <  2:145468/279‑1 (MQ=14)
                                                                                                                                                                                                                                cccTGTGAAGTGCTCCACCCCGGCGGTCGGATAAGGCGGTTCCGCCGCATCCGACATCTAACGCCCAGGCCGGTTGCCTGCTGCGACGCTGTCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                               <  1:696771/287‑1 (MQ=255)
                                                                                                                                                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTCCGCCGCATCCGACATCTAACGc                                              >  1:101162/1‑288 (MQ=255)
                                                                                                                                                                                                                                  cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                              >  1:81727/1‑288 (MQ=255)
                                                                                                                                                                                                                                   ttggaatgcggcccccccGAAGGTGGGTTAAGGCGGTTAGGCCGCATCCGACATCTAACTCCCACGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATAAGGCATAAACCGCGGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGcc                                             <  2:205148/274‑1 (MQ=14)
                                                                                                                                                                                                                                       aaGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCA‑Ag                                        >  1:519583/1‑288 (MQ=255)
                                                                                                                                                                                                                                            cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGCTAACTGCTCCACCCCTTTGGTCGGAACAGCCGGGTAGGCCGCATCCGAACTCCACCGGCCA‑AGCCGGtt                                  >  2:17059/1‑289 (MQ=11)
                                                                                                                                                                                                                                                         gTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                                                                                >  1:616798/1‑215 (MQ=255)
                                                                                                                                                                                                                                                         gTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                                                                                <  2:616798/215‑1 (MQ=255)
                                                                                                                                                                                                                                                         gTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGACGCATCCGACATCTAACGCACA‑AGCCGGTCGCCTGATGCGAc                      >  1:491917/1‑288 (MQ=255)
                                                                                                                                                                                                                                                           cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCTCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGCCCCGTCTTATCAGGCCCCCACCGCTGTTAAGTGCTCACCCCCATCGGTCGGATCAGGCGTTTACGCTCACTCCCACTTC‑CACGCCCACAGCCGGCTGCCTTATGTGACGCt                   >  2:99157/1‑289 (MQ=14)
                                                                                                                                                                                                                                                                            ccgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCA‑AGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAgg   >  1:728864/1‑288 (MQ=40)
                                                                                                                                                                                                                                                                            ccgcATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCA‑AGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAgca  <  1:17059/289‑3 (MQ=18)
                                                                                                                                                                                                                                                                             cgcATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGCACAGGCGGTTACGCCGCCTCCGACATCTAACGCCCA‑CGCCGGTTGCCTGATGCCACGATGgccccgtttctccggc  >  2:723342/1‑274 (MQ=14)
                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                   
CGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCC‑GGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCA‑AGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGC  >  NC_000913/2304397‑2304952

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑