Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,307,909	G→A	7.9%	intergenic (+182/+533)	eco → / ← mqo	serine protease inhibitor ecotin/malate dehydrogenase (quinone)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,307,909	0	G	A	7.9%	137.7 / 3.3	64	intergenic (+182/+533)	eco/mqo	serine protease inhibitor ecotin/malate dehydrogenase (quinone)
Reads supporting (aligned to +/- strand):  ref base G (27/30);  new base A (2/2);  total (32/32)
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.49e-01

GTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCC  >  NC_007779/2307671‑2308120                                                                                                                                                                                                                                               |                                                                                                                                                                                                                    ggaaaGAACCGCGTCTGGAAGGCGAAAGAGAAAAGGGCCAGCGCGGTGGTTCGCAAACCTGCCGTGAAGGGCGGCACCCGGTGGGTCAGCCAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCGGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCg                                                                                                                                                                                                                     <  1:1428903/236‑1 (MQ=17)  tgAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCTGCCGACATGCAACCCCCGa                                                                                                                                                                                                                    >  2:393074/1‑239 (MQ=255)    aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCCCGCCGCATCCGACATCCAACGCCCGAg                                                                                                                                                                                                                   >  2:1679859/1‑238 (MQ=255)    aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAg                                                                                                                                                                                                                   >  1:1627228/1‑238 (MQ=255)    aaGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAA‑‑GCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCTAACCCCCAAGccc                                                                                                                                                                                                                >  1:944958/1‑238 (MQ=17)      gTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCg                                                                                                                                                                                                                >  1:587482/1‑239 (MQ=255)      gTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTCCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCg                                                                                                                                                                                                                >  1:1385866/1‑239 (MQ=255)         ccGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGtt                                                                                                                                                                                                             >  1:190838/1‑239 (MQ=255)            cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCTCATCCGACATCCAACGCCCGAGCCGGTTGc                                                                                                                                                                                                           >  1:503244/1‑238 (MQ=255)            cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGCGCCGGTTGc                                                                                                                                                                                                           >  2:868924/1‑238 (MQ=255)            cgTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCAAGCCGCATCCGACATCCAACCCCCGACCCGGTTGc                                                                                                                                                                                                           >  2:704210/1‑238 (MQ=255)              tCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTg                                                                                                                                                                                                        >  1:1405572/1‑239 (MQ=255)                              aaaaaCTGCCAAGGCGGAAGTGCGCAAAACTGCCG‑GGAGCGCGGCCCCCGGTAGCTCAGACCGGGGCGTCACGCCGAGCCCGACAGCAACCGCCCGGGCCGGTTGCCGGATGCGACGCTGGCGCGCCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                                                                                                                       <  2:1973591/238‑1 (MQ=14)                                 aaTTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGCGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCCAGCCGGCTGCCTGATGCGACGCTGGCGCGTCt                                                                                                                                                                                     >  1:615293/1‑239 (MQ=255)                                    tGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCCGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGCCATCCAACCCCCGACCCGGTTGCCTGATGCGACGCTGGCGCGACTTa                                                                                                                                                                                   >  1:205210/1‑238 (MQ=255)                                         aCGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCGGCTAAGGCGGTCACGCCGCCTCCGACATCCAACGCCCGAGCCGGTTGCCCGATTCCCCCCTGgcgccccctaccaag                                                                                                                                                                              >  2:1665513/1‑227 (MQ=14)                                         aCGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAATGCGGTCACGCCGCATCCGACATCCAACGCCCCAGCCGGTTGCCTGATGCGACGCCGGCGCGTCTTATCAga                                                                                                                                                                              >  1:1788992/1‑237 (MQ=255)                                             ggTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGACTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGCGAAGAGCGGCACCCCGCAGGTCAGACAAGGCGGTAACGCCGCCTCCGACATCCAACGCCGGAACCGGTTGCCGGATGCGACGCCGCCGCGCCTCATCAGGCCTa                                                                                                                                                                          >  1:980273/1‑238 (MQ=14)                                               tggtgcGCTAACCCGCGGTAGAGTGCGCCCCCCCGTAGGGCGGACAAGGCGGCCCCCCCGCATCCGCCCTCCAACGCCGGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                        <  2:980273/233‑1 (MQ=14)                                                  gtCGCAAAACTGCCGGGAAGGGCGGCCCCCCGAGGTTGAGTGGGGGGGGCAACGCCGCGTCCGACTTCCAACGCCCGAGCCGGTTGCGTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCg                                                                                                                                                                     <  2:1627228/237‑1 (MQ=14)                                                              ccGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCAACCGACATCCAACGCCCCAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTCATCAGGGCTACACCGCCGTGAACTGcg                                                                                                                                                         >  1:383877/1‑237 (MQ=255)                                                                      aggggCACCCGGTAGGTGAGACAAGGCGGCCAGGCCGCGTCCAACAGCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTAATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCg                                                                                                                                                 <  2:663668/235‑1 (MQ=14)                                                                                          gagaGGCGGCCGGGCCGATTCGGACCTCCAACGCCCGGGCCGGTTGCCTGTGGCGACGCTGGGGCGTCTTAGCAGGCCTACACCGCTGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGtt                                                                                                                              <  2:2066788/234‑1 (MQ=14)                                                                                           cAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCGGACAAGGCGGTCACGCCGCATCCGACATCCCACGCCCGTGCCGGTTGCCTGATGCGCCGCCCGCGCGTCTTCTCCGGCCTACACCGCTGTTAAGTGCTCCACCCCCTAGGGCGGACAAGGccggtcc                                                                                                                            >  2:850279/1‑232 (MQ=11)                                                                                             aGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCTCCCCGTAGGTCGGATAAGGCGGTTAcgcc                                                                                                                         <  1:1924935/239‑1 (MQ=255)                                                                                                cGGCCACGCCGCATCCGCCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgccgc                                                                                                                       <  1:850279/238‑1 (MQ=255)                                                                                                    cgcgccgcATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTTAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCTCATccc                                                                                                                  >  1:2171096/3‑238 (MQ=255)                                                                                                    cACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                                                                  <  1:704210/239‑1 (MQ=255)                                                                                                            cgccGACGTCCAGCGCCCGAGCGGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAAGTCGGATAAGGCGGTTACGCCGCATCCGACATCTa                                                                                                           <  2:588457/236‑1 (MQ=255)                                                                                                                     ccAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                        <  1:2116340/168‑1 (MQ=255)                                                                                                                               aGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGaa                                                                                                                                                              >  2:2116340/1‑168 (MQ=255)                                                                                                                                   ggTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                                                                                                                                  <  1:1572588/144‑1 (MQ=255)                                                                                                                                    gTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGc                                                                                   <  1:2053159/238‑1 (MQ=40)                                                                                                                                     ttGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCCCATCCCACACCCAACACCCCACCCGGGTTCCAGATGCg                                                                                  >  2:2053159/1‑238 (MQ=14)                                                                                                                                             tGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                                                                                                                                        >  2:1572588/1‑144 (MQ=255)                                                                                                                                                     tGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCCCCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTAt                                                                 <  1:734329/239‑1 (MQ=255)                                                                                                                                                       gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                            <  2:1695689/114‑1 (MQ=255)                                                                                                                                                       gcgcGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTgg                                                                                                                                                                                            >  1:1695689/1‑114 (MQ=255)                                                                                                                                                            tCTTATCAGGCCGACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCt                                                                                                            <  1:2081910/189‑1 (MQ=255)                                                                                                                                                               tATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGGAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                       <  1:1679859/239‑1 (MQ=255)                                                                                                                                                               agtCAGGCCTCCCCCGCTGTGAAGTGCGCCCCCCCGTAGCTCGGACAAGGCGGTCACGCCGCATCCGACATCCACCGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                       <  2:115695/237‑1 (MQ=255)                                                                                                                                                                ttgggcgcTACACCGGCGTGACGGCCGGCCCCCGGGGGGTCAGGCAGGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                       <  2:1819957/231‑1 (MQ=11)                                                                                                                                                                gtcgggCGTACACCGCGGTGAGGTGCGGAGCCCCGTAGGTCGGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTaca                                                       <  2:1187871/234‑1 (MQ=255)                                                                                                                                                                 tCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGCTAAGGCGGTTACGCCGCATCCGACACCTAACCCCCAAGCCCGTTGCCCGATGCGACGCTGGGCCGGCTTACCAGGCCTacac                                                      >  1:115695/1‑238 (MQ=14)                                                                                                                                                                     gCCTCCCCCGCGGTGAGGTGCGGCCCCCCGTAGGTCGGCAAAGGCGGGCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCg                                                    <  2:408016/236‑1 (MQ=11)                                                                                                                                                                            ccGCTGTGAACTGCGCCACCCCGTAGGTCAGAAAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTg                                          <  1:1541507/239‑1 (MQ=25)                                                                                                                                                                               cggggaggTGGGGCGCCCCATGGGTCAGACAGGGCGGTCACCCCCCATCCGACATCCAACGCCGGAGCCGGTTGCCGGAGGCGACGCTGGCGCGTCTTATCAGGCCTACTCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                        <  2:615293/231‑1 (MQ=11)                                                                                                                                                                               cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACCTCTCACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAgctgt                                        >  1:1991891/1‑234 (MQ=255)                                                                                                                                                                               cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACCACGCCGTGAAGGGcc                                        >  1:305675/1‑237 (MQ=255)                                                                                                                                                                               cTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGCTGGCCTGATGCGACGCTCGCGCGGCTTATCAGGCCTACCCCGCTGTGAAGGGccc                                       >  1:408016/1‑237 (MQ=255)                                                                                                                                                                                    aaggccGGCCCCCGGTAGGTCAGCCAGGGCGGTCAGGCCGCAACCGACATCCAACGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAA‑‑GCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                 <  2:944958/233‑1 (MQ=11)                                                                                                                                                                                      gTGCGGCCCCCCTTAGGTCAGACCAGGCGGTCACGCCTCAGCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCGGGCGCGTCTTATCAGGCCTACCCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCAcccc                                 <  2:1216423/238‑1 (MQ=255)                                                                                                                                                                                       tGCGGCACCCCGTAGGGCAGACAAGGCGTGAACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCg                                <  2:1788992/238‑1 (MQ=255)                                                                                                                                                                                         cGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAACACGGTTGCCTCATGCGACGCTGGCGCGTCTTATCAGGCCCACACCGCTGTGAAGAGCTCCAACCCGTAg                             >  1:1216423/1‑239 (MQ=255)                                                                                                                                                                                            cACCCCGGAGGCCAGGCGTGGCGGCCACTCAGCGTCCAACATCCGACGCCCGAGCGGGTGGCCTGAGGCGAAGCGGGCGCGTCTTATGAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCTTACACCGCTGTGAAGTGCTCCACCCCGTAGGt                           <  2:2060441/238‑1 (MQ=11)                                                                                                                                                                                                  gtggggcAGAAAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGAGGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATaa                    <  1:393074/233‑1 (MQ=16)                                                                                                                                                                                                     ggccggaaaaGGCGGTCACGCCGCATCCGCCATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAgg                  <  2:724085/230‑1 (MQ=12)                                                                                                                                                                                                         aGACAGGGCGGTCACGCCGCCTCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGt              <  2:1039646/238‑1 (MQ=18)                                                                                                                                                                                                         aGACAAGGCGGTCACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCCTCCGACATCTCACCCCCCAACCGGTTGCCTGATGCGACGCTGGCGCGGCTTCTCAGGCCTACCCCGCTGGGCAGTCCCCCCCCCCGTCGGTTGGATACGGCGGtt             >  2:619375/1‑239 (MQ=11)                                                                                                                                                                                                         aGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGGTGCCTGATGCGACGCTTGCGCGTCTTATCTGGCCTACACCGCTGTGAAGGGCTCCCCCCCGTAGGTTGGGTGAGGCggg              >  1:2116221/1‑237 (MQ=255)                                                                                                                                                                                                                gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCt                                                                                                                                                         >  1:239575/1‑92 (MQ=31)                                                                                                                                                                                                                gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCATGATGCGACg                                                                               <  1:372638/166‑1 (MQ=18)                                                                                                                                                                                                                gCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCATGATGCGACg                                                                               >  2:372638/1‑166 (MQ=18)                                                                                                                                                                                                                    ccccgccccATCCGACTTCCAACGCCCGAGCCGGTGGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATcc  <  2:1872099/236‑1 (MQ=9)                                                                                                                                                                                                                                  aCATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATg                                                                                    <  2:779706/143‑1 (MQ=17)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                    GTGAAGTACCGCGTCTGGAAGGCGGAAGAGAAAATTGACAACGCGGTAGTTCGCTAAACTGCCGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGGCACCCCGTAGGTCAGACAAGGCGGTCACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCC  >  NC_007779/2307671‑2308120

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