Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	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,380	1	.	C	96.0%	143.1 / ‑1.1	50	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (26/22);  total (27/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	4,296,380	2	.	G	96.0%	139.4 / ‑1.4	50	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (26/22);  total (27/23)
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.98e-01

TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCG  >  NC_000913/4296153‑4296664                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                                                                tacgCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGATAAGGCGTTTACGC‑‑CGC‑AGCCGCCATCAATTCCTGACGCGACGCTTGCCGCCTCCTGTCAGGCCTATCTTAACCTTTGGTTAAt                                                                                                                                                                                                                       >  1:141078/4‑301 (MQ=11) ctattgcgctgtctggcgtaggCCTGATAAGACGCGGCAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCAGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTCAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGCACTCAATGCCTGATGCGACGCTTGCCGCGTGTTATCAGGCCTATCTTAACCGTTGGTGAAt                                                                                                                                                                                                                       >  2:244181/22‑301 (MQ=11)     ccGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCCTTGGTTAAttttt                                                                                                                                                                                                                   >  2:485899/1‑301 (MQ=255)      ggaaccgctgagagaggCCTGACGGGTCGCTTGACCCCGCTTACCCGGCCTAGGCCGGAAGGCGAATGTGCCTGATTTAGCGGAATTTGCTAGGTCTGAGAAGTGGTTTATGCCGCATCCGTCATCATCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                                                                                                                                                                                                                  <  1:246554/285‑1 (MQ=11)         aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCCACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTCTCCGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                                                                >  2:471780/1‑300 (MQ=255)          tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGA‑TTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                                                                                                                                                                                                             >  2:106541/1‑301 (MQ=255)                aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTcc                                                                                                                                                                                                        >  2:166340/1‑300 (MQ=255)                       ccggATGCGACGCTGGCGGCGTGTTATCAGGCCTACGGCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCgg                                                                                                                                                                                                 <  2:233751/298‑1 (MQ=14)                           aTGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGGCATAACAGGCCTATCTTAACGGTTGGTTAATTTTTTAAACCATTCCCACCCgggttc                                                                                                                                                                                             >  2:328109/1‑297 (MQ=11)                               gcgtcttgccgcgTCTTATCAGTCGTACGCCAGACAGCGGAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                                                                         <  2:195119/292‑1 (MQ=14)                               gACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGCCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc                                                                                                                                                                                         <  2:47239/301‑1 (MQ=255)                                      gccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                                                                                <  2:30249/269‑1 (MQ=25)                                      gccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAAACATCGCAATAGCTTGATTTAGCGTGATTTTGTAGGTTGGACAAGGCGGTTA‑‑CG‑CGCTATCTGACCTCAATGCTTGTTGCAAAGCTTGCCGCTTCTTTTCAGGCCTATCTTGACCTATTGTTAATCTTTCtaaattcggaatagaggcgtgttgggtaag                                                                                                                                                                                    >  1:30243/3‑272 (MQ=11)                                          cgTCTTAGCACGCCTAGCGAAGACAGCGAAATAGCCTGATTCAGCGTGATTTGGTGGGGCGGATAAGGCGTTAATGCGGAAGCAGACATCAACGCCTGATGCGACGCTTACCGCTTCTTATCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTATAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACa                                                                                                                                                                              <  2:124367/301‑1 (MQ=17)                                                aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGCCGCTTGCCTCTTCTTATCCGGCCTATCTTAACCGTTGGGTAATTTTTCAAACTCTTCAACCCCGGCATCAAACCCTCCCATAccaccca                                                                                                                                                                         >  1:210077/1‑294 (MQ=14)                                                 cgaGGCCTCCGCCTGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                                                                                                                                                                       <  1:243164/299‑1 (MQ=255)                                                  cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACCGCTCAGGt                                                                                                                                                                      >  2:293709/1‑301 (MQ=255)                                                      cctgcgcaaaccAGCGACATGGCCTAATTTAGCGGAACTTTGTAGGCCGACTAGGGCTTTTAGGCCGAATCCGACATCAAGGCCTGATGCGACGCTTAACTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGA‑TTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAAc                                                                                                                                                                 <  1:106510/290‑1 (MQ=17)                                                              aGACATCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCa                                                                                                                                                          <  2:182238/301‑1 (MQ=255)                                                               cccAGCGCAATGCCCTGGTTTATAGTAGTTTTGTAGTTCGGATAAGCGGTTTATGCCACATCCAACATCAACGCCTGATTCGACGCTTACCCCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAg                                                                                                                                                         <  1:527602/299‑1 (MQ=18)                                                                 cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACACGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCATTGCAg                                                                                                                                                         >  2:527751/1‑299 (MQ=255)                                                                                aTTTAGTGTGATTTTGTAGGTCGGATAATGCTTTAATCCCGCTTCCGACAGCACCGCCTAATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATGGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATAAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc                                                                                                                                         <  1:485767/300‑1 (MQ=18)                                                                                    ggCGTGATCTTGTGCGTCGGATAACGCGTTTACGACGCATCCGACATCAACGCCTGATGCGACGCGTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                                                                    <  1:3582/300‑1 (MQ=18)                                                                                     gcggatattttcgtgcttgatcgtccGTTTATGCCGCCACCGCACTCCACGCCGTGTGCGAGGCTTACCGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGAGAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTtctc                                                                                                                                    <  1:328007/275‑1 (MQ=21)                                                                                     gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCCGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATAc                                                                                                                                                                               >  2:535489/1‑257 (MQ=255)                                                                                         gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTcc                                                                                                                                                                                   >  2:192994/1‑249 (MQ=255)                                                                                         gATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTcc                                                                                                                                                                                   <  1:192941/249‑1 (MQ=255)                                                                                                           aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTa                                                                                                             <  2:141116/301‑1 (MQ=255)                                                                                                            ggCGTTGATGTCGGATGCGGCATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTAAAACCTTCTTCACCGTTTa                                                                                                             >  2:665651/1‑300 (MQ=21)                                                                                                             gCGTTGATGTCGGATGCGGCATAAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCACCCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGTTaaaa                                                                                                                             >  2:611855/1‑283 (MQ=21)                                                                                                                  tctGCCCCATCCGGCGTAAACGCCTGAAGCAAGGCTTAAACGGTCTTATTAGGGCTACGGCAGACAGCGCAATAGCCTTATTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCGGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTAGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAca                                                                                                       <  1:619841/298‑1 (MQ=25)                                                                                                                     gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGTTAAGGCGTTTACGC‑‑CGC‑ATCCGACATCCCTGCCTGATGCGACGCATGCCGCGTATAATCAGGCCTATCTTAACCGTTGGTTAATTTTTCTACCTCTTCACACCCGGGCTGAAACCCTTCAATACCACTCAGCTTCAACCTGTGACGCGCCTTTGGATTGATCTTCTCGACAATCActtcttagactctatagatcattat                                                                                                   >  2:310121/1‑282 (MQ=17)                                                                                                                               aCATCAACGACTTGTGCGACGCTTGACGCGTCTTATCGAGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATTCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAg                                                                                         <  1:306864/301‑1 (MQ=255)                                                                                                                                 aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttc                                                                                                                                         >  2:528034/1‑251 (MQ=255)                                                                                                                                              gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCTGCCCAGTACTCGGAATAACCg                                                                          >  2:654464/1‑301 (MQ=255)                                                                                                                                                gtcgtattcccgttCTTGTCAGGCCTACGACAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTA‑‑CGACGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGt                                                                        <  1:463164/288‑1 (MQ=37)                                                                                                                                                gccGCTTAACGCTTCTTATCAGGCCTACGACAAACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGt                                                                        <  1:665466/299‑1 (MQ=255)                                                                                                                                                       aaCGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGa                                                                  <  2:210134/300‑1 (MQ=255)                                                                                                                                                                 aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCATTACTCGGAATACCCGGTTCGGGAAATTGCagaac                                                       >  1:131527/1‑299 (MQ=255)                                                                                                                                                                           cgcccagcagttcAATAGCCTGATTTAGCGTGATTGTGTAGGTCTGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAGTGCCTGATGCGAAGCTTGCCGAGTCTTATCAGGCCTATCTTAACCGTTGGTTAAATTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGaa                                             <  1:293615/289‑1 (MQ=255)                                                                                                                                                                           cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAAATCATTCCCGCCCAGTAATCGGAATAACCGGCTCGGGAAATTGCAAATCTGCGTTTGa                                              >  1:518430/1‑300 (MQ=255)                                                                                                                                                                           cGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCCCGATAAAACCTTCTTCACCGTTTAAGAACATAATTCCCGCCCAGTAATCGGAATAACCGGTTCGGGAAATTGCAGAGCTCCGTTTGa                                              >  2:352414/1‑300 (MQ=255)                                                                                                                                                                                             ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTGTAAGAACATCATTCCCGCTCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCTTTTGAAATACCTGGTGTCCttttt                           >  1:53445/1‑298 (MQ=255)                                                                                                                                                                                                        gTGATTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGACTAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCt                >  2:230090/1‑301 (MQ=255)                                                                                                                                                                                                           aTTTTGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGGTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTat             >  1:324819/1‑300 (MQ=255)                                                                                                                                                                                                             tttGTAGTTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAAt            <  2:518570/300‑1 (MQ=255)                                                                                                                                                                                                             tttGTAGGTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTGTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCTTTCCCGCCCATTACTCGGCATAACCGGTTCGGGAAATTGCAAAGCTGCGTTGGAATTACCAGGTTGCCTTGTCTTCTTCTGCCTTAAtt           >  1:559228/1‑300 (MQ=38)                                                                                                                                                                                                                   ggTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAAt                                                               <  2:311883/243‑1 (MQ=255)                                                                                                                                                                                                                    gTCGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGc    >  1:111836/1‑301 (MQ=255)                                                                                                                                                                                                                      cGGATAAGGCGTTTACGC‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCg  >  2:133989/1‑301 (MQ=255)                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                                                                TATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTA‑‑C‑‑CGC‑ATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGTTTGAAATACCAGGTTGCCTTGTCATCGTCGGCCTTAATGCCCACGCCG  >  NC_000913/4296153‑4296664

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