Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,185,915 A→G 39.4% 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,9150AG39.4% 18.2 / 35.9 33M36T (ATG→ACG) potBpolyamine transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (14/6);  new base G (9/4);  total (23/10)
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.61e-01

GTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTG  >  W3110S.gb/1185851‑1185974
                                                                |                                                           
gTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCAtttt                                                         <  1:1601393/69‑1 (MQ=255)
      tAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCGTTTTGACaaa                                                   <  1:2545685/69‑1 (MQ=255)
             tCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGaaa                                                     >  1:1771033/1‑58 (MQ=255)
             tCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                    >  1:104116/1‑61 (MQ=255)
             tCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                    >  1:3230691/1‑61 (MQ=255)
             tCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                    >  1:1896416/1‑61 (MQ=255)
             tCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACaa                                                    >  1:1557657/1‑61 (MQ=255)
              cAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGcgt                                           >  1:2430549/1‑69 (MQ=255)
                     gagCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTcgcg                                    <  1:2807994/69‑1 (MQ=255)
                                 cAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCt                        <  1:3524721/69‑1 (MQ=255)
                                 cAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCt                        <  1:260232/69‑1 (MQ=255)
                                        gtgtAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCa                 <  1:683983/69‑1 (MQ=255)
                                        gtgtAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCa                 <  1:1314266/69‑1 (MQ=255)
                                        gtgtAGTTATCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCa                 <  1:3233352/69‑1 (MQ=255)
                                        gtgtAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCa                 <  1:43534/69‑1 (MQ=255)
                                        gtgtAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCa                 <  1:3205248/69‑1 (MQ=255)
                                                aTCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatc           >  1:1198348/1‑67 (MQ=255)
                                                aTCCAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatc           >  1:1142856/1‑67 (MQ=255)
                                                aTCCAGCGTAAAGACCGTTTTCACAAAACTGGCGTCGTCGCGGGTCaaa                             >  1:2267324/1‑49 (MQ=255)
                                                aTCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatc           >  1:2429140/1‑67 (MQ=255)
                                                aTCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGatc           >  1:3018243/1‑67 (MQ=255)
                                                  ccAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCt                        >  1:724545/1‑52 (MQ=255)
                                                   cAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAg      >  1:970462/1‑69 (MQ=255)
                                                   cAGCGTAAAGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAg      >  1:1820112/1‑69 (MQ=255)
                                                        tAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGCTCATCAGGtt   >  1:2011576/1‑67 (MQ=255)
                                                          aaGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:2939345/1‑66 (MQ=255)
                                                          aaGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:2476125/1‑66 (MQ=255)
                                                          aaGACCGTTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:1246557/1‑66 (MQ=255)
                                                          aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:2588682/1‑66 (MQ=255)
                                                          aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:1980651/1‑66 (MQ=255)
                                                          aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:1859412/1‑66 (MQ=255)
                                                          aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:156438/1‑66 (MQ=255)
                                                          aaGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTg  >  1:1054892/1‑66 (MQ=255)
                                                                |                                                           
GTGCAATAGCACTTCAAAATAGAGCGGATCGAGCAGACGCGTGTAGTTATCCAGCGTAAAGACCATTTTGACAAAACTGGCGTCGTCGCGGGTCAAAAAGCTGGTGCCAATGATCATCAGGTTG  >  W3110S.gb/1185851‑1185974

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