Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,878,535 T→C 35.3% G210G (GGT→GGC cycA → D‑alanine/D‑serine/glycine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,878,5350TC35.3% 39.8 / 36.8 65G210G (GGT→GGCcycAD‑alanine/D‑serine/glycine transporter
Reads supporting (aligned to +/- strand):  ref base T (8/34);  new base C (14/9);  total (22/43)
Fisher's exact test for biased strand distribution p-value = 1.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTAC  >  minE/2878470‑2878603
                                                                 |                                                                    
tCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:2394616/69‑1 (MQ=255)
tCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:1684022/69‑1 (MQ=255)
tCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:177210/69‑1 (MQ=255)
tCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:2897973/69‑1 (MQ=255)
tCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:2851234/69‑1 (MQ=255)
tCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:2540032/69‑1 (MQ=255)
 cACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTaa                                                                  >  1:463178/1‑69 (MQ=255)
 cACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTaa                                                                  >  1:1943691/1‑69 (MQ=255)
 cACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTaa                                                                  >  1:1268271/1‑69 (MQ=255)
 cACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTaa                                                                  >  1:1867537/1‑69 (MQ=255)
    cGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGt                                                                <  1:1660343/68‑1 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGCCGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:106567/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1900426/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1869252/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1834548/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:222567/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:3023958/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1083427/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1017319/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1262167/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:447825/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:3252109/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:3122264/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgcg                                                               >  1:1526883/1‑67 (MQ=255)
    cGACCGGTGTGGAAGCATCATTTGCACATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGCTTAAgc                                                                >  1:2180391/1‑67 (MQ=255)
     gACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGt                                                                <  1:2525710/67‑1 (MQ=255)
               gAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGc                                                     >  1:3390497/1‑68 (MQ=255)
               gAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGc                                                     >  1:2242116/1‑68 (MQ=255)
               gAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGc                                                     >  1:2508767/1‑68 (MQ=255)
                       aTTCGCGCATTTGTTGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:45246/46‑1 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGtta                                                                    <  1:2350335/45‑2 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:2862187/46‑1 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:104689/46‑1 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:842079/46‑1 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:968055/46‑1 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:1072870/46‑1 (MQ=255)
                       aTTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTa                                                                   <  1:1504640/46‑1 (MQ=255)
                          cgcgCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGAt                                         >  1:2592278/1‑69 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:2241476/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:999063/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:1221632/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:1290777/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:3366104/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:3130887/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:23491/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:1899721/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:2022135/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:2175665/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:2668213/69‑1 (MQ=255)
                                      gAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGGTTTCGc                             <  1:99118/69‑1 (MQ=255)
                                       aaTGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGc                             <  1:31310/68‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:430770/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:987771/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:1048044/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:2534636/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:1698337/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:3434927/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:3279041/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:15018/56‑1 (MQ=255)
                                         tGACGGCGGCTGGTTCCCGAAAGGCTTAAGCGGCTTCTTTGCCGGATTCCAGATAg                                       <  1:1974244/56‑1 (MQ=255)
                                               cggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTggg                    <  1:1254460/69‑1 (MQ=255)
                                               cggcTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTggg                    <  1:757826/69‑1 (MQ=255)
                                                    ggTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTg               <  1:64707/69‑1 (MQ=255)
                                                    ggTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTg               <  1:3230824/69‑1 (MQ=255)
                                                    ggTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTg               <  1:3027327/69‑1 (MQ=255)
                                                    ggTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTACAGATAGCGGTTTTCGCTTTCGTGGGGATTg               <  1:1310653/69‑1 (MQ=255)
                                                                 tttAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTac  <  1:787347/69‑1 (MQ=255)
                                                                 |                                                                    
TCACCGACTGGTGTGGAAGCGTCATTCGCGCATTTGTGGAATGACGGCGGCTGGTTCCCGAAAGGTTTAAGTGGCTTCTTTGCCGGATTCCAGATAGCGGTTTTCGCTTTCGTGGGGATTGAGCTGGTAGGTAC  >  minE/2878470‑2878603

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: