Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,304,535	A→G	8.7%	intergenic (+144/+571)	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_000913	2,304,535	0	A	G	8.7%	140.5 / 3.8	61	intergenic (+144/+571)	eco/mqo	ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base G (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 = 6.39e-01

ATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGC  >  NC_000913/2304270‑2304809                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                   atgcgtgAAGGCTGCGTTACAAAGGGGAGCGGCCGTTCGGGGTGTATATGCCAGACAATTGAGCTGGGAAGAACCGCGTCTGGAAGGCGGAATGGAAAAGTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGGCAAGTCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATcc                                                                                                                                                                                                                                                             <  1:549248/283‑1 (MQ=18)            ctgcGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAAc                                                                                                                                                                                                                                                   >  1:326154/1‑288 (MQ=255)                  tACAACAGCAAGCTGCCGATCGTGGTTTATACGCCAGACAATGTAAATGTGAAGTACCGCGTGTGGAAGGCGGAAAAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAATTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGa                                                                                                                                                                                                                                             <  2:960176/288‑1 (MQ=255)                         gCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGtt                                                                                                                                                                                                                                      >  2:699463/1‑288 (MQ=255)                           aaGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGc                                                                                                                                                                                                                                    >  1:167167/1‑288 (MQ=255)                            aGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTg                                                                                                                                                                                                                                     >  2:211588/1‑286 (MQ=255)                            aGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGATGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCAGCAGCCGACATCCAACGCCCGAGCCGGTTg                                                                                                                                                                                                                                     >  2:1030933/1‑286 (MQ=255)                                 ccGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCGACCCACCGTAGGTCCGACACGTCGGTTCCGCCGCCCCCGCCATCCAACGCCCACGCCGGTCGCctgctga                                                                                                                                                                                                                             >  1:743539/1‑288 (MQ=21)                                  cGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGACGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGCAGGTCAGATAAGGCGGTCACGCCCCCTCCGCCCTCAAACCCCCGCGCCGGTTGCCTGCTGCg                                                                                                                                                                                                                            >  1:673603/1‑289 (MQ=21)                                        tggtgTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCAGTCACGCCGCATCCGACATCCAACGCCCGCACCGGTTCCCTCATACGACGCTg                                                                                                                                                                                                                      >  2:338613/1‑289 (MQ=255)                                         ggtgTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                                                     <  2:67541/289‑1 (MQ=255)                                          gtgtagagcgcAGGCAATGTGGATGGGCAGTACCGAGGATGGAAGGCGGCAGGGACACGTGACAACGCGGTAGTGCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACAGCAAACCCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                                                     <  1:840649/278‑1 (MQ=17)                                          gtgtATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                                                     <  2:167167/288‑1 (MQ=255)                                              ataCGCCAGACAAAGGAGATGGCAAGTGCCGGGACTGGAAGGCGGAAGAGACAAGTGACAACGCGGTAGTGCGCGAAGCTGCCGTGAAGTGCGGCACCCGCTAGGCCAGGCAAGGCGGTCACGCCGCATCCGAGAGCCAACGCCCGAGCGGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                                                <  1:814376/289‑1 (MQ=18)                                               taCGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGcgcg                                                                                                                                                                                                                 >  2:98973/1‑287 (MQ=255)                                                cggCCAGACAAAGTGGGTGTGAAGGACCGCGGCGGGCAGTCGGAGGAGAAAATGGACAACGGGGGAGGGCGCGAAGCTGCCGGAAGGTGCGGCACCCCGGAGGTCGGGAAAGGCGGGCACGCGGCAACCGACATCCAGCGCCCGAGCCGGGGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCTCCCCGTAGGTCAGACAAGGCGGACACGCGGCAGCCGACATCCAACGCCCGGGCCGGTTGCCTGATGGGACGCTGGCGCGt                                                                                                                                                                                                                <  1:973746/285‑1 (MQ=17)                                                    cgagaaagggagggggtaaGTACCAAGTATGGAAGGCAAAAGAGAGAATGGGCCACTGGGTGGTTCGGGAAACTGCCGGTAAGTGCGGCAGCTCGTGGTGCAGACAAGGCGGTCAGGCCGCATCCGACAGCCAATGCGCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCTGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                           <  1:383331/271‑1 (MQ=14)                                                    cAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                                           <  2:124442/288‑1 (MQ=255)                                                     aGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCCGTTGCCTGCTGCGACGCTGGCGCGTCGTCTCa                                                                                                                                                                                                         >  2:938606/1‑289 (MQ=255)                                                       aCAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTTCGCCACCCCGTAGGTCGGATAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCCGGcccggtgtctcag                                                                                                                                                                                                        >  1:202183/1‑279 (MQ=18)                                                       aCAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGCCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTCCCTGATGCGACGCTGGCGCGTCTTATCAg                                                                                                                                                                                                        >  2:798943/1‑288 (MQ=255)                                                              ggaggcgAAGTACCGCGCCAGGAAGGCGAAAGAGAAAATTGACACGGCGGGAGTTCGCGAAAGTGCCGTGAAGGCCGGAACCCGCTGGGTGAGAGAAGGCGGTCACGCAGCATCCGACAGCCAACCCCCGGGCCGGTTGCCGGATGCGACGCTGGCGCGTCTGATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                 <  1:430178/282‑1 (MQ=17)                                                              agagGTGAAGTACCGCGGCTGGAAGGGGGAAGAGCAAAGTGACCGCGCGGTAGTTGGCTAAGCTGCCGTGAAGTGCGGCACCCCGTAGGTCGGACAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                                                 <  1:1006427/288‑1 (MQ=17)                                                                aTGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAcc                                                                                                                                                                                               >  2:830938/1‑288 (MQ=255)                                                                aTGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGAAGCTGGCGCGTCTTATCAGGTCAACAcc                                                                                                                                                                                               >  2:638154/1‑288 (MQ=255)                                                                   tgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACAcc                                                                                                                                                                                               >  2:1006427/1‑285 (MQ=255)                                                                   tgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCAGCACCACGTAGGTCAGACAAGGCGATCACGCCGCCTCCGACATCGAACGCCCGAGACGGTTGCCTGATGCGACGCCGGCGCGTCTTGTCAGGCCTCCAcc                                                                                                                                                                                               >  2:712082/1‑285 (MQ=255)                                                                     aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGCCAGACAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATTAGACGCTGGTGCGTCTGCTGAGGCCGCAACCGCttga                                                                                                                                                                                         >  1:145741/1‑286 (MQ=18)                                                                     aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGCGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCtgt                                                                                                                                                                                          >  2:84860/1‑288 (MQ=255)                                                                             cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCt                                                                                                                                                                                  >  2:1068612/1‑288 (MQ=255)                                                                                    aaGGCGGGAGATAAAATTGATGACGCGGGAGTCAGGTAAAGGCCCGGGAAGTGCGGCACCGGGTAGGCCGTCCAAGGCGGGCAGTCCGCATGCGACTGACCACGCCCGAGCGGGGTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                                                                                                                                                           <  1:556608/288‑1 (MQ=14)                                                                                        cGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGtt                                                                                                                                                                                                                                      >  2:634878/1‑225 (MQ=255)                                                                                        cGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAgg                                                                                                                                                                      >  2:721163/1‑289 (MQ=255)                                                                                                 aaaTGGACAACGCGGGAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGAAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGACCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATaa                                                                                                                                                              <  2:673603/288‑1 (MQ=255)                                                                                                                          aaCTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                                                                                                                    >  2:709075/1‑289 (MQ=9)                                                                                                                            ctcgcgggAAGGGCGGCGCCCCGTAGGGTAGACAAGGCAGTGAGGCAGCGTCGGACATCCACCGCGCGGGCCGGTTGCCGAATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAAc                                                                                                                                  <  1:830938/282‑1 (MQ=14)                                                                                                                                gtgaaataggggcgtctgcgggcggggACAGGGTGGTCAGGCAGCGTCGGACAGCAGACGCCCGGGCAGATTGCGGGAGGCGACGCTGGCGCGTCGTAGGAGGCCTAAACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGcc                                                                                                                               <  1:302230/262‑1 (MQ=14)                                                                                                                                        gcggacccccGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGGGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCgg                                                                                                                       <  2:77296/282‑1 (MQ=255)                                                                                                                                                      ctcatccaacgcgtcccgcccccacccGACATCCAACACCCGGCCCGGTTTCCTGAGGCGACGCTGGCGTGTCTTATCAGGCCTACACCGCTGTGACGTGCGGCACCCCGTAGGTCAGACAGGGCGGTCACGCCGCATCCTACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTCTCAGGCCTACACCGCTGTCAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                                        <  2:995915/264‑1 (MQ=14)                                                                                                                                                                gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGATGGcgcg                                                                                                >  2:88219/1‑287 (MQ=255)                                                                                                                                                                     cgcgcccgatgcGAAATCGAGCGCCCGAGCGGGTGGCCGAGTGGGAAGGTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACAACCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                          <  1:721163/277‑1 (MQ=14)                                                                                                                                                                     cACGCCCCATCCGACATCCAACGCCCGAGCGGGTTGCCGGATGGGACGCTGGCGCGTCTTATCAGGCCTAAACAGCTGTGAAGGGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATc                                                                                         <  2:872496/289‑1 (MQ=255)                                                                                                                                                                            cctCCTACATCCAACCCCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGCCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAATTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTa                                                                                  <  2:125066/287‑1 (MQ=255)                                                                                                                                                                            cATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTTCCCCACCCCGTAGGCCGGATAAGGCGGTTACGCCGCATCCGACCTCCAACGCCCCAACCGGTTGCCCGGTGCGCCGCTGGCCCCTGTTAGCAGGCCt                                                                                   >  1:1088125/1‑288 (MQ=14)                                                                                                                                                                              tcgcagtaccAACGACCGCGCCGGTGGCCTGCCGCGCGGCTGGCGCGTCTTACCAGGCCTACACCGCTGGGAGGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                <  2:890181/281‑1 (MQ=14)                                                                                                                                                                              tCCGACTTCCAACGCCGGAGCCGGTGGCCTGATGCGACGCTGTCGCGTCTAAGCAGGCCGAGACGGCGGTGAAGGGAGGCACCCGGTAGGTCGGACGAGGCGGTCACGCCGCATCCGACATCCAACGCCCGGGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                <  1:98973/289‑1 (MQ=11)                                                                                                                                                                              tCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                <  2:106829/289‑1 (MQ=40)                                                                                                                                                                                               gggCCGGTGGCCTAATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCCGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACACCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTcc                                                               <  2:52460/287‑1 (MQ=11)                                                                                                                                                                                               gagaCGGTTGCCTGATGCGACTCTGGCGCATCTTATCAGGCCTAAACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTcc                                                               <  2:280040/289‑1 (MQ=255)                                                                                                                                                                                                                                       ccTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGACGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCCTCCGACa                       >  2:31995/1‑289 (MQ=255)                                                                                                                                                                                                                                        cTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                          >  2:1008126/1‑221 (MQ=255)                                                                                                                                                                                                                                           cacCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                   <  2:326154/289‑1 (MQ=255)                                                                                                                                                                                                                                            acCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCa                   >  2:50300/1‑288 (MQ=40)                                                                                                                                                                                                                                                cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGTTAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCCGGCCTACCCCGCTGCGCAATGCTCCACCCCCTAGCTTCGACAAAGCCGTTACGCCCAAGCCGACATcaatcgc               >  1:249514/1‑282 (MQ=11)                                                                                                                                                                                                                                                        tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGCTGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGTTAAGTGCCCCACCCCGTAGTTTCGATTAGGCGGCTACCCCGCCTCCCActaccacgcaccgcgccgg      >  1:164201/1‑271 (MQ=11)                                                                                                                                                                                                                                                        tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                               >  1:1101/1‑200 (MQ=255)                                                                                                                                                                                                                                                        tGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                               <  2:1101/200‑1 (MQ=255)                                                                                                                                                                                                                                                          cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTACCAGGCCTACACCGCTGTGAAGTGCCCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGATATCCAACGTCCGATCCCGtt    >  2:938422/1‑289 (MQ=255)                                                                                                                                                                                                                                                            gcatctgttAGGTGAGGTAACGCGGTCAGTCTGCACCGTACATTCGACGCGCGAGCCGGTGGCCCAATGCCACGGGGGCGCGGCTTAGCGGTTCTAAACCGCTGTGAAGTGCGCCACCCCTTAGGTCTGATAAGGGGGTGCCGCCGCGTCCAACATCTCCCGCCCAAGCTGGTTGCGTGGTGCAACGCTGGCGCGTCTGCTCAGGGTTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACGCCCGAGCCGGTTGc  <  1:211588/281‑1 (MQ=2)                                                                                                                                                                                                                                                                  cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                          <  1:596940/195‑1 (MQ=38)                                                                                                                                                                                                                                                                  cGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                          >  2:596940/1‑195 (MQ=38)                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                   ATGCTGGAATGCTGCGTTACAACAGCAAGCTGCCGATCGTGGTGTATACGCCAGACAATGTAGATGTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGC  >  NC_000913/2304270‑2304809

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