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	94.9%	111.1 / NA	41	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/1);  major base C (17/20);  minor base A (2/0);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 2.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
*	NC_000913	4,296,380	2	.	G	95.1%	109.5 / ‑1.8	41	intergenic (+586/+56)	gltP/yjcO	glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (19/20);  total (20/21)
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

GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTG  >  NC_000913/4296193‑4296605                                                                                                                                                                                               ||                                                                                                                                                                                                                                     gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGACTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCAGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                                                                                                                                                      >  1:491871/1‑251 (MQ=2)       tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTGGCGTGATTTTGTAGGTCGGGTAAGGCGTTGATGCCGCATCCGACATCAACGCCTGATGCGACGCTGAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                                                                                                                                                   <  1:355444/248‑1 (MQ=11)       tATCAGGACTACGCCAGACAGCGCAATGGCCTGATGTAGCGTGATTTTGTAAGTCGGATAAGGCGTTTATGCCGCATCCGACATCAGCGCCTGATTCGACGCTTAACGCTTCTTATCAGGGCTACGCTAGACAGCGCAATTGCCTGATTTAGCGTGATTGTTTGTGTCGGATTTGGCGTGT‑‑ACGC‑‑CGCATCCGACAGTA‑ATGCCGTAGTCGTCGCTTGCCGCGTCTTATCAGCCCTATCTTAACCGTTg                                                                                                                                                                  <  1:48382/249‑1 (MQ=2)         ccgtgCCTACGCCAGACAGCGCAATAGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTTATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                                                                                                                                              <  2:229109/247‑1 (MQ=11)           aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                                                                                                                                            >  2:172284/1‑251 (MQ=255)                   cctatctaGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                     <  2:429721/243‑1 (MQ=255)                       gACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                                                                                                                                                >  1:530984/1‑251 (MQ=255)                        aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                                                                                                                                               <  1:141421/251‑1 (MQ=255)                          agagcAATAGCCTGAGTTAGAGTGATTTTGTAGGTTGCATAAGGCGTTTATGCCGCATCGGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATAAGGCCTACGTCAGCCAGTGGTATCGCTTGCTTTAGCGCGATTTTGTGGGTCGGATAAGGCGTTTCCG‑‑C‑‑CGCCTCTGACA‑CATATGCCTGATGTGACGCTTGTCGCGTCTTATAAGGCCTATCTTAACCGTTGGTTCATTTTTCACACTCTTc                                                                                                                                              >  2:288984/4‑250 (MQ=2)                          aGCGCTATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCTCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGTCTATCTTAGCCGTTGGTTAATTTTTCAACCTCTTca                                                                                                                                             <  1:414032/251‑1 (MQ=11)                          aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                                                                                                                                                     >  1:429721/1‑243 (MQ=255)                          aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                                                                                                                                              >  1:499847/1‑250 (MQ=255)                             gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACCTCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCacac                                                                                                                                          >  1:358453/1‑251 (MQ=255)                                    ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt                                                                                                                                   <  1:172284/251‑1 (MQ=255)                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGCTTTAGCGTGATTTTGTACGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑TTGACTGATGCGACGCTTGCCGCGTCTTTTCAGCCTTATCTTAACCGTTGGATAATTTTTCCAACTCTTCACACCCGGTATCa                                                                                                                                 >  2:57856/1‑250 (MQ=14)                                       gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACCTCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa                                                                                                                                 >  1:13001/1‑250 (MQ=255)                                          ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCCGGCCTACGCCAGACAGCGAAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGATT‑‑AAGC‑‑CGCATCTGACCTCA‑ATGCCAGATGCGACGCTTTCCGCCTCTTATCAGGCCTATCTTAAGCTTCGGTACATTTTTCCAACCCATCACACCCGGTATccagcc                                                                                                                             >  2:248513/1‑246 (MQ=11)                                                 aaaatagtagattggaaggaacGTTTATGCCGCAGCGTACACCAACGCCAGATACAACACATAACTCGTCTTATCACACCTACGCCATTCACCGCAAAATCCAAATATGGCGTGATTTTGTATGTCGTATAAGTCGTTT‑‑A‑‑CGACGTCCCCGACATCA‑ATATCTTAATCGACCCTTGCCGCGTCTTAGCGAGCCTATCATATCCGTTGTTTAATTTTTCAATCTCTTCACACCCGGTATCAAACCCTTCCAt                                                                                                                      <  1:50057/230‑1 (MQ=255)                                                   ttttGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATa                                                                                                                     <  2:146203/250‑1 (MQ=255)                                                              ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                                                                                                        >  2:108660/1‑204 (MQ=255)                                                              ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                                                                                                                                                        <  1:108660/204‑1 (MQ=255)                                                                           aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCCATAGCCTGACTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCAGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAgcg                                                                                            >  1:328038/1‑251 (MQ=18)                                                                                     cGACATCAACGCCTGATGCGACGCTTAACACGTCTTATCAGGCCTACGCCAGACAGCGCAATCGACTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATACGACATCA‑ATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATATTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAACCCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCtt                                                                                  >  1:388495/1‑251 (MQ=18)                                                                                        cATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGtt                                                                               <  1:40632/251‑1 (MQ=255)                                                                                           cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt                                                                            >  2:266898/1‑251 (MQ=255)                                                                                           cAACGCCTGATGCAAAGCTTACCGCGTCTTATCCGGCCTACGCCAAACAGCGCAATAGCCTGATTTAGCGAGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGAATCCGACATCC‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCACGCCTATCTTACCCGTTGGTTACTTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGttctt                                                                            >  2:369721/1‑251 (MQ=21)                                                                                               gccggATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                        <  2:311104/247‑1 (MQ=255)                                                                                               gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCTCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCg                                                                        >  2:169473/1‑251 (MQ=255)                                                                                                 cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCAAGACCACGCAATAGCCTGATTTAGCGAGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGAt                                                                      >  1:33408/1‑251 (MQ=25)                                                                                                           gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACcttct                                                            >  2:165735/1‑251 (MQ=255)                                                                                                                   cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCGTGTTCTGCTCGATAAAACCGTCTTCACCGtt                                                    <  2:339093/251‑1 (MQ=25)                                                                                                                        tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTGTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTGCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCGTCACCGTTTaa                                                 <  2:530984/249‑1 (MQ=255)                                                                                                                                  ccGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCGTCACCGTTTAAGAACATCATTcc                                     <  2:349488/250‑1 (MQ=255)                                                                                                                                            gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTAc                           <  2:6497/251‑1 (MQ=255)                                                                                                                                            gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTAc                           <  1:380271/251‑1 (MQ=255)                                                                                                                                                aaTAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                                                                                                               >  1:398405/1‑163 (MQ=255)                                                                                                                                                aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTc                                                                                                               <  2:398405/163‑1 (MQ=255)                                                                                                                                                           ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCTGAATAACCGGt             <  1:445981/250‑1 (MQ=255)                                                                                                                                                             aGCGTGATTTTGTAGGTCGGCTAAGGCGTTT‑‑AAGC‑‑AGCATCCGCCATCA‑ATGCCTGCTGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGCTAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg          >  2:314089/1‑251 (MQ=255)                                                                                                                                                             aGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGGCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACATCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTTTTCGCGATAAAACCTTCTGCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCg          <  1:506545/251‑1 (MQ=38)                                                                                                                                                                      ttGTAGGTCGGATAAGGCGTTT‑‑ACGC‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTg  <  1:326939/250‑1 (MQ=255)                                                                                                                                                                                               ||                                                                                                                                                                                                                                     GCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑C‑‑CGCATCCGACATCA‑ATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACCCTTCCATACAGCTCAGGTTCAACCAGTGCAGCGCCTTTTGCTTGTTCTTCTCGATAAAACCTTCTTCACCGTTTAAGAACATCATTCCCGCCCAGTACTCGGAATAACCGGTTCGGGAAATTG  >  NC_000913/4296193‑4296605

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: ‑