Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,277,106	A→G	14.6%	T84A (ACT→GCT)	sohA →	predicted regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,277,106	0	A	G	14.6%	112.3 / 13.3	48	T84A (ACT→GCT)	sohA	predicted regulator
Reads supporting (aligned to +/- strand):  ref base A (17/24);  new base G (2/5);  total (19/29)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

GGATCATACTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATG  >  W3110S.gb/3277042‑3277172                                                                 |                                                                   ggaTCATACTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAGCTCGTc                                                              <  1:3560944/71‑1 (MQ=255)         cTATGAATGCATTTTTGCGTTTTCTGGATGCAGTTATCCAGAACAACCCGCAAAAAGCTCGTCcat                                                           >  1:326377/1‑66 (MQ=255)         cTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAGCTCGTCcattcaa                                                       >  1:1025300/1‑70 (MQ=255)         cTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaa                                                       >  1:2564155/1‑70 (MQ=255)         cTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaa                                                       >  1:2559439/1‑70 (MQ=255)         cTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaa                                                       <  1:3543231/70‑1 (MQ=255)         cTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaa                                                       >  1:3906641/1‑70 (MQ=255)               aTGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattca                                                        <  1:3078121/63‑1 (MQ=255)                tGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaacattcaa                                                >  1:2936382/1‑70 (MQ=255)                tGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCattcaacattcaac                                               <  1:1491705/71‑1 (MQ=255)                tGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCattcaacattcaac                                               <  1:1546260/71‑1 (MQ=255)                  cATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattc                                                         <  1:2555026/59‑1 (MQ=255)                  cATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCattcaacattcaac                                               >  1:3656426/1‑69 (MQ=255)                  cATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCattcaacattcaac                                               >  1:1313232/1‑69 (MQ=255)                  cATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCAttcaacattcaaca                                              <  1:1593845/70‑1 (MQ=255)                   aTTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaa                                                       >  1:1797643/1‑60 (MQ=255)                    tttttGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCcattcaacattc                                                  <  1:2285253/64‑1 (MQ=255)                      tttGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCca                                                            >  1:2605657/1‑52 (MQ=255)                      tttGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGga                                          >  1:1738025/1‑70 (MQ=255)                       ttGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGga                                          >  1:1015364/1‑69 (MQ=255)                          cGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAGCTCGTCCATTCAACATTCAACAAGgaaaga                                      <  1:460122/70‑1 (MQ=255)                          cGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGgaaaga                                      <  1:1382361/70‑1 (MQ=255)                           gTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGgaaagaaa                                    <  1:3385647/71‑1 (MQ=255)                               tCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACtt                                 <  1:3101237/70‑1 (MQ=255)                               tCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACtt                                 <  1:2404204/70‑1 (MQ=255)                               tCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACtt                                 <  1:4024250/70‑1 (MQ=255)                               tCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACtt                                 <  1:1431645/70‑1 (MQ=255)                               tCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACtt                                 <  1:95768/70‑1 (MQ=255)                                        aGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTgg                         <  1:2056142/69‑1 (MQ=255)                                           tatCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTCGCATg                     <  1:2532924/70‑1 (MQ=255)                                            atCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGa                   >  1:1828905/1‑71 (MQ=255)                                             tCCAGAACAACCCGCAAAAAACTCGTCcattcaaca                                                     >  1:973030/1‑36 (MQ=255)                                               cAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGgaaagaaa                                    >  1:301122/1‑51 (MQ=255)                                               cAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGgaaagaaa                                    >  1:2015491/1‑51 (MQ=255)                                               cAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGgaaagaaa                                    >  1:533367/1‑51 (MQ=255)                                                 gAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTc               <  1:1907267/70‑1 (MQ=255)                                                 gAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTc               <  1:3322372/70‑1 (MQ=255)                                                 caacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTc               <  1:1935027/69‑1 (MQ=255)                                                  aacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCaa             <  1:2888824/71‑1 (MQ=255)                                                  aacaacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCaa             <  1:3502786/71‑1 (MQ=255)                                                    caacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAAc            >  1:1023130/1‑70 (MQ=255)                                                    caacCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACa           >  1:1655277/1‑71 (MQ=255)                                                         cGCAAAAAGCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATtgat      <  1:2291607/71‑1 (MQ=255)                                                         cGCAAAAAGCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATtgat      <  1:1093376/71‑1 (MQ=255)                                                         cGCAAAAAGCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATtgat      <  1:93997/71‑1 (MQ=255)                                                         cGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATtgat      <  1:3925324/71‑1 (MQ=255)                                                           cAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTc                              <  1:3244723/45‑1 (MQ=255)                                                           cAAAAAACTCCTCCATTCAACATTCAACAAGGAAAGAAACTTGTc                              <  1:2836568/45‑1 (MQ=255)                                                             aaaaaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATtgatgatg  <  1:3448638/71‑1 (MQ=255)                                                               aaaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATtgatgat   >  1:974110/1‑68 (MQ=255)                                                                aaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATg                     <  1:887186/49‑1 (MQ=255)                                                                aaCTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATg                     <  1:964386/49‑1 (MQ=255)                                                                 |                                                                   GGATCATACTATGAATGCATTTTTGCGTTTTCTGGATGCAGATATCCAGAACAACCCGCAAAAAACTCGTCCATTCAACATTCAACAAGGAAAGAAACTTGTCGCTGGCATGGACGTCAACATTGATGATG  >  W3110S.gb/3277042‑3277172

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