Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,037	A→G	I485M (ATA→TTG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,037	0	A	G	91.3%	28.8 / ‑1.3	23	I485M (ATA→TTG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (8/13);  total (9/14)
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

ACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTA  >  NC_000913/730903‑731111                                                                                                                                       |                                                                           acaacaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCg                                                                       <  2:827059/140‑1 (MQ=16)    acaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGa                                                                    <  2:2374142/140‑1 (MQ=16)    acaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGa                                                                    <  1:2468999/140‑1 (MQ=16)    acaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGa                                                                    <  1:1098356/140‑1 (MQ=16)      aGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAat                                                                  >  2:548659/1‑140 (MQ=255)        aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                                                                >  1:1400621/1‑140 (MQ=21)        aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                                                                >  1:1998143/1‑140 (MQ=21)        aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                                                                >  1:2370078/1‑140 (MQ=21)        aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCG‑TTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  1:2139495/140‑1 (MQ=255)         aTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  2:947931/140‑1 (MQ=9)         aTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  2:1998143/140‑1 (MQ=9)         aTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  2:1604555/140‑1 (MQ=9)         aTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  2:1558236/140‑1 (MQ=9)         aTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  2:1372555/140‑1 (MQ=9)         aTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatg                                                               <  2:1155069/140‑1 (MQ=9)          tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                                                                <  2:755478/138‑1 (MQ=2)          tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAatat                                                                >  1:755478/1‑138 (MQ=18)               tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCg                                                                       >  2:125539/1‑126 (MQ=11)               tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCg                                                                       <  1:125539/126‑1 (MQ=0)               tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt                                                         >  2:2467216/1‑140 (MQ=11)               tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt                                                         >  1:828751/1‑140 (MQ=11)               tCCGGATATGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAAt                                                         >  2:1417457/1‑140 (MQ=255)                                                                      tCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTa  <  1:169195/140‑1 (MQ=40)                                                                                                                                       |                                                                           ACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATGGGGCCGTCTGATTCAGGAAACTGCCCCTGACGGCGATATCACCCGCTACCGTTA  >  NC_000913/730903‑731111

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