Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,783,813	G→A	L106L (CTG→TTG)	secB ←	protein export chaperone SecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,783,813	0	G	A	90.0%	18.2 / ‑3.4	10	L106L (CTG→TTG)	secB	protein export chaperone SecB
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (6/3);  total (6/4)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.37e-01

CCTCGTGGCCATTTCTTGCCAGGGTGATGGCAAGAGCGGTGCCGTACGAGCCGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGAATGTTCGGGCAGTATGCTCCCAGGCAATGCGCCATCTGGGTGCCTTCGATACCCGCGATGGAGAAAATACCGCCCTGCTGAACTTCACACAGGAACGCGGTTTCTTCGCCCAAAGAGGCCGTTACGGTAACACGCAGTACCACTTCGTATACGTCATCTGCCAGTTGGGAAGATGCCGTATCCAGATCAAGTTTAACTTCTGGTTGCCAATCTTTCTGGAAAACGTGCGGCGCGTTCGGCGCTT  >  NC_000913/3783571‑3784058                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      ccTCGTGGCCATTTCTTGCCAGGGTGATGGCAAGGGCGGTGCAGTACTAGCCGGCACCGAGCACAGTCATTGAAGCATTACGTGGGTTCATCAGGCATCCGGATTTTCTTCATTACCCTCGCCATCCTTCTGCGGCAAAGATTTCATGAACAGCGCACCGAAGTTAACCGGCGCAATGTCCAGTTTCGGTAATTTACCGCTGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAAAATGTTCg                                                                                                                                                                                                                                               <  1:235437/251‑1 (MQ=255)                        gTGATGGCAAGAGCGGTGCCGTACGAGCCGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAAAATGTTCGGGCAGTATGCTCCCAGGCAATgc                                                                                                                                                                                                                        <  1:150934/251‑1 (MQ=255)                                                                                            caGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAAAATGTTCgggn                                                                                                                                                                                                                                            <  1:204869/163‑2 (MQ=255)                                                                                                    cTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCCTCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGACTGTACCGAGGGCTACCATGCTGGTGATGCACTCACGAGCATACGGGAACAACATGTTCGGGCAGTATGCTCCCAGGCAAt                                                                                                                                                                                                                          >  1:182885/1‑173 (MQ=255)                                                                                                                                                                              cAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGCGATGCACTCACGAGCCTACAGGAACAAAATGTTCGGGCAGTATGCTCCCAGGCAATGCGCCCTCTGGGTGCCTTCGATACCCGCGCTGGAGAAAATACCGCCCTGCTGAACTTCACACAGGAACGCGGTTTCTTCGCCCAAAGAGGCCGTTACGGTAACACGCAGt                                                                                                            >  1:45567/1‑209 (MQ=255)                                                                                                                                                                              cAAGGTTCAGTTGCGGGAATGGACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGAATGGTCGGGCAGTATGCTCCCAGGCAATGCGCCATCTGGGTGCCTTCGATACCCGCGATGGAGAAAATACCGCGCTGCTGAACTTCACACAGGa                                                                                                                                                        <  1:117285/165‑1 (MQ=255)                                                                                                                                                                                                     ccGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGCACAAAATATTCGGCCAGTATGCTCCCAGGCAATGCGCCATCTGGGTGCCTTCGATACCCGCGATGGAGAAAATACCGCCCTGCTGAACTTCACACAGGAACGCGGTTTCTTCTCCCAAAGAGGCCGTTACGGTAACACGCAGTACCACTTCGTATACGTCATCTGCCAGTTGGGAAGATGCCGTTTCCAGATCAAGTTTAACTTCTgg                                           >  1:45743/1‑251 (MQ=255)                                                                                                                                                                                                      cgcgGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAAAATGTTCGGGCAGTATGCTCCCAGGCAATGCGCCATCTGGGTGCCTTCGATACCCGCGATGGAGAAAATACCGCCCTGCTGAACTTCACAAAGGACCGGGGTTTTTTCGCCCAAAAAGCCCGTTACGGAACCACGCAGAACCACTTCGTATACGTCATCGGCCAGTTGGGAAGATGCGGTATCCAGATCAGGTTTAACTTCTGGt                                          >  1:258134/1‑251 (MQ=255)                                                                                                                                                                                                                        ttaTGCACTCACGAGCATACGGGAACAAAATGTTCGGGCAGTATGCTCCCAGGCAATGCGCCCTCTGGGTGCCTTCGATACCCGCGATGGAGAAAATCCCGCCCTGCTGAACTTc                                                                                                                                                                >  1:301308/3‑115 (MQ=255)                                                                                                                                                                                                                                              gAACAAAATGTTCGGGCAGTATGCTCCCAGGCAATGCGCCATCTGGGTGCCTTCGATACCCGCGATGGAGAAAATACCGCCCTGCTGACCTTCACACAGGAACGCGGTTTCTTCGCCCAAAGAGGCCGTTACGGTAACACGCAGTACCACTTCGTATACGTCATCTGCCCGTTGGGAAGATGCCGTATCCAGATCAAGTTTAACTTCTGGTTGCCAATCTTTCTGGAAAACGTGCGGCGCGTTCGGCGCtt  >  1:70736/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      CCTCGTGGCCATTTCTTGCCAGGGTGATGGCAAGAGCGGTGCCGTACGAGCCGGCACCGATCACAGTCATTGAAGCATTACGTTGGTTCATCAGGCATCCTGATGTTCTTCAGTACCTTCGCCAGCCTGCTGCTGCAAATAGTTCATGAACAGCGCATCGAAGTTAACCGGCGCAAGGTTCAGTTGCGGGAATGTACCGCGGGATACCATGCTGGTGATGCACTCACGAGCATACGGGAACAGAATGTTCGGGCAGTATGCTCCCAGGCAATGCGCCATCTGGGTGCCTTCGATACCCGCGATGGAGAAAATACCGCCCTGCTGAACTTCACACAGGAACGCGGTTTCTTCGCCCAAAGAGGCCGTTACGGTAACACGCAGTACCACTTCGTATACGTCATCTGCCAGTTGGGAAGATGCCGTATCCAGATCAAGTTTAACTTCTGGTTGCCAATCTTTCTGGAAAACGTGCGGCGCGTTCGGCGCTT  >  NC_000913/3783571‑3784058

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