Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,296,191	A→C	21.6%	intergenic (+397/+245)	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,191	0	A	C	21.6%	56.6 / 19.6	41	intergenic (+397/+245)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base A (16/14);  new base C (4/5);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

TGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC  >  NC_000913/4295907‑4296471                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            tGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATAAACGCCTGATGGGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGAGAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc‑g                                                                                                                                                                                                                                                                                        <  2:405373/289‑1 (MQ=19)   tgagTCGGATAAGTCGGTGATGCCCCATACGAAATCAACGCATGATGCGACGCTTACCGCGTCTTATCAGGCCTACGGCAGACGGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑Gt                                                                                                                                                                                                                                                                                       <  2:254700/285‑1 (MQ=11)        cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTAt                                                                                                                                                                                                                                                                                  >  2:104805/1‑288 (MQ=37)        cGGATAAGGCGGTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTATGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTA‑TGCCGCATCCGACATCCACGCTTGATGGGACGATTACCGC‑GTCTTAt                                                                                                                                                                                                                                                                                  >  1:706332/1‑288 (MQ=11)        cGGAAAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTAt                                                                                                                                                                                                                                                                                  <  1:123811/288‑1 (MQ=37)            tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACTCTTATCGGGTCTTATCGGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  2:95668/289‑1 (MQ=255)              cggggTTTATGCCGAGACCGACATCAACTCTTGATGTGACGATTACGGCGACTGATCAGGGCTACGCGAGACAGGGCAATGGTCTGATTTAGCGGGATTTTCTAGGTCGGATAAGGGGTTTACGCCGCATCCGAGATCAACGGCTGATGCGACGCTTGCGGCGTCTTGTCATGCCTACGCCAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  2:352705/284‑1 (MQ=11)                gtgTTTTTGCCCAACCCGATATCAATGGCTGATCCAACGCATAACGCGTCTTATCAGGACGACGCCAGACAGCGCAATAGCCTAATTTAGCGTTATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  2:661862/283‑1 (MQ=11)                      aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGa                                                                                                                                                                                                                                                                   <  1:104805/289‑1 (MQ=37)                          cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCCACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             >  2:285754/1‑275 (MQ=25)                          cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGc                                                                                                                                                                                                                                                                             <  1:285754/275‑1 (MQ=35)                            cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATATCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTA‑GGCCGCATCTGACATCAACGCCTGACGCCACGCTTTACGC‑GTCTTATCCGGCCTACGCCAGACAgccc                                                                                                                                                                                                                                                             >  2:181018/1‑287 (MQ=255)                            cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACCTCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAgcgc                                                                                                                                                                                                                                                             >  1:99009/1‑289 (MQ=25)                              tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                                                                                                                                           <  1:370589/289‑1 (MQ=35)                                      cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTATGTCGGATCAGGCGTTTA‑TGCCGCATCCGACCTCAACGCCTGATTCGACGCTTAACGC‑GTCTTACCAGGCCTACCCCAGACAGCGCAATATCCTg                                                                                                                                                                                                                                                    >  2:720649/1‑288 (MQ=255)                                                    aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCTCCGCTTAACGT‑GTCTTATCCGGCCTACGCTCGACAGCGCAATAGCCTGATTTAGCGTGAgttt                                                                                                                                                                                                                                     >  2:629602/1‑289 (MQ=255)                                                    aCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtttt                                                                                                                                                                                                                                     >  1:302861/1‑289 (MQ=35)                                                                            tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCCTGATTTTGTAGGTCCGATAAAGCTTTTacg                                                                                                                                                                                                              >  2:553792/1‑286 (MQ=12)                                                                              cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTAGGTCGGATAAGTCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCCAACAGCGCAATAGCCTGTATTAGCTTGCTTTTGTATGTCGGATCAGGCGTTTACGcc                                                                                                                                                                                                            >  2:440833/1‑288 (MQ=11)                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACATCGCACTCGCCTGATTTAGCTTTATTTCGTAGGTCGGATAAGGCGTTTACGCCGCACCCGAcc                                                                                                                                                                                                  >  2:528212/1‑287 (MQ=255)                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAc                                                                                                                                                                                                                                                                                            >  1:53445/1‑199 (MQ=34)                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACTCCAGACCTCGCAATAGCCTGATTTAGCCTGATTTTGTAGTTCGGATAAGGCGTTTCCGCCCCATCCGaac                                                                                                                                                                                                  >  2:350646/1‑286 (MQ=255)                                                                                        cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCTCGACAGCGCAATAGCCTGATTGCGCGTTATTTTGTATGTCCGATCATGCCTTTACGCCGCATTCGACAt                                                                                                                                                                                                 >  2:540342/1‑289 (MQ=11)                                                                                            aTAGCCTGATTTAGCGGGAGTTTGTGGGTCGGTTAAGGCGTTTTCGCCGGATCCTACGTCAACGCCTGATGCGACGCTTACCGGGTCTTATCCGGCATACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑CGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCaat                                                                                                                                                                                             <  2:587549/289‑2 (MQ=11)                                                                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                                                                    >  1:223955/1‑287 (MQ=17)                                                                                                                           gATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACa                                                                                                                                                                                                                                                                 <  1:507251/190‑1 (MQ=35)                                                                                                                                     ttACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGc                                                                                                                                                                                                                                                       >  2:507251/1‑190 (MQ=35)                                                                                                                                                             cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                                                                  <  1:390839/235‑1 (MQ=14)                                                                                                                                                             cTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGAATCCGACATCAACGCCTGATGCGa                                                                                                                                                                                  >  2:390839/1‑235 (MQ=14)                                                                                                                                                                                               cgcgagagcgcgcAATAGACTGATTTAAGGGAATTGTGTGGGTCGGATACAGCGTTTA‑TGCCCCGACCGACATGAACTCCTGAGGCGACGCTTATCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑cgccg                                                                                        <  2:296055/280‑5 (MQ=11)                                                                                                                                                                                                         cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                          >  2:44176/1‑231 (MQ=33)                                                                                                                                                                                                         cgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCaa                                                                                                                                          <  1:44176/231‑1 (MQ=33)                                                                                                                                                                                                                     aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑cg                                                                                           <  2:623955/264‑2 (MQ=14)                                                                                                                                                                                                                                              tAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                           <  1:609434/289‑1 (MQ=14)                                                                                                                                                                                                                                                        tacgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                  <  1:553792/286‑1 (MQ=2)                                                                                                                                                                                                                                                              gCATCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                            >  1:581922/1‑289 (MQ=12)                                                                                                                                                                                                                                                                aTCCGACATCAACGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACTCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                          >  2:700199/1‑289 (MQ=255)                                                                                                                                                                                                                                                                            ggCCTGGTGCGACGCTT‑GCGCGGTCTGATTAGGCCGAAGCCAAGTAGCTAAATAGCCTGATTTAGGGTGATTTCGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCGGATGCGACGCTTAACGCGTCTTTTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa              <  2:430554/288‑1 (MQ=11)                                                                                                                                                                                                                                                                            cGCCTGATGCGACGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa              <  2:302861/289‑1 (MQ=255)                                                                                                                                                                                                                                                                            cGCCTGATGCGACGCTGGCCGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa              <  1:528212/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                        cGCTTGCCGC‑GTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACAccc  <  1:181018/289‑1 (MQ=255)                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                            TGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑TGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGC‑GTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCC  >  NC_000913/4295907‑4296471

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