Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,185,915 A→G 48.8% M36T (ATG→ACG)  potB ← polyamine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,185,9150AG48.8% ‑5.3 / 59.7 41M36T (ATG→ACG) potBpolyamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (7/14);  new base G (16/4);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 4.37e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

CAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGA  >  W3110S.gb/1185854‑1185984
                                                             |                                                                     
cAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                              <  1:2584153/71‑1 (MQ=255)
cAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                              <  1:3180569/71‑1 (MQ=255)
cAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                              <  1:3660795/71‑1 (MQ=255)
   tAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCGTTTTGACAAAAc                                                           >  1:2081059/1‑71 (MQ=255)
    aGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCGTTTTGACAAAACt                                                          >  1:1270946/1‑71 (MQ=255)
    aGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCGTTTTGACAAAACt                                                          >  1:1553967/1‑71 (MQ=255)
    aGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACt                                                          >  1:2243479/1‑71 (MQ=255)
       aCTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGc                                                       <  1:2105542/71‑1 (MQ=255)
                 agagCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCgg                                             >  1:1650340/1‑71 (MQ=255)
                        aTCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCaaaa                                      <  1:2512613/71‑1 (MQ=255)
                         tCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCaaa                                       >  1:2375983/1‑69 (MQ=255)
                                gACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGGGTCGTCGCGGGTCAAAAAGCTg                                 <  1:1593717/68‑1 (MQ=255)
                                gACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTg                                 <  1:401774/68‑1 (MQ=255)
                                  cgcgTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGc                             >  1:2104195/1‑70 (MQ=255)
                                        tAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGat                      >  1:34202/1‑71 (MQ=255)
                                        tAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGa                       >  1:2704446/1‑70 (MQ=255)
                                        tAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGat                      <  1:925892/71‑1 (MQ=255)
                                        tAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGat                      <  1:2649060/71‑1 (MQ=255)
                                         aGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGat                      <  1:1483042/70‑1 (MQ=255)
                                            tATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatcat                   >  1:2316298/1‑70 (MQ=255)
                                            tATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatcat                   >  1:3591728/1‑70 (MQ=255)
                                            tATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatcat                   >  1:3306076/1‑70 (MQ=255)
                                             aTCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatcat                   >  1:3164238/1‑69 (MQ=255)
                                              tCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCGGGTGCCAATGAtcatca                 <  1:183280/70‑1 (MQ=255)
                                                  gCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCATGTTg            >  1:2140386/1‑71 (MQ=255)
                                                     tAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCaaaaa                                     >  1:2378096/1‑43 (MQ=255)
                                                     tAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCaaaaa                                     >  1:841027/1‑43 (MQ=255)
                                                     tAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCaaaaa                                     >  1:3229746/1‑43 (MQ=255)
                                                     tAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCaaaaa                                     >  1:281681/1‑43 (MQ=255)
                                                     tAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGat                      <  1:2864198/58‑1 (MQ=255)
                                                     tAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGat                      <  1:2779319/58‑1 (MQ=255)
                                                       aaGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatca                    >  1:1152492/1‑58 (MQ=255)
                                                       aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGTCAATGATCATCAGGTTggg          >  1:1220071/1‑68 (MQ=255)
                                                       aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTggg          >  1:3357043/1‑68 (MQ=255)
                                                       aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTggg          >  1:3593840/1‑68 (MQ=255)
                                                          aCCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGa      <  1:3107109/69‑1 (MQ=255)
                                                           ccATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGa      >  1:2583777/1‑68 (MQ=255)
                                                            cATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGa  <  1:198862/71‑1 (MQ=255)
                                                            cATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGa  <  1:3430059/71‑1 (MQ=255)
                                                            cATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGa  <  1:347191/71‑1 (MQ=255)
                                                            cATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGa  <  1:72756/71‑1 (MQ=255)
                                                             |                                                                     
CAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTGGGCAGAAAGA  >  W3110S.gb/1185854‑1185984

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