Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,273,362	T→A	24.4%	A170A (GCA→GCT)	uvrA ←	excision nuclease subunit A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,273,362	0	T	A	24.4%	20.6 / 4.6	16	A170A (GCA→GCT)	uvrA	excision nuclease subunit A
Reads supporting (aligned to +/- strand):  ref base T (8/4);  new base A (2/2);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

AAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGACGTC  >  NC_000913/4273221‑4273502                                                                                                                                              |                                                                                                                                             aaGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTg                                                                                                                                             >  1:533437/1‑143 (MQ=255) aaGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTCTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATCCGAGCACGGATGTAGCCCGGGCTTg                                                                                                                                             >  2:224397/1‑143 (MQ=255)       ttGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAggg                                                                                                                                       >  1:585176/1‑142 (MQ=255)        tGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                      <  1:194411/143‑3 (MQ=255)        tGGGGAAGATCGTCACGCACCTTGGAGCGATCAACTACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGGTTGGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCAAGCCAgggc                                                                                                                                      <  1:286862/143‑3 (MQ=255)                   gTCACGCACCTTGAAGCGATCAACCACCACTTAAATGGTATGTTTCTTTTGAAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGtat                                                                                                                           >  2:351038/1‑141 (MQ=255)                                                 ttCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCTGGTCCTCCAGCGTTTTGGTGTGTTCGCCTTGGCGCTCTTCAATGa                                                                                             >  2:438508/1‑143 (MQ=255)                                                               tCTTTTACCGTTCCCTTTTTGTCGGTTCGGAATGTTCGCAGTCTTCGCCATCTATACGAGCACGGATGTGGCCCTTGCTTGCCAGGTTTTTCAGCGTTTTGGTGTGTTCGCTTTTGCGCTCTTTAATGATTGGCGCGAGTAGc                                                                               <  2:209648/143‑1 (MQ=255)                                                                    tGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTCTAATGATTGGCGCGAGTAGCATCAg                                                                          <  1:342661/143‑1 (MQ=255)                                                                                 tgagaggAGCGGAAAGACCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTcc                                                             <  1:43063/138‑1 (MQ=255)                                                                                    gCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGc                                                          >  2:245897/1‑143 (MQ=255)                                                                                         tCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGTCTAGCCAGGTTCTCCATCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTGGCTTCAGACGCTTGCCTTCCGCCAGCGa                                                     >  2:336389/1‑143 (MQ=255)                                                                                                                 tCAATACGAGCACGGATGTAGCCCTGGCTTGCCCGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTCTCCTGCTTGGCGCGAGTAGCTTCCCACGCTTGCCTTCCGGCTGCGTCCGCTCGTCGTCCACCGTcagcct                             >  1:78403/1‑138 (MQ=255)                                                                                                                           cACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGCTTGTCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGCTATCCCCCACCTCGCTGCCgctccgc                   >  2:426083/1‑137 (MQ=255)                                                                                                                                  gagaaccTGGCAAGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGAATGCGGCGGGGAGCACCAGACGCTTGCCTTCGGGCTGCGACAGCACGTTACCCACCATCTGGCTCACGGTTCGCGCCGCCa            >  2:110943/6‑143 (MQ=255)                                                                                                                                            cTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATTAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGAcgtc  <  1:351038/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             AAGACGTTGGGTAAGATCGTCACGCACCTTGAAGCGATCAACCACCACTTCAATGGTATGTTTCTTTTGCAGTTCCAGTTTTGGCGGATCGGAAAGATCGCAGACTTCGCCATCAATACGAGCACGGATGTAGCCCTGGCTTGCCAGGTTCTCCAGCGTTTTGGTGTGTTCGCCTTTGCGCTCTTTAATGATTGGCGCGAGTAGCATCAGACGCTTGCCTTCCGGCTGCGACAGCACGTTATCCACCATCTGGCTGACGGTTTGCGCCGCCAGCGGGACGTC  >  NC_000913/4273221‑4273502

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