Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,703,000	G→T	16.6%	G322G (GGC→GGA)	lepB ←	leader peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,703,000	0	G	T	16.6%	138.9 / 10.3	48	G322G (GGC→GGA)	lepB	leader peptidase
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base T (5/3);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 7.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

TATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATC  >  W3110S.gb/2702932‑2703063                                                                     |                                                                tATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgc                                                                <  1:2763452/70‑1 (MQ=255)          aataatTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGc                                                       >  1:667947/1‑70 (MQ=255)                     ttGTCGCCATAACGGCGACAACGTGAACGAAGATTGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCa                                           <  1:4047631/71‑1 (MQ=255)                     ttGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                            <  1:3142452/54‑1 (MQ=255)                             aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGa                                         >  1:2536802/1‑65 (MQ=255)                             aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTc                                   >  1:4418993/1‑71 (MQ=255)                             aTAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTc                                   >  1:3418548/1‑71 (MQ=255)                               aaCGGCGACAACGTGAACGAAGATGGCTATTAATGGATgccgcc                                                            <  1:1812664/44‑1 (MQ=255)                                aCGGCGACAACGTGAACGAAGATGGCTATTAATGGATTCCGCCaa                                                          >  1:1233083/1‑45 (MQ=255)                                aCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                 <  1:3952712/70‑1 (MQ=255)                                aCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                >  1:11776/1‑71 (MQ=255)                                aCGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGc                                >  1:2465754/1‑71 (MQ=255)                                    cGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAtt                            <  1:880082/71‑1 (MQ=255)                                       cAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                           >  1:2668519/1‑69 (MQ=255)                                       cAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGc                         >  1:3846298/1‑71 (MQ=255)                                        aaCGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGCCTTAAGCGCAGACCAGTCGGCCATTCGc                         >  1:764359/1‑70 (MQ=255)                                        aaCGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCa                                      >  1:2085770/1‑57 (MQ=255)                                        aaCGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCa                                      >  1:3232275/1‑57 (MQ=255)                                          cGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCa                              <  1:1493228/63‑1 (MQ=255)                                            tGAACGAAGATGTCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAg                                          <  1:1291245/49‑1 (MQ=255)                                            tGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                    >  1:4326673/1‑71 (MQ=255)                                             gAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                   >  1:66875/1‑71 (MQ=255)                                             gAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                   >  1:3581472/1‑71 (MQ=255)                                             gAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttct                    >  1:184004/1‑70 (MQ=255)                                              aaCGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcg                                             <  1:3705812/44‑1 (MQ=255)                                              aaCGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                   <  1:1424076/70‑1 (MQ=255)                                              aaCGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                   <  1:2579504/70‑1 (MQ=255)                                              aaCGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCcttctt                   <  1:1549500/70‑1 (MQ=255)                                                cGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAgcg                                             <  1:3414520/42‑1 (MQ=255)                                                cGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCt                >  1:1116730/1‑71 (MQ=255)                                                  aaGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTa              >  1:3244021/1‑71 (MQ=255)                                                   aGATGGCTATTAATGGATTCCGCCAATGCGACTTAAGCGCAGa                                         >  1:2564345/1‑43 (MQ=255)                                                   aGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTAt             <  1:3482084/71‑1 (MQ=255)                                                   aGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTAt             <  1:1753335/71‑1 (MQ=255)                                                   aGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTAt             <  1:1758871/71‑1 (MQ=255)                                                     aTGGCTATTAATGGATTCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCAATCGCCTTCTTGCTTATCg           <  1:2290994/71‑1 (MQ=255)                                                       ggCTATTAATGGATTCCGCCAATGCGACTTAAGCGCAGa                                         >  1:1821358/1‑39 (MQ=255)                                                        gCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa         <  1:683365/70‑1 (MQ=255)                                                        gCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGaa         <  1:2635421/70‑1 (MQ=255)                                                         cTATTAATGGATTCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCctt                      <  1:4369707/56‑1 (MQ=255)                                                           aTTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                 <  1:2501134/43‑1 (MQ=255)                                                            ttAATGGATTCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTc     <  1:1620671/70‑1 (MQ=255)                                                            ttAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt      <  1:3455877/69‑1 (MQ=255)                                                             tAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCt      <  1:2038372/68‑1 (MQ=255)                                                             aaaTGGATGCCGCCAATGGGACTTAAGCGCAGTGCAGTCGGCCa                              <  1:4397416/43‑1 (MQ=37)                                                              aaTGGATTCCGCCAATGCGACTTAAGCGCAGACCAGTCgg                                 >  1:3485846/1‑40 (MQ=255)                                                              aaTGGATTCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATc            >  1:45550/1‑61 (MQ=255)                                                              aaTGGATGCCGCCAATGCGACTTAAGCGCAGACCAGt                                    >  1:208412/1‑37 (MQ=255)                                                              aaTGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATc  >  1:3103507/1‑71 (MQ=255)                                                                  gATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTc                           >  1:628511/1‑42 (MQ=255)                                                                    tGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCtt               <  1:3350884/52‑1 (MQ=255)                                                                     |                                                                TATCTCATAAATAATTCACGTTGTCGCCATAACGGCGACAACGTGAACGAAGATGGCTATTAATGGATGCCGCCAATGCGACTTAAGCGCAGACCAGTCGGCCATTCGCCTTCTTGCTTATCGAAGCTCATC  >  W3110S.gb/2702932‑2703063

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