Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,578,721	T→C	8.6%	T25T (ACT→ACC)	glnA →	glutamine synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,578,721	0	T	C	8.6%	218.8 / 9.6	81	T25T (ACT→ACC)	glnA	glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (42/32);  new base C (4/3);  total (46/35)
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

CGCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCGAAGAAGGCAAAA  >  W3110S.gb/3578652‑3578791                                                                      |                                                                       cGCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACt                                                                        <  1:2315087/70‑1 (MQ=255) cGCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTa                                                                       <  1:2609813/71‑1 (MQ=255)  gCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACCaa                                                                      <  1:1889521/71‑1 (MQ=255)  gCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACCaa                                                                      <  1:1505141/71‑1 (MQ=255)  gCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACCaa                                                                      <  1:2625699/71‑1 (MQ=255)      aaCACGTACTGACGATGCTGAAAGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAgg                                                                   <  1:1613610/70‑1 (MQ=255)         acGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTa                                                                 >  1:1828880/1‑69 (MQ=255)          cGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAg                                                              <  1:2445856/71‑1 (MQ=255)           gTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAg                                                              >  1:146833/1‑70 (MQ=255)           gTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAg                                                              >  1:2840542/1‑70 (MQ=255)            tACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTa                                                                       <  1:1036902/60‑1 (MQ=255)              cTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTaaa                                                                     <  1:1886343/60‑1 (MQ=255)              cTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGaa                                                            >  1:2528407/1‑69 (MQ=255)              cTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGaa                                                            >  1:1983268/1‑69 (MQ=255)                gACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGATTCACCGATACTAAAGGTAAAGAACAg                                                         >  1:2817800/1‑70 (MQ=255)                 aCGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAg                                                         <  1:2247450/69‑1 (MQ=255)                 aCGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAg                                                         <  1:3009393/69‑1 (MQ=255)                 aCGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCa                                                       >  1:2999352/1‑71 (MQ=255)                    aTGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAATGTAAAGAACAGCACgg                                                    <  1:1430712/71‑2 (MQ=255)                    aTGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGt                                                    >  1:156336/1‑71 (MQ=255)                     tGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGt                                                    >  1:423528/1‑70 (MQ=255)                      gCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCa                                                  >  1:1656673/1‑71 (MQ=255)                       cTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCAc                                                 <  1:3817024/71‑1 (MQ=255)                       cTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCAc                                                 <  1:2267632/71‑1 (MQ=255)                       cTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCAc                                                 <  1:972385/71‑1 (MQ=255)                       cTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCAc                                                 <  1:1946173/71‑1 (MQ=255)                        tGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCAc                                                 >  1:741435/1‑70 (MQ=255)                         gAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTa                                               >  1:751501/1‑71 (MQ=255)                          aaCGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGaa                                                            <  1:1911801/57‑1 (MQ=255)                          aaCGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGaa                                                            <  1:1696576/57‑1 (MQ=255)                          aaCGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACt                                                <  1:2637415/69‑1 (MQ=255)                          aaCGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTAt                                              <  1:1204267/71‑1 (MQ=255)                          aaCGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTAt                                              <  1:1640661/71‑1 (MQ=255)                          aaCGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACAAt                                              >  1:1502008/1‑71 (MQ=255)                            cGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACt                                                >  1:1204869/1‑67 (MQ=255)                            cGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATcc                                            >  1:979764/1‑71 (MQ=255)                            cGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATcc                                            >  1:1869489/1‑71 (MQ=255)                              aGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATccc                                           >  1:2632436/1‑70 (MQ=255)                               gCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTaaa                                                               >  1:760837/1‑49 (MQ=255)                               gCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTg                                         >  1:910669/1‑71 (MQ=255)                               gCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTg                                         <  1:744437/71‑1 (MQ=255)                                  cGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCtc                                      <  1:3042880/71‑1 (MQ=255)                                   gaagtgaagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACg                                                     >  1:2947822/1‑55 (MQ=255)                                   gaagtgaagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCtca                                     >  1:64542/1‑71 (MQ=255)                                    aagtgaagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCtcat                                    <  1:2517156/71‑1 (MQ=255)                                      gtgaagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCtcat                                    >  1:1404561/1‑69 (MQ=255)                                       tgaagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCtcatca                                  <  1:2955432/70‑1 (MQ=255)                                         aagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGt                               >  1:1437804/1‑71 (MQ=255)                                         aagtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGt                               <  1:2359713/71‑1 (MQ=255)                                          agtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGt                               >  1:1425826/1‑70 (MQ=255)                                          agtTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTg                              <  1:2672849/71‑1 (MQ=255)                                           gtTCGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTg                              >  1:2131364/1‑70 (MQ=255)                                              tgttgATTTGCGCTTCACCGATACTAAAGGTAAAGAAc                                                           >  1:1907220/1‑38 (MQ=255)                                              tgttgATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATg                          >  1:2477362/1‑71 (MQ=255)                                                 tgATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCt                        >  1:3664925/1‑70 (MQ=255)                                                 tgATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCt                        >  1:503939/1‑70 (MQ=255)                                                 tgATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCt                        >  1:2191493/1‑70 (MQ=255)                                                 tgATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCt                        >  1:2062569/1‑70 (MQ=255)                                                 tgATTTCCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCt                        >  1:499253/1‑70 (MQ=255)                                                  gATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGc                         <  1:1198474/68‑1 (MQ=255)                                                  gATTTGCGCTTCACCGATACCAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGccg                       >  1:2435718/1‑68 (MQ=255)                                                  gATTTGCGCTTCACCGATACCAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGccg                       >  1:1516327/1‑68 (MQ=255)                                                  gATTTGCGCTTCACCGATACCAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGccg                       >  1:2192603/1‑68 (MQ=255)                                                  gATTTGCGCTTCACCGATACCAAAGGTAAAGAACAGCACGTCACCATCCCTGCTCATCAGGTGAATGccg                       >  1:866853/1‑68 (MQ=255)                                                   aTTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGaa                     <  1:1507082/71‑1 (MQ=255)                                                    tttGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGa                             >  1:1111309/1‑62 (MQ=255)                                                    tttGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGa                             >  1:1991058/1‑62 (MQ=255)                                                     ttGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAAtt                   >  1:350617/1‑71 (MQ=255)                                                      tGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAAttc                  >  1:3296713/1‑71 (MQ=255)                                                      tGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAAttc                  >  1:1668031/1‑71 (MQ=255)                                                       gcgcTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCt                                          >  1:1554404/1‑46 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaaga          <  1:3228209/71‑1 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaaga          <  1:862843/71‑1 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaaga          <  1:174508/71‑1 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaaga          <  1:739925/71‑1 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaaga          <  1:560204/71‑1 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaaga          <  1:2797967/71‑1 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaag           >  1:1699475/1‑70 (MQ=255)                                                              aCCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCgaag           >  1:2789857/1‑70 (MQ=255)                                                                    aCTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCGAAGAAGGCa     >  1:2224654/1‑70 (MQ=255)                                                                     cTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCGAAGAAGGCaa    >  1:686804/1‑70 (MQ=255)                                                                     cTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCGAAGAAGGCaa    <  1:2508459/70‑1 (MQ=255)                                                                      tAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCGAAGAAGGCaaaa  >  1:2964067/1‑71 (MQ=255)                                                                      |                                                                       CGCTGAACACGTACTGACGATGCTGAACGAGCACGAAGTGAAGTTTGTTGATTTGCGCTTCACCGATACTAAAGGTAAAGAACAGCACGTCACTATCCCTGCTCATCAGGTGAATGCTGAATTCTTCGAAGAAGGCAAAA  >  W3110S.gb/3578652‑3578791

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