Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,325,152	A→G	I5T (ATC→ACC)	phnC ←	phosphonate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,325,152	0	A	G	95.2%	68.1 / ‑6.4	21	I5T (ATC→ACC)	phnC	phosphonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (7/13);  minor base T (1/0);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 3.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAG  >  NC_000913/4325021‑4325279                                                                                                                                    |                                                                                                                                gACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGt                                                                                                                          <  2:28939/139‑1 (MQ=255)              ggATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAg                                                                                                             <  1:106524/139‑1 (MQ=255)                aTTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGcc                                                                                                           <  1:200080/139‑1 (MQ=255)                         aaCCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTaaaa                                                                                                  <  2:702102/139‑1 (MQ=255)                               aCGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGTTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTgg                                                                                            >  1:513021/1‑139 (MQ=255)                                             gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                                          >  2:209792/1‑111 (MQ=255)                                             gCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCt                                                                                                          <  1:209792/111‑1 (MQ=255)                                                 ccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCat                                                                                                                    >  1:450133/1‑97 (MQ=255)                                                 ccaTTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCat                                                                                                                    <  2:450133/97‑1 (MQ=255)                                                          cGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGa                                                                 >  1:484563/1‑139 (MQ=255)                                                           gTGATGAATGTTCAGATCAACCGCATGCAGCGCCTCATGCTGATTGAAGGTTTTGGCGAGCTTATCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGAt                                                                <  2:484563/139‑1 (MQ=255)                                                                         gATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTGTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAATAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTaaaa                                                  <  1:651748/139‑1 (MQ=255)                                                                          aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAg                                                 <  1:418492/139‑1 (MQ=255)                                                                          aTCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAg                                                 >  2:269423/1‑139 (MQ=255)                                                                              aCCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGAc                                             >  2:686712/1‑139 (MQ=255)                                                                               ccGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACt                                            >  2:235381/1‑139 (MQ=255)                                                                                 gCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTaa                                          <  1:235381/139‑1 (MQ=255)                                                                                                    gATTGAAGGTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACggg                       <  2:343861/139‑1 (MQ=255)                                                                                                            gTTTTGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGAt               <  2:513021/139‑1 (MQ=255)                                                                                                                tGGCGAGCTTCTCGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACa           <  1:269423/139‑1 (MQ=255)                                                                                                                         tctcGACACGGGTAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAg  >  1:105279/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                GACGTAAAAGGGTGGATTTTCCGGAACCCGACGGCCCAAGCAGAGCCACCATTTCACCGTGATGAATGTTCAGATCAACCGCATGCAGCGCCTGATGCTGATTGAAGGTTTTGGCGAGCTTCTCGACACGGATAATCGTTTGCATGATGCAGCCTCCCTAAAAAAGTGGCTCCATCGTGGCGGATTATTGTGACGATTTGGTTAATTAAAAGTGACTAACAGATGAAGAGTTAACGGGAATTCGATGACAGTGCGGGAG  >  NC_000913/4325021‑4325279

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