Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	W3110S.gb	4,174,189	C→T	intergenic (+20/+9)	bfr → / ← gspO	bacterioferritin, iron storage and detoxification protein/bifunctional prepilin leader peptidase and methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,174,189	0	C	T	100.0%	113.0 / NA	40	intergenic (+20/+9)	bfr/gspO	bacterioferritin, iron storage and detoxification protein/bifunctional prepilin leader peptidase and methylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/27);  total (13/27)

GATGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAC  >  W3110S.gb/4174121‑4174255                                                                     |                                                                   gaTGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAAtcc                                                                  <  1:271734/71‑1 (MQ=255)    gggCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAAtcctcc                                                               >  1:771158/1‑71 (MQ=255)    gggCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAAtcctcc                                                               >  1:3969521/1‑71 (MQ=255)    gggCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAAtcctcc                                                               >  1:4502009/1‑71 (MQ=255)      gCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCg                                                             >  1:498211/1‑71 (MQ=255)       cgtGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCgg                                                            <  1:2450428/69‑1 (MQ=255)           cAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTtaa                                                        >  1:3049099/1‑71 (MQ=255)            aaaaTTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTtaat                                                       <  1:642428/71‑1 (MQ=255)              aaTTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTtaataa                                                     <  1:3787822/71‑1 (MQ=255)                 tATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTtaataa                                                     >  1:3584795/1‑68 (MQ=255)                  aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTtaataa                                                     >  1:2518628/1‑67 (MQ=255)                  aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:2260040/70‑1 (MQ=255)                  aTCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAagag                                                  <  1:3534188/70‑1 (MQ=255)                   tCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAgaga                                                 <  1:153478/70‑1 (MQ=255)                       cAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATc                                               <  1:3113151/68‑1 (MQ=255)                         aGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGc                                          >  1:2842188/1‑71 (MQ=255)                           gacaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:3202285/70‑1 (MQ=255)                           cacaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:3084046/71‑1 (MQ=255)                            acaGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAgaga                                                 <  1:4271256/61‑1 (MQ=255)                             caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:789915/71‑1 (MQ=255)                             caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:3048565/71‑1 (MQ=255)                             caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:2785589/71‑1 (MQ=255)                             caGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:2790978/71‑1 (MQ=255)                              aGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:1059865/70‑1 (MQ=255)                                    gcgAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:3048664/62‑1 (MQ=255)                                    gcgAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:2616662/62‑1 (MQ=255)                                    gcgAAGAAGGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCa                                        <  1:4155625/62‑1 (MQ=255)                                         gaagGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGtt                           >  1:1545493/1‑70 (MQ=255)                                           agGTTGAGCCGTGTCGGGATGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:3273112/71‑1 (MQ=255)                                           agGTTGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTAt                        <  1:585293/71‑1 (MQ=255)                                               tGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCAtg                    >  1:4611006/1‑71 (MQ=255)                                               tGAGCCGTGTCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCAtg                    >  1:4627477/1‑71 (MQ=255)                                                      tgtCGGGAGGGAGAATCCTCCCGGTTAATAAGAGACCCAGGCCa                                        <  1:4259644/44‑1 (MQ=255)                                                        tCGGGAGGGAGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAgc           >  1:2436626/1‑71 (MQ=255)                                                          gggagggaGAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAgcg          >  1:4160546/1‑70 (MQ=255)                                                              gggagAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGc     <  1:780428/71‑1 (MQ=255)                                                              gggagAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGc     <  1:1066908/71‑1 (MQ=255)                                                              gggagAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGc     <  1:4047865/71‑1 (MQ=255)                                                              gggagAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGc     <  1:1844946/71‑1 (MQ=255)                                                               ggagAATCCTCCCGGTTAATAAGAGATCCAGGCCAGa                                      <  1:4313524/37‑1 (MQ=255)                                                               ggagAATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCa        <  1:2529532/67‑1 (MQ=255)                                                                 agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCa   >  1:2583608/1‑70 (MQ=255)                                                                 agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCa   >  1:2950390/1‑70 (MQ=255)                                                                 agaATCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAc  >  1:1089451/1‑71 (MQ=255)                                                                     |                                                                   GATGGGCCTGCAAAATTATCTGCAAGCACAGATCCGCGAAGAAGGTTGAGCCGTGTCGGGAGGGAGAACCCTCCCGGTTAATAAGAGATCCAGGCCAGATACCCCAGGTTTATCATGCTGCCGAGCGCAAGCCAC  >  W3110S.gb/4174121‑4174255

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