Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	730,389	A→T	16.0%	E269D (GAA→GAT)	rhsC →	Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	730,389	0	A	T	16.0%	62.6 / 3.9	30	E269D (GAA→GAT)	rhsC	Rhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base A (10/15);  new base T (2/3);  total (12/18)
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

CGTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTCCTG  >  NC_000913/730290‑730405                                                                                                    |                 cGTGACGGATGGTGCCGGACGTCACTACCGGCTGGTACTGACGACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAAGCCCTGTCCGGCGGGACGGaa                  >  1:3487006/1‑100 (MQ=11)  gTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAAc                 <  1:581246/100‑1 (MQ=33)   tGACGGATGGTGCCGGGCGTCACTGCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAAcc                <  2:362611/100‑1 (MQ=17)    gACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAAcc                <  1:4292943/99‑1 (MQ=33)    gACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCg               >  1:3290277/1‑100 (MQ=33)    gACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCg               <  1:3500276/100‑1 (MQ=33)     aCGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCg               >  1:340834/1‑99 (MQ=33)      cGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTc             <  2:3364149/100‑1 (MQ=33)      cGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTc             <  2:1941365/100‑1 (MQ=33)       ggATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTcc            <  1:751983/100‑1 (MQ=33)        gATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGTCGCAGCGGGCAGCAGACTCCCGGCAGCAGGCCATTTCCGGCGCGACCGATCCGTc             >  1:1264156/1‑98 (MQ=11)         aTGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGc          <  2:619507/100‑1 (MQ=33)         aTGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGc          <  2:1690901/100‑1 (MQ=33)          tGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGc          <  1:1656920/99‑1 (MQ=32)           ggTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgcc        <  1:744893/100‑5 (MQ=11)           ggTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCtt        <  2:3676315/100‑1 (MQ=32)           ggTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCtt        >  2:3760450/1‑100 (MQ=33)           ggTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGC‑GGACGGTTCCGTcccgccg       <  1:2522110/100‑6 (MQ=11)            gTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCttt       >  1:3524408/1‑100 (MQ=33)            gTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCttt       >  2:468820/1‑100 (MQ=33)             tGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCtttt      >  1:3867199/1‑100 (MQ=32)             tGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCtttt      >  2:2802953/1‑100 (MQ=32)              gCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTc     >  2:1485135/1‑100 (MQ=32)              gCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTc     >  2:4225495/1‑100 (MQ=32)               ccGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTcc    >  1:3127531/1‑100 (MQ=32)               ccGGGCGTCACTTCCGGCTGCTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTcc    <  1:4101099/100‑1 (MQ=17)               ccGGGCGTCACTGCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTcc    <  1:137289/100‑1 (MQ=17)                 gggCTTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTCCTg  <  1:4246105/100‑1 (MQ=11)                 gggCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGTTCCGTcccgccggaaa   <  2:2500143/99‑10 (MQ=11)                 gggCGTCACTGCCGGCGGGTACTGACCGCGCAGGCGCAGCTGGCAGAAGAGGCTCTGCAGCAGGCCATTTCCGGCGGGACGTATCCGTCCGCTTTTCCTg  >  1:2570477/1‑100 (MQ=11)                 gggCGTCACGTCCGGCTGGTACTGACCACGCAGGCGCAGCGTGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGATGGCACCGTCCGCTTTTCGTg  <  1:236390/100‑1 (MQ=11)                                                                                                    |                 CGTGACGGATGGTGCCGGGCGTCACTTCCGGCTGGTACTGACCACGCAGGCGCAGCGGGCAGAAGAAGCCCGGCAGCAGGCCATTTCCGGCGGGACGGAACCGTCCGCTTTTCCTG  >  NC_000913/730290‑730405

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