Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,035	A→T	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,035	0	A	T	96.8%	57.1 / NA	32	I485M (ATA→TTG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base T (23/7);  minor base G (1/0);  total (25/7)
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

GGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATG  >  NC_000913/730901‑731057                                                                                                                                       |                       ggACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCATCACAGCCAATTAGCGTTAGCCACCGGGCCTGACGGGCTGGAATTgcgc                   >  2:683387/1‑140 (MQ=11) ggACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTgcgc                   <  1:1722762/140‑1 (MQ=21) ggACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTgcgc                   >  2:1470736/1‑140 (MQ=255) gcaaaacaGAATACGGTCCGGATGTGGTGACGGACCTCATTACGCTCATCACCACGCCGCATCGCAGGGCATCGGCGTTTTACTATGAGCACCACAGACAGATAACGTCAGCCGCCGGGCCTGACGGGCTGGAAATgcgc                   <  2:2084887/136‑1 (MQ=14)  gACAACAGAATACAGGCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGcc                  <  1:1870367/140‑1 (MQ=12)  gACAACAGAATACAGATCGGATGTGATGACGGGCCTCGTCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGcc                  <  1:1038698/140‑1 (MQ=11)   acaacaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCg                 <  1:2207852/140‑1 (MQ=21)      acaGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGa              <  1:2462066/140‑1 (MQ=16)       caGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGaa             >  1:1755507/1‑140 (MQ=255)       caGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGaa             >  2:1538397/1‑140 (MQ=21)        aGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAat            >  2:1664394/1‑140 (MQ=255)          aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGcc                  <  1:139822/132‑1 (MQ=16)          aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGcc                  >  2:139822/1‑132 (MQ=21)          aaTACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGACACCGGGGCTGACGGGCTGGAGTTGCGCCGGGAatat          >  1:719383/1‑140 (MQ=11)            tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGGCACGGGGGCTGACGGGCTGGAATTGCGCCGGGAATatga        >  2:780596/1‑140 (MQ=11)            tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCgg                >  2:1447916/1‑132 (MQ=18)            tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCgg                <  1:1447916/132‑1 (MQ=2)            tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatga        >  1:2082883/1‑140 (MQ=18)            tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTCAACGTCAGGCACCGAGGCTGACGGGCTGGACGTGCGCCGGGAATatga        >  2:1200850/1‑140 (MQ=11)            tACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCAGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatga        >  1:1844535/1‑140 (MQ=16)             aCAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatgat       >  1:2348297/1‑140 (MQ=18)              cAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATatgatg      >  2:2422362/1‑140 (MQ=18)               aGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCTGGCCTGACGGGCTGGAATTGCGCCGGGAATAtgatga     >  2:1792005/1‑140 (MQ=12)               aGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATAtgatga     >  2:2426628/1‑140 (MQ=17)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt   >  1:62561/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt   >  1:38132/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt   >  2:73176/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt   >  1:247/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt   >  2:1519718/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAATTGCGCCGGGAATATGATGAAt   >  1:2024729/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGCCCTGACGAGCTTGAATTGCGCTGGGAATATGATGAAt   >  1:876250/1‑140 (MQ=11)                 tCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGACAGCCACCGGGCCTGATCGGCTGGAGTTGCGGCGGGAATGTGATGAAt   >  2:1184206/1‑140 (MQ=11)                  ccGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCAGACGGGCTGGAAATGCGCAGGGAATATGATGAAtt  >  2:1079138/1‑139 (MQ=255)                  ccGGATGCGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTGAGCCACCGGGCCAGACGGGCTGGAAATGCGACGGGAATAGGATGAAtt  >  2:785364/1‑139 (MQ=11)                                                                                                                                       |                       GGACAACAGAATACAGTCCGGATGTGGTGACGGGCCTCATCACGCGCATCACCACGCCGGATGGCAGGGCATCGGCGTTTTACTATAACCACCACAGCCAGTTAACGTCAGCCACCGGGCCTGACGGGCTGGAAATACGCCGGGAATATGATGAATG  >  NC_000913/730901‑731057

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