Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,296,381	+GC	intergenic (+587/+55)	gltP → / ← yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,296,380	1	.	C	90.5%	121.4 / 2.8	42	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base C (16/22);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 2.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	4,296,380	2	.	G	90.2%	116.9 / 2.6	42	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base G (16/21);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑A‑‑C‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGT  >  NC_000913/4296193‑4296616                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                gcgTCTTATCAGGCCTAAGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                                                                 <  2:152275/251‑1 (MQ=2)             gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                                                                                                                                                                     >  1:242198/1‑251 (MQ=25)                        aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGTCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                                                                                                                          <  1:118189/251‑1 (MQ=255)                         cAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                                         <  1:443754/251‑1 (MQ=255)                         cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGAACTCAATGCCTGATGCGACGCGTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                                         >  1:218632/1‑251 (MQ=255)                                aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAAACATCGCAATAGCCTTCTTTATCGTGATTTTCTATGTCGGATAAGGC‑‑GATA‑‑ACGC‑C‑G‑AATCCGACATCATTGCGTGATGCTTCGCTTGCCTCCTCTTCTCAAGCCTCTCTTAACCGTTGGATAATTTTTCAAGCTCTTCACACCCg                                                                                                                                                  >  2:128084/1‑251 (MQ=14)                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGCCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATACATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCATCTCTTAACCGTTGGTTAATTTTTCTAACTCTTCACAcctagtctca                                                                                                                                            >  2:486617/1‑242 (MQ=14)                                             gCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCtt                                                                                                                                     <  1:172146/251‑1 (MQ=255)                                                  aTTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCAt                                                                                                                                 >  1:108394/1‑250 (MQ=255)                                                          ggtcggttagtgCGTTTACGCCCCATTCTACAACAACGCCTGATGTGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTTGCGTGATTTTGTAGGGCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAAGGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAg                                                                                                                        <  2:343897/240‑1 (MQ=17)                                                          ggTCGGATAAGGCGTTTATGCCGCATCCGACATCCACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCACTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATACGACATCAATGCCTGATTCAACGCTTGCCGCGTCTCATCAGGCATATCTTACCCGTTGGTTAATTTTTCACACTCTTCACCCCCGGTCTCAAACCCTTCCATACCGCTc                                                                                                                          >  2:144023/1‑249 (MQ=14)                                                           gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAgg                                                                                                                       >  1:428714/1‑251 (MQ=255)                                                                aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCaa                                                                                                                  <  1:159088/251‑1 (MQ=255)                                                                aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCa                                                                                                                   >  1:306413/1‑250 (MQ=255)                                                                  aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTCTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGACCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCACACCCTTCCATACAGCTCAGGTTCAAcc                                                                                                                >  2:331465/1‑251 (MQ=255)                                                                  aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTNTCAAACCCTTCCATCCAGCTCCGGTTCAAcc                                                                                                                >  2:260744/1‑251 (MQ=255)                                                                              ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCt                                                                                                    <  1:490168/251‑1 (MQ=255)                                                                               cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGccgt                                                                                                   >  1:464741/1‑249 (MQ=255)                                                                                     cGACATCACCGCCTGATGCGACGCTTAACGCGTCTGATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGGGCAGCGCCTTTTGCtt                                                                                             <  2:401085/251‑1 (MQ=21)                                                                                     cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                             >  1:208906/1‑251 (MQ=255)                                                                                             aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTtcttct                                                                                     <  1:211251/251‑1 (MQ=255)                                                                                                  tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGAt                                                                                 >  2:84425/1‑250 (MQ=255)                                                                                                       cGACGCGTACCGCGTCTGATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAAcc                                                                           <  2:531769/251‑1 (MQ=255)                                                                                                              tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCttcttc                                                                     <  1:425117/250‑1 (MQ=255)                                                                                                               aactcgtttagataagCCTACGCAAGACAGCGCAATAGCCTGCTTTAGCGTGGTTTTGTAGGTCGTATAAGTCTAGTAT‑‑G‑‑C‑C‑G‑GAACTGACGTCAATGCCTGAGGCGACGCTTGCCGCGTCTTGTCAGGCCTATCTTAACCGTGGATTAATTTTTCAAACTCTTCACGCCCGGTATCAAGCCCTTCCATACAGCTCAGTTTCAACCGGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                                                   <  2:145162/236‑1 (MQ=255)                                                                                                               aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCAc                                                                   <  1:216620/251‑1 (MQ=255)                                                                                                                         aTCAGGCCTACGCCAGACAGCGCAATAGCCTGATTGAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTGTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTGTAAGa                                                          <  2:164142/250‑1 (MQ=255)                                                                                                                              gCCTACGCCAGCCAGCGCATTTGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTTCCG‑‑C‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCGTGCTGCGTCTTATCAGGGCTATCTTAACCGTTGGTTAATGTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCATGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCGTCTCGATAAACCCTTCTTCCCCGTTTAAGAAcatca                                                    <  1:510344/251‑1 (MQ=25)                                                                                                                              gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAAcatc                                                     >  1:398821/1‑250 (MQ=255)                                                                                                                                  aCGCCAGACAGCGCAATAGCCTGATTTAGCGTGCTTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTAACACCCGGTATCAAACCCTTCCATACAGATCAGGTTCAACCAGCGCAGCGCCTTTTGCTTGTTCTTCTCGTTAAAACCTTCTTCACCGTTTAAGAACATCCTTcc                                                >  2:492646/1‑251 (MQ=25)                                                                                                                                               cAATATCCTGATTTAGCGTGATTTTGTAGATCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACCCCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCg                                   <  2:224637/251‑1 (MQ=255)                                                                                                                                                aaTAGCCTGATTTAGCGTGATTTTTTAGGTCGTAGAAGGC‑‑GTTT‑‑A‑‑CGC‑GGCAGCTGACATGAAGGCCTGATGCGAGGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAATCCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCTCCTTTTGCTTAGTCGTCTCGATAAAACCTTCTTCACCGGTTAGGAACATCATTCCCGCCCAggacgcgg                                  <  2:213815/251‑8 (MQ=255)                                                                                                                                                        gATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAcc                           <  2:208906/250‑1 (MQ=255)                                                                                                                                                        gATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAAc                            <  1:451259/249‑1 (MQ=255)                                                                                                                                                           ttAGGGTGATTTTGTAGGTCTGATGAGGC‑‑GTTT‑‑A‑‑CGCTG‑CATCCGTCATCAATGCCTGATTGGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAGTTTTCAATCTTTTCACGCCCGGTATCAACCCTTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTTTTCTTCTCGATAAAACCTTCTTCACTGTTTAAGATCATCATTCCCGCCCAGTACTCGGAATAACCGGt                        <  2:95197/250‑1 (MQ=255)                                                                                                                                                            tAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGtt                       <  1:105987/250‑1 (MQ=255)                                                                                                                                                             aGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCTCCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTc                      <  1:464462/250‑1 (MQ=255)                                                                                                                                                             aGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCGTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTc                      <  1:250898/250‑1 (MQ=255)                                                                                                                                                                     tttGTAGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAAtt              <  2:398142/250‑1 (MQ=255)                                                                                                                                                                        gtagggctgATAAGGC‑‑GTTT‑‑ACTC‑C‑G‑GATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACGCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGTTTCAACCAGTGCATCGCCTTTTGCTTGTTCTTCTCGATAAAACCTGCTTCACCGTTTGAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCag          <  1:300170/243‑1 (MQ=38)                                                                                                                                                                          aGGTCGGATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCag          >  2:383286/1‑249 (MQ=255)                                                                                                                                                                                gATAAGGC‑‑GTTT‑‑ACGC‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGt  <  1:104643/251‑1 (MQ=255)                                                                                                                                                                                                 ||                                                                                                                                                                                                                                                GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGC‑‑GTTT‑‑A‑‑C‑C‑G‑CATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTGCAGAGCTGCGT  >  NC_000913/4296193‑4296616

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