Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,017,703	T→C	N1175D (AAC→GAC)	rhsB ←	rhsB element core protein RshB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,017,703	0	T	C	100.0%	230.3 / NA	71	N1175D (AAC→GAC)	rhsB	rhsB element core protein RshB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/37);  total (34/37)

TGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGT  >  W3110S.gb/4017645‑4017748                                                           |                                              tgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGc                                             >  1:212661/1‑61 (MQ=25) tgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGc                                             >  1:1174889/1‑61 (MQ=25) tgcCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGc                                             >  1:4643548/1‑61 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCa                                         <  1:3948322/61‑1 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCa                                         <  1:1479311/61‑1 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCa                                         <  1:3525845/61‑1 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCa                                         <  1:240762/61‑1 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCa                                         <  1:2275567/61‑1 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCa                                         <  1:192025/61‑1 (MQ=25)     ccGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGGCGCCCCa                                         <  1:655568/61‑1 (MQ=255)      cGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCAtt                                       <  1:32165/62‑1 (MQ=25)      cGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCAtt                                       <  1:3622624/62‑1 (MQ=25)      cGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCAtt                                       <  1:4224569/62‑1 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:2131136/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:2048028/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:3255024/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:317472/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:3349067/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:3498359/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:3025263/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:3644333/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:3649764/1‑62 (MQ=25)           gCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATtcatca                                  >  1:4251553/1‑62 (MQ=25)             cGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATtata                                <  1:4449173/62‑1 (MQ=255)              ggATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                <  1:3994612/61‑1 (MQ=25)                 tAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTc                             <  1:4621376/61‑1 (MQ=25)                      tgctgcAGCTTATGCGGGTTCTCTTCATTCAGCAGTTCGCCCCATTCATCATATTCTGCGca                       <  1:3350207/62‑1 (MQ=11)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:4367154/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:3173799/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:3580165/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:3097599/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:3758268/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:4264108/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:1142372/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:1138122/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:753263/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:1358046/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:1402540/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:1959450/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:2008036/62‑1 (MQ=25)                      tgctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:3077245/62‑1 (MQ=25)                       gctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCTCCCCATTCATCATATTCTGCGcac                      <  1:176383/62‑1 (MQ=255)                       gctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGcac                      <  1:2081275/62‑1 (MQ=25)                       gctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:2145099/61‑1 (MQ=25)                       gctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:773935/61‑1 (MQ=25)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCGCATTCATGata                                >  1:1612742/1‑51 (MQ=11)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTcatcat                                 >  1:4117320/1‑50 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:4551717/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:4687076/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:931393/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:3930076/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:3701879/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:1617618/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:3247244/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                >  1:1743909/1‑51 (MQ=255)                        ctgcAGCTGATGCGGGTTCTCTTCATTCAACAGGTCGCCCCATTCATCATATTCTGCGcac                      <  1:3813244/61‑1 (MQ=255)                         tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGcacca                    >  1:1351830/1‑62 (MQ=25)                         tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGcacca                    >  1:3213004/1‑62 (MQ=25)                         tgcAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGcacca                    >  1:2504327/1‑62 (MQ=25)                           cAGCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCata                                <  1:2185341/48‑1 (MQ=255)                             gCTGATGCGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGca                       <  1:3085571/55‑1 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTg                           >  1:1735097/1‑44 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGccc         >  1:1454082/1‑62 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGccc         >  1:1794877/1‑62 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGccc         >  1:1727821/1‑62 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGcc          >  1:374509/1‑61 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGcc          >  1:4054040/1‑61 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGcc          >  1:332845/1‑61 (MQ=255)                                    cGGGTTCTCTTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGATGccc         >  1:4048666/1‑62 (MQ=255)                                            ctTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGt  <  1:1297372/61‑1 (MQ=255)                                            ctTCATTCAGCAGGTCGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGt  <  1:1211800/61‑1 (MQ=255)                                                           |                                              TGCCCCGGCAGCCGGATAAGCTGCTGCAGCTGATGCGGGTTCTCTTCATTCAGCAGGTTGCCCCATTCATCATATTCTGCGCACCACGCTGTTGCCCCTTCCGT  >  W3110S.gb/4017645‑4017748

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑