Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,796 T→C 10.3% intergenic (+405/+310) 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,7960TC10.3% 132.9 / 3.8 58intergenic (+405/+310)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base T (30/21);  new base C (3/4);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGG  >  NC_000913/2304509‑2305079
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                           
cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGt                                                                                                                                                                                                                                                                                             >  2:120757/1‑288 (MQ=40)
cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCCCCACCCCGTAGGTTGGATCAGGCGGTTACGCCGCATCCGACATCCAACGt                                                                                                                                                                                                                                                                                             >  2:401148/1‑288 (MQ=255)
             cACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGCTGCGACGC‑TGGCGCGTCTTATCAGGTCTACACCCATGTGAAGTGCTCCACCCCGCAGGTTGGATCAGGCGGTTACGTCGCCTCCGCCCTCCCACGTCCGAGCCGGTTGct                                                                                                                                                                                                                                                                               >  2:176090/1‑288 (MQ=14)
             cACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                                                                                                                                                                                                                                                                                 >  2:620052/1‑287 (MQ=39)
             cACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTg                                                                                                                                                                                                                                                                                 >  2:272119/1‑287 (MQ=39)
                  cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATCAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGCCCTACGCCGCTGTCAAGTGCCCCCCCCCGTAGGTTTGCTTAGGACGCTCCGCCCGACCCGACATCGAAGGCCCCAGCGCGTTGcatcat                                                                                                                                                                                                                                                                           >  1:385497/1‑283 (MQ=11)
                            aCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGgc                                                                                                                                                                                                                                                                >  2:58062/1‑289 (MQ=37)
                                      taCGCCGCATCCGACATCCACCTCCCGAGCTGTTTGCCTGATGCGTCGCTGGCCCGTCTTACAAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       <  2:749612/287‑1 (MQ=255)
                                        gcgcgcaggcggcgtCCAACCGCCGGGCCGGTGGCCGGAGGGGAACGTGGCGCCTCTTGTCGGGCCTACAGCGCGTGGAAGGGCGCCACCCCGTGGGGCGGAGAAGGCGGTTAGGCAGCATCCGAGATCAACCGCCCGAGCGGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCGGGCCTACACCGCCGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                                                                                                                                                                                    <  1:527666/275‑1 (MQ=11)
                                             cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGTTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGGAGGTTGGATAAGGCGGTTACGCCGCATCCGACACCCAACGTCCGAGCCGCTTGCCTGATGCGACGCCGGCGCGTCTTATCAGGCCt                                                                                                                                                                                                                                                >  2:425297/1‑288 (MQ=255)
                                               tccgacacgccccccccGAGCCGGTTGGCGGATGCGAGGGTGGGGCGTCTTACCCGGCCTACACGGCTGTGAAGGGCTCCACCCCTTGGTTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGGTGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                             <  1:120757/275‑1 (MQ=11)
                                               tcagactgacAACGCCGGGGCGGGTGGCCTGCTGGGACGCTGGCGTGCTTTAGCAGGCCTACGCCGCTGAGATGTGCGCCCCCCCGTAGGTCGGAGAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCGGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                             <  1:58062/280‑1 (MQ=11)
                                                                       gtGCCTGGCGCGCCGCTGGCGCGTGTTCGCAGGCCTACACCGCTGTGCAGTCGTCCACCCCGGGGGGCGGACAAGGCGGTTACGCCGCATCGAACAGCCAACGCCCGAGCCGGTGGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCGGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCg                                                                                                                                                                                                                      <  1:404949/287‑1 (MQ=11)
                                                                                 tcatgtgtgtacaatcgcaatggCATACACCGCGGTGAAGTGCTCCACCCCGTATGTCGGTTAGGGCGTTTCCGCCGCATCCGACATCTACCGCCCAAGCCGGTTGCCTGTTGCG‑CGCTTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGTCGGTTACGCCTCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGAt                                                                                                                                                                                                            <  1:759934/267‑1 (MQ=11)
                                                                                          cgcgTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTAAACCGCTGTGAAGTACTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGAATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCTCCACCCCGTAGGTCGGCTAAGGCGGtt                                                                                                                                                                                                   >  2:501807/1‑288 (MQ=255)
                                                                                           gcgTCGTAGCTGTCCTACACCCCTGTGAAGTGCGCCGCCCCGTAGGTCAGATAATGCGGTCAAGCCGCAGCCGCCATCGACCGCCCAGGCCGGTTGCCTGATGCGCCGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTa                                                                                                                                                                                                  <  1:272119/288‑1 (MQ=11)
                                                                                                       gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                                                                                                 <  1:298713/213‑1 (MQ=21)
                                                                                                       gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                                                                                                 >  2:298713/1‑213 (MQ=34)
                                                                                                        cctcaacCGCTGTCAAGCGCCCCCCCCCGTAGGTCGGAGAAGGGGGTTACGCCGCATCCGACATCTACCGCCCTAGCCTGTTGCCTGATGCAACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCGGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                                                                                                                                       <  2:320371/281‑1 (MQ=14)
                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGCTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       >  1:744570/1‑213 (MQ=21)
                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       <  2:744570/213‑1 (MQ=34)
                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGCTa                                                                                                                                                                                                  >  2:759934/1‑266 (MQ=16)
                                                                                                                 cTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTCAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGCCTTATCAGGCCTACACCTCTGTGAAGTCATCCACCCCGTAGTTCGGCTGAGGAGGTTATCGCGCATCAGCCATCTAACGc                                                                                                                                                                            >  2:690814/1‑288 (MQ=14)
                                                                                                                    tgAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGCCGGATAAGGCGGTTACGGCGCATCCGACATCGAACGCCCaa                                                                                                                                                                        >  2:813264/1‑289 (MQ=255)
                                                                                                                                   ctacggTCGAATAAGGCGGTACAGCCGCATCCCACATCTAACGCCCCGGCCGGTTGCCTGATGCGACGC‑TTGCGTGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                                         <  2:766784/285‑1 (MQ=14)
                                                                                                                                   cGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       >  2:425759/1‑195 (MQ=34)
                                                                                                                                           ggATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCCACGCTg                                                                                                                                                 >  2:755942/1‑289 (MQ=255)
                                                                                                                                                               agtcGACATCAAGAGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGGCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTTGGAAAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                                                                                              <  2:550915/285‑1 (MQ=11)
                                                                                                                                                                               ccAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                                                                                                                                                                                            <  1:283360/162‑1 (MQ=40)
                                                                                                                                                                               ccAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTacac                                                                                                                                                                                                                                            >  2:283360/1‑162 (MQ=40)
                                                                                                                                                                                  aGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTcc                                                                                                           >  2:669561/1‑288 (MQ=38)
                                                                                                                                                                                        tggccgattGCGAGCC‑TGCCGCGTCTGATCAGTCCTAACCCGATGTTACGTACTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCg                                                                                                     <  2:528420/280‑1 (MQ=11)
                                                                                                                                                                                         tGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCTCCCCCGCTGTCACGTGCCACACCCCGTAGGTCCGATCACGCGGTTACGCCGCCTCCGACATCTAACGCCCCAGCCGGTTGCCTGAAGCGACGCTTGCGCGTCTTATCCGGCCTAGACCCCtggtcatcgcccccccccgg                                                                                                    >  1:440940/1‑270 (MQ=11)
                                                                                                                                                                                                 gCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                      >  2:660465/1‑246 (MQ=39)
                                                                                                                                                                                                 gCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                      <  1:660465/246‑1 (MQ=39)
                                                                                                                                                                                                                tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGACGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCCGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCTCCACCCCCTAGGTTGGATAAGGCGGGTAcgc                                                                              >  2:53246/1‑288 (MQ=255)
                                                                                                                                                                                                                tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       >  2:375853/1‑119 (MQ=17)
                                                                                                                                                                                                                tCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       <  1:375853/119‑1 (MQ=17)
                                                                                                                                                                                                                     cgaGGCCAACGCCCCGGTGAAGGGCTCCCCCCCGTAGGTTGGATAAGGCGGTTTCGCGGCATCCGCCATCCAACGCCCGATCCTGTTGCCTGATGAGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATc                                                                        <  1:620052/287‑1 (MQ=11)
                                                                                                                                                                                                                             acacGGCTGGGAAGTGCTCCACCCCGGAGGTGGGTGAAGGCGGTTACGCCGGAGCCGAGGTCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAGCAGGCCTAGACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCAGCCAACATCt                                                                <  2:383603/289‑1 (MQ=11)
                                                                                                                                                                                                                             acacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCGACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTTCTCCACCCCGTAGGCTGGATAAGGCGGTTACTCCGCATCCCAAATc                                                                 >  1:45769/1‑288 (MQ=255)
                                                                                                                                                                                                                              cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                             >  1:318606/1‑115 (MQ=255)
                                                                                                                                                                                                                              cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                             <  2:318606/115‑1 (MQ=255)
                                                                                                                                                                                                                              cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                                         >  1:740844/1‑199 (MQ=255)
                                                                                                                                                                                                                              cacCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                                                                                         <  2:740844/199‑1 (MQ=255)
                                                                                                                                                                                                                                   cTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCTACCCCGTAGGCCGGATAAGGAGGTTACGCACCATCCGACCTCTAACGCCCACGCCGGTTGCCTCATGCGCCGCTGGCGCGACCTATCCGGCTCACACCGCTCTGACTTGCTTCACCCCGCCGGTTGGCTCAGGTGGTTCCACTGCATCCAACATCTacctc                                                           >  1:558221/1‑284 (MQ=14)
                                                                                                                                                                                                                                             cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTCCGCCCCCTCCACCATCCAACCCCCGACCCGGtt                                                >  1:454961/1‑289 (MQ=255)
                                                                                                                                                                                                                                                  ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGTCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTCAAGTGCTCCACCCCGTCGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACCCCCGATCCGGTTGccgg                                           >  2:376400/1‑287 (MQ=255)
                                                                                                                                                                                                                                                     cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTCGGTCGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGAt                                         >  2:121953/1‑288 (MQ=255)
                                                                                                                                                                                                                                                     cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATccga                                                                     <  1:27016/260‑3 (MQ=25)
                                                                                                                                                                                                                                                     cGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATccga                                                                     >  2:27016/1‑258 (MQ=25)
                                                                                                                                                                                                                                                       tAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGc                                       >  2:665927/1‑288 (MQ=39)
                                                                                                                                                                                                                                                             ggATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTg                                >  2:111279/1‑289 (MQ=37)
                                                                                                                                                                                                                                                               tgagtgCGGTTCCGCGCCAGCCGACATCCAACGTCCGAGCGGGTTGCCGAATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                               <  2:151234/283‑1 (MQ=255)
                                                                                                                                                                                                                                                                    gCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                                                                       <  1:423534/67‑1 (MQ=40)
                                                                                                                                                                                                                                                                    gCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTTCGCCGCATCCGGCATCTAACGCCCACGCCGGTTGCCTGATGCGACGCTGGCGCGTCCTATCCGGCCTCCACCGCTGTCAAGTGCACCACCCCGTAGGTTGGTTAAGACGGTTCCCCCGCACCCCACATCTCACCACCGCGCTCGTTCGCTGATTGTACGCTggttca                           >  1:151234/1‑283 (MQ=11)
                                                                                                                                                                                                                                                                         tccgcctcagccgcaagccAACGCCCGTGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCATGCCTACACCGCTGTGAAGGGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAGATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                     <  1:401148/271‑1 (MQ=11)
                                                                                                                                                                                                                                                                                cATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATAc             >  2:691917/1‑289 (MQ=37)
                                                                                                                                                                                                                                                                                      aCATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                                                             <  1:608915/59‑1 (MQ=35)
                                                                                                                                                                                                                                                                                            aaCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCgg  >  2:82269/1‑288 (MQ=40)
                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                           
CTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGC‑TGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCATACGCCACATCCGG  >  NC_000913/2304509‑2305079

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