Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,549,754	G→A	25.9%	intergenic (+43/‑24)	fimC → / → fimD	chaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,549,754	0	G	A	25.9%	42.9 / 14.7	27	intergenic (+43/‑24)	fimC/fimD	chaperone, periplasmic/outer membrane usher protein, type 1 fimbrial synthesis
Reads supporting (aligned to +/- strand):  ref base G (10/10);  new base A (3/4);  total (13/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 = 7.92e-01

ATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATG  >  W3110S.gb/4549691‑4549821                                                                |                                                                    aTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTg                                                                 <  1:2690014/68‑1 (MQ=255)    aCGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGAtt                                                              >  1:2596884/1‑68 (MQ=255)            aaTGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCTGGGTGATTTTAAGCCgg                                                     <  1:3206621/69‑1 (MQ=255)               ggAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCggagg                                                  <  1:2624580/69‑1 (MQ=255)                   tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaata                                              >  1:2326392/1‑69 (MQ=255)                   tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGaata                                              >  1:2400416/1‑69 (MQ=255)                     aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGATTGCCGGGTGATTTTAAGCCGGAGGaataat                                            >  1:2106641/1‑69 (MQ=255)                     aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGATTGCCGGGTGATTTTAAGCCGGAGGaataa                                             >  1:1325500/1‑68 (MQ=255)                       gCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                          <  1:3472416/69‑1 (MQ=255)                       gCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                          <  1:358202/69‑1 (MQ=255)                                 tttttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGATGGAGGAATAATGTCATATCTGaa                                <  1:1735797/69‑1 (MQ=255)                                   tttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAAt                               >  1:985257/1‑68 (MQ=255)                                   tttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAAt                               >  1:2984713/1‑68 (MQ=255)                                   tttCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAAt                               >  1:443279/1‑68 (MQ=255)                                         cTGAATAAAAAGAATTGATTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGAc                        <  1:3272264/69‑1 (MQ=255)                                         cTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGa                         <  1:1601793/68‑1 (MQ=255)                                            aaTAAAAAGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACttt                     >  1:3432262/1‑69 (MQ=255)                                              tAAAAAGAATTGACTGCTGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACttt                     <  1:2747717/67‑1 (MQ=255)                                              tAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGTAGGAATAATGTc                                         <  1:501495/47‑1 (MQ=39)                                              tAAAAAGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGTc                                         >  1:2869744/1‑47 (MQ=255)                                               aaaaaGAATTGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTAcc                  <  1:813075/69‑1 (MQ=255)                                               aaaaaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                          >  1:2941259/1‑45 (MQ=255)                                               aaaaaGAATTGACTGCCAGGTGATTTTAAGCCGGAGGAATAATGt                                          >  1:1475581/1‑45 (MQ=255)                                                        tGACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca         <  1:260072/69‑1 (MQ=255)                                                         gACTGCCGGGTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAaca         >  1:3100604/1‑68 (MQ=255)                                                           cTGCCGGGTGATATTAAGCTGGATGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAacaca       <  1:393310/68‑1 (MQ=255)                                                             gCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACaa    <  1:1947349/69‑1 (MQ=255)                                                                gggTGATTTTAAGCTGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATg  <  1:2497964/68‑1 (MQ=255)                                                                |                                                                    ATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATG  >  W3110S.gb/4549691‑4549821

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: ‑