Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	4,509,018	A→G	6.5%	G93G (GGA→GGG)	yjhB →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	4,509,018	0	A	G	6.5%	258.3 / 8.5	92	G93G (GGA→GGG)	yjhB	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (32/54);  new base G (2/4);  total (34/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGG  >  W3110S.gb/4508952‑4509086                                                                   |                                                                     tGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAg                                                                  <  1:619632/71‑1 (MQ=255)        gCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTa                                                            >  1:3605413/1‑70 (MQ=255)        gCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAg                                                           >  1:1068510/1‑71 (MQ=255)        gCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAg                                                           >  1:2899458/1‑71 (MQ=255)         cTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGc                                                          >  1:3164277/1‑71 (MQ=255)         cTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGc                                                          >  1:514385/1‑71 (MQ=255)          tGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCg                                                         >  1:3747895/1‑71 (MQ=255)            gtAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:1224077/69‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACa                                                                   <  1:329379/59‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:2846395/70‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:1165437/70‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:2781091/70‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:3189275/70‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:1517037/70‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:1565330/70‑1 (MQ=255)            aTAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:1648661/70‑1 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGGACAGGCCTTAGCGGTAt                                                     >  1:268504/1‑71 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGc                                                                >  1:2832930/1‑60 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGc                                                                >  1:1696579/1‑60 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        >  1:3860300/1‑68 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTAt                                                     >  1:3554773/1‑71 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTAt                                                     >  1:2959943/1‑71 (MQ=255)              aaaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTAt                                                     >  1:630440/1‑71 (MQ=255)               aaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTAt                                                     <  1:1554055/70‑1 (MQ=255)               aaTATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTAt                                                     <  1:1495610/70‑1 (MQ=255)               aaTATGGTCGTAAGCCAATGATGATATGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTAtt                                                    >  1:2925734/1‑71 (MQ=255)                   tGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTa                                                >  1:3665627/1‑71 (MQ=255)                        gTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAAc                                           >  1:1492579/1‑71 (MQ=255)                         tAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACt                                          <  1:2352883/71‑1 (MQ=255)                             ccAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGGACAGGCCTTAGCGGTATTGCTACAAACTtata                                      <  1:2071941/71‑1 (MQ=255)                              cAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:110426/52‑1 (MQ=255)                              cAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:4572/52‑1 (MQ=255)                              cAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTtata                                      >  1:1259036/1‑70 (MQ=255)                              cAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTtata                                      >  1:2070617/1‑70 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAg                                                                  <  1:323402/40‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAg                                                                  <  1:2132776/40‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAg                                                                  <  1:158718/40‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAg                                                                  <  1:3207124/40‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:574446/50‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:746088/50‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACt                                          <  1:1782039/64‑1 (MQ=255)                                atgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACt                                          <  1:1790998/64‑1 (MQ=255)                                 tgatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATg                                  >  1:292546/1‑71 (MQ=255)                                  gatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGc                                                                >  1:991121/1‑40 (MQ=255)                                  gatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGc                                                                >  1:2163222/1‑40 (MQ=255)                                  gatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGc                                 >  1:1999658/1‑71 (MQ=255)                                  gatgatgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGc                                 >  1:3647636/1‑71 (MQ=255)                                      atgTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCg                              <  1:850697/70‑1 (MQ=255)                                        gtgGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAg                           >  1:88886/1‑71 (MQ=255)                                        gtgGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAg                           >  1:1666459/1‑71 (MQ=255)                                         tgGGCAATTTTCATTTACTCAGTCGGAACAGGCCtt                                                             <  1:1461413/36‑1 (MQ=255)                                         tgGGCAATTTTCATTTACTCAGTCGGAACAGGCCtt                                                             <  1:2064482/36‑1 (MQ=255)                                         tgGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCg                              <  1:2875846/67‑1 (MQ=255)                                         tgGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCg                              <  1:2608764/67‑1 (MQ=255)                                               aTTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:2914800/35‑1 (MQ=255)                                               aTTTTCATTTACTCAGTCGGAACAGGCCTTAGCgg                                                        <  1:859832/35‑1 (MQ=255)                                                ttttCATTTACTCAGTCGGGACAGGCCTTAGCGGTAt                                                     <  1:1792528/37‑1 (MQ=255)                                                ttttCATTTACTCAGTCGGGACAGGCCTTAGCGGTAt                                                     <  1:2145718/37‑1 (MQ=255)                                                ttttCATTTACTCAGTCGGGACAGGCCTTAGCGGTAt                                                     <  1:2003382/37‑1 (MQ=255)                                                ttttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACaa                                             <  1:3397551/45‑1 (MQ=255)                                                ttttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACaa                                             <  1:1953288/45‑1 (MQ=255)                                                ttttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTtata                                      <  1:883130/52‑1 (MQ=255)                                                 tttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGc                                 >  1:411933/1‑56 (MQ=255)                                                  ttCATTTACTCAGTCGGGACAGGCCTTAGCGGTATTg                                                   >  1:1001272/1‑37 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTta                                        <  1:3239788/48‑1 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGttt                 >  1:1975897/1‑71 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGttt                 >  1:1741256/1‑71 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACT‑‑TATATGCTCGCAGTTTGCCGtttt                <  1:2590603/70‑1 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACT‑‑TATATGCTCGCAGTTTGCCGtttt                <  1:2693395/70‑1 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACT‑‑TATATGCTCGCAGTTTGCCGtttt                <  1:776914/70‑1 (MQ=255)                                                  ttCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACT‑‑TATATGCTCGCAGTTTGCCGtttt                <  1:1363584/70‑1 (MQ=255)                                                   tCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGttt                 <  1:2588826/70‑1 (MQ=255)                                                    cATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCa                            <  1:1901371/58‑1 (MQ=255)                                                    cATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCa                            <  1:2149108/58‑1 (MQ=255)                                                      tttACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGt                          >  1:3147845/1‑58 (MQ=255)                                                      tttACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAtt             <  1:1501705/71‑1 (MQ=255)                                                        tACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAttgt           >  1:3667192/1‑71 (MQ=255)                                                          cTCAGTCGGAACAGGCCTTAGGGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAttgtt          <  1:1857249/70‑1 (MQ=255)                                                          cTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAttgtt          <  1:373112/70‑1 (MQ=255)                                                           tCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAttgttg         <  1:3277927/70‑1 (MQ=255)                                                            cAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTtat                                       <  1:1584965/39‑1 (MQ=255)                                                              gTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAttgt           <  1:1038374/65‑1 (MQ=255)                                                              gTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCtt     <  1:1290999/71‑1 (MQ=255)                                                               tCGGAACAGGCCTTAGCGGTATTGCTACAAACTtata                                      <  1:3351009/37‑1 (MQ=255)                                                               tCGGAACAGGCCTTAGCGGTATTGCTACAAACTtata                                      <  1:2588953/37‑1 (MQ=255)                                                               tCGGAACAGGCCTTAGCGGTATTGCTACAAACTta                                        >  1:1310684/1‑35 (MQ=255)                                                                cGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTAttgt           >  1:3934706/1‑63 (MQ=255)                                                                 ggAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAgg  <  1:722850/71‑1 (MQ=255)                                                                 ggAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAgg  <  1:2156122/71‑1 (MQ=255)                                                                 ggAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAgg  >  1:2293165/1‑71 (MQ=255)                                                                 ggAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAgg  <  1:3024650/71‑1 (MQ=255)                                                                  gAACAGGCCTTAGCGGTTTTGCTACAAACT‑‑TATATGCTCGCAGTTTGCCGtttt                <  1:535154/54‑1 (MQ=255)                                                                   |                                                                     TGCCATGGCTGATAAATATGGTCGTAAGCCAATGATGATGTGGGCAATTTTCATTTACTCAGTCGGAACAGGCCTTAGCGGTATTGCTACAAACTTATATATGCTCGCAGTTTGCCGTTTTATTGTTGGCTTAGG  >  W3110S.gb/4508952‑4509086

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