Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,190	A→C	43.9%	intergenic (+397/+247)	gltP → / ← yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,190	0	A	C	43.9%	‑3.2 / 40.0	36	intergenic (+397/+247)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (9/11);  new base C (12/4);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 9.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT  >  NC_000913/4295910‑4296464                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                    ggTCGGATAAGGTGTTTATGCCGCATACGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTAGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                                                                                                                                                                                                                                                                                     <  2:1149394/281‑1 (MQ=37)     ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGt                                                                                                                                                                                                                                                                                <  2:1789639/282‑1 (MQ=37)           tgcggataggccccatcaggaaatcaCGCTTGATGCGAAGCTACACGCGGCTTATCAGGTCTACGCCAGACGGCGTAATAGCATGATTTCGGGTGATTTTGTAGGGCGGATAAGGGGTTTCCGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                                                                                                                                                                           <  1:836760/256‑1 (MQ=11)           ggCGGTTGTGCCGAAGCAGACGTCAAGGCTGGAGGCGACGGTAACCGCGGCTGATCAGGGTTATGCCAGACAGCGCAATAGCATAATTTGGAGTGATTATGTAGGTCGAATAGTGCGTTTACGCCGCGTCCGGAATTAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGTAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt                                                                                                                                                                                                                                                                           <  1:338921/281‑1 (MQ=11)               tttATGCCGCGTCCGCCATCACCGCCTGATGCGACGCTTAACGGGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGTCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                      <  2:1573518/282‑1 (MQ=255)               tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                      <  2:787068/282‑1 (MQ=37)               ttattGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                      <  2:639960/278‑1 (MQ=25)                ttATGCCCCGACCGACATCAACGCCTGAGGCTACGCTTAGCGCGTCTTGACAGGCCTACGCCAGACAGCGCAATGGCCGAATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                      <  1:1140106/281‑1 (MQ=11)                 tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACTGCGTTTACGCCGCATCCGACATCAACCCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGATTTTTTAGGTCGAATAATCCGTTTATGCCCCCTCCGACATCAACGCCTGATGCCAAGCTTTAcgtctccatccaggcat                                                                                                                                                                                                                                                                    >  1:44209/1‑267 (MQ=11)                      cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGCTGCGACGCTTAACGCGTCTTATCAGGCCTACGc                                                                                                                                                                                                                                                                >  2:1331332/1‑281 (MQ=25)                                     cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGACAAGGCGTTTATGCAGCATCCGACACAAACCCCTGATGCGACGCTGCACGCTTATTTTCAGGCCTACACAAGaggccgccaacgc                                                                                                                                                                                                                                                >  1:576397/1‑270 (MQ=11)                                       ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGcc                                                                                                                                                                                                                                              >  2:1918506/1‑281 (MQ=35)                                                aCGCTTACCGCGTATTATCAGGGCTAAGACAGACAGAGCAGTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGACGCATCCGAAATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCg                                                                                                                                                                                                                                    <  2:1916351/282‑1 (MQ=255)                                                            tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACATCGCAATAGCCTGAATTAGCGTGATTTTGTAGGTCGGATCATGAGTTTATGCCGCATCCGACTACAACGCTTTATGCGACGCTTAACGCTTCTTCTACGGCATACCTCAGACAGCGC‑CACACCCTTATTTAgacggtatattgagg                                                                                                                                                                                                                        >  1:1198959/1‑268 (MQ=11)                                                                          cGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGATTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACAGCGTCTTATCAGGCCTACGCCAGACAGCGC‑CATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAGAAGGCGt                                                                                                                                                                                                           >  2:1730324/1‑281 (MQ=255)                                                                          cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                                                                                                                                                                                           >  2:1411700/1‑281 (MQ=33)                                                                                    cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTCTAGGTCGCATAAGTCGTTTACGCCCCCTTCGACACATATGCCTGATGCGACGCTTTCCGCGTCTTATCAGTCCGACGCCAGCCAGCCCTTAAACCATAATTTA‑CGTCATTTTTTACGACGGACAAGTCGTTTaacgcgca                                                                                                                                                                                                 >  1:107710/1‑274 (MQ=11)                                                                                           gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                                                                                                                                      <  2:1377318/206‑1 (MQ=34)                                                                                              tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTa                                                                                                                                                                                                                                                                   >  1:1377318/1‑206 (MQ=21)                                                                                                                               ttttACGCGGCATCGCAAATCAACGCCTGAGGAGAGGCTTGCCGGGTCTAATCAGGCCTACCCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGAGAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCa                                                                                                                                                                                                                                                              <  1:224477/177‑1 (MQ=11)                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAc                                                                                                                                                                                                                                                           >  2:224477/1‑178 (MQ=35)                                                                                                                                  tACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAACGCGGTTATGCCGCATCCGACATCCAGGCCTGATGCGCCGCTTCACGCCGCTTGGCAGGTCTACACCAGACAGCTG‑AACAGCCTGTATTAGCGTCATTTTGGAGTCCGGATAAGGCTATTCCGCCGCTTCCGACCTCAACGCCTTACGCAACTCttgccgccgcctctccggcct                                                                                                                                                  >  1:628637/1‑263 (MQ=2)                                                                                                                                  cccaccgcggaagtggggatggCCTGAGGCGACGCTTGCCGCTGTTTATCAGGCCTACGCCGGACAGCGCATTGGCCTGATTTGGCGTGATTTTGTGGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc                                                                                                                                                    <  1:1907248/259‑1 (MQ=11)                                                                                                                                                                                                                                                 tttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTATACGCCGCATCCGACATCaatgcctgatgcgacgcttgccgcgtcttatcaggc                                   >  2:1803756‑M2/1‑246 (MQ=255)                                                                                                                                                                                                                                                 tttACGCCGCATCCGACATCAACGCCGTATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGA‑AATAGCCTGATTTGGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatgcctgatgcgacgcttgccgcgtcttatcaggc                                   >  2:587182‑M2/1‑246 (MQ=255)                                                                                                                                                                                                                                                  ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaatgcctgatgcgacgcttgccgcgtcttatcaggcc                                  >  2:197698‑M2/1‑245 (MQ=255)                                                                                                                                                                                                                                                   tacgCCGCATCCGACATCAACGCCTGATGCGCCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTTGCGTGCTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACCTCAACCCCTGATGCGACGCTTAACGCGTCTTACTCGGCCTACGCCAAACCGCCCAACCACCTGATTTCGCGTGATTTTGTATGTCGGATACGGACTTTACGCCGCCTCCGACATCCATCCCTGTTTCGACGCTTGGCGCGTTTTAGCAGGcc                                  >  1:1964760/4‑281 (MQ=1)                                                                                                                                                                                                                                                        cGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTACGGCGTTTACGCCGCATCCGACCTCAATGCCTGATGCGACCCTTGCCGGGTCTTATCAGGCCTATCt                             >  1:357935/1‑281 (MQ=255)                                                                                                                                                                                                                                                                    cagcGCCTGATGCGGCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACGGCGC‑AATTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                <  2:1198959/279‑1 (MQ=255)                                                                                                                                                                                                                                                                             aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTGATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       <  2:1862442/282‑1 (MQ=255)                                                                                                                                                                                                                                                                              ttcGACGCTGGCCGCGTCTTATCAGGCCTAAGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAATTTTCaaa      <  2:1063487/280‑1 (MQ=255)                                                                                                                                                                                                                                                                              tGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa       >  2:771663/1‑281 (MQ=38)                                                                                                                                                                                                                                                                               gCGACGCTTGCCGCGTCTTATCATGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa        <  2:1310987/279‑1 (MQ=255)                                                                                                                                                                                                                                                                                  aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGCCACCAATGCCTCATGCGAACCTTGCCGCGCCTTTTCAGACCTATCTTAACCGTGGGTTAATTTTTCAAActtt  >  1:1310987/1‑280 (MQ=11)                                                                                                                                                                                                                                                                                  aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   >  2:1429425/1‑281 (MQ=255)                                                                                                                                                                                                                                                                                  aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc   >  2:1299477/1‑281 (MQ=255)                                                                                                                                                                                                                                                                                  aCGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAAGCTTAACGCGTCTAATCAGGCCTACGCCAGCCAGCGCCATAGCCTGATTTAGAGTGAGTTTGTAGGTAGGATAAGGCGTCTCCGCCGCATCCGACATCACAGCCTGATGCGAGTCTTGCCGCCTCCGATCAGGCCCATGTTAGccatcgccatatcctcatgtcgtg  >  1:1799660/1‑260 (MQ=11)                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                    GGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGC‑AATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCT  >  NC_000913/4295910‑4296464

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