Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,092,566	(G)6→7	coding (689/951 nt)	gshB →	glutathione synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,092,560	1	.	G	95.0%	65.4 / ‑3.1	20	coding (683/951 nt)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (13/6);  total (14/6)
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

GACTGAGCATGGCACTCGCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCA‑GG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTG  >  NC_000913/3092438‑3092685                                                                                                                            |                                                                                                                               gactgaGCATGGCACTCGCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCA‑gg‑gggggcga                                                                                                                     >  1:581471/1‑129 (MQ=255)    tgaGCATGGCACTCGCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGTTACCGTACTGCCAGGCGCGTATTCCTCA‑GGTGGGGCGAAAc                                                                                                                   >  2:249801/1‑133 (MQ=255)      aGCATGGCACTCGCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCA‑gg‑gggggcgaa                                                                                                                    <  2:56913/129‑6 (MQ=255)      aGCATGGCACTCGCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCA‑gg‑gggggcgaa                                                                                                                    >  1:56913/1‑124 (MQ=255)                  gCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTgg                                                                                                     >  2:590059/1‑133 (MQ=255)                  gCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTgg                                                                                                     >  2:477565/1‑133 (MQ=255)                        gCATGGCGAAAAATTGCCGGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCg                                                                                               <  1:249801/133‑1 (MQ=255)                                 aaaaTTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTcgcg                                                                                      >  2:632238/1‑133 (MQ=255)                                   aaTTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGt                                                                                    <  2:573325/133‑1 (MQ=255)                                             cagccaTTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTc                                                                          <  1:110982/133‑1 (MQ=255)                                                           ggCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGt                                                            >  2:150484/1‑133 (MQ=255)                                                           ggCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGt                                                            >  2:330117/1‑133 (MQ=255)                                                                   aCGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGGGGTCGCGGGGAAGCACATCCGCTGAAGGAAAGTGACTGGaa                                                    >  2:311007/1‑133 (MQ=255)                                                                        tGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCg                                               <  1:590059/133‑1 (MQ=255)                                                                             gtggtggATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGGCTGAAAGTGAGTGGAAGATCGCCCGt                                          >  2:290144/1‑133 (MQ=255)                                                                                       gCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGACTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGc                                >  2:319367/1‑133 (MQ=255)                                                                                           gCCGGTACCGTACTGCCTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGTCTGGAAACTCGGCCGTCAGATAGGGCCGAc                            >  1:210618/1‑133 (MQ=255)                                                                                                          ccTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAgg             >  2:521181/1‑133 (MQ=255)                                                                                                          ccTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAgg             <  1:150484/133‑1 (MQ=255)                                                                                                          ccTGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAgg             >  1:139102/1‑133 (MQ=255)                                                                                                            tGGCGCGTATTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGACAATCGCCCGTCACATCTGGCCGACGCTGAAAGAAAAAGGGc           >  2:254999/1‑133 (MQ=255)                                                                                                                    aTTCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGAttttt   >  1:337821/1‑133 (MQ=255)                                                                                                                     ttCCGCAGGG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTttg  <  1:632238/133‑1 (MQ=255)                                                                                                                            |                                                                                                                               GACTGAGCATGGCACTCGCTACTGCATGGCGCAAAATTACCTGCCAGCCATTAAAGATGGCGACAAACGCGTGCTGGTGGTGGATGGCGAGCCGGTACCGTACTGCCTGGCGCGTATTCCGCA‑GG‑GGGGCGAAACCCGTGGCAATCTGGCTGCCGGTGGTCGCGGTGAACCTCGTCCGCTGACGGAAAGTGACTGGAAAATCGCCCGTCAGATCGGGCCGACGCTGAAAGAAAAAGGGCTGATTTTTG  >  NC_000913/3092438‑3092685

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