Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 824,680 T→C 35.9% D172D (GAT→GAC chaC → regulatory protein for cation transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE824,6800TC35.9% 37.9 / 37.1 64D172D (GAT→GACchaCregulatory protein for cation transport
Reads supporting (aligned to +/- strand):  ref base T (28/13);  new base C (4/19);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.86e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ATGGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCGCTGGGAACCAACGCAC  >  minE/824616‑824744
                                                                |                                                                
atgGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTc                                                              >  1:1138881/1‑69 (MQ=255)
atgGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTc                                                              >  1:2541956/1‑69 (MQ=255)
atgGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTc                                                              >  1:330061/1‑69 (MQ=255)
atgGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTc                                                              >  1:3183669/1‑69 (MQ=255)
atgGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTc                                                              >  1:1902199/1‑69 (MQ=255)
           gTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGt                                                    >  1:888436/1‑68 (MQ=255)
           gTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGt                                                    >  1:2061953/1‑68 (MQ=255)
           gTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGt                                                    >  1:2588018/1‑68 (MQ=255)
           gTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGt                                                    >  1:3072837/1‑68 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCAt                                                 >  1:1675614/1‑68 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCAt                                                 >  1:1537371/1‑68 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:3145946/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:975251/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:280221/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:3149769/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:3373784/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:874849/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:216239/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:1249547/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:1251357/1‑69 (MQ=255)
              aaCGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATc                                                >  1:1521826/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCgcgc                                            >  1:2459037/1‑65 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGcc                                           >  1:3237487/1‑66 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGtt                                        >  1:9491/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGtt                                        >  1:2346830/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGtt                                        >  1:3187329/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGtt                                        >  1:1367911/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGtt                                        >  1:2904044/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGtt                                        >  1:2873567/1‑69 (MQ=255)
                      tgtgTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGt                                         >  1:3378950/1‑68 (MQ=255)
                       gtgtTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTg                                       <  1:210613/69‑1 (MQ=255)
                       gtgtTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTg                                       <  1:2005435/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTTAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:257282/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:2466179/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:1001801/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:2121950/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:924288/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:2222011/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:769747/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:751745/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:670634/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:604066/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:43017/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:3394129/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:1253500/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:327010/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:2780881/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:1902027/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:2905748/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:2492591/69‑1 (MQ=255)
                                     cccGCGACATCCAGAATATGAATCTGACACTCGCGCTCAGGTCATAGCGCCGTTGATTGCGGCGGCGAg                         <  1:1275311/69‑1 (MQ=255)
                                           aCATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:252658/69‑1 (MQ=255)
                                           aCATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:703908/69‑1 (MQ=255)
                                           aCATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:867596/69‑1 (MQ=255)
                                           aCATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:2714164/69‑1 (MQ=255)
                                            cATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:1679424/68‑1 (MQ=255)
                                            cATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:902837/68‑1 (MQ=255)
                                            cATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTcc                   <  1:2597721/68‑1 (MQ=255)
                                            cATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCg                  <  1:1711553/69‑1 (MQ=255)
                                            cATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCg                  <  1:3001366/69‑1 (MQ=255)
                                                cAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCGCTg               >  1:2456573/1‑68 (MQ=255)
                                                cAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCGCTg               >  1:963883/1‑68 (MQ=255)
                                                     tatGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCg                  <  1:46929/60‑1 (MQ=255)
                                                            cTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCGCTGGGAACCAACGcac  <  1:1925011/69‑1 (MQ=255)
                                                                |                                                                
ATGGACGCACAGTAAACGCCATTGTGTTTATTATGGACCCGCGACATCCAGAATATGAATCTGATACTCGCGCTCAGGTCATCGCGCCGTTGATTGCGGCGGCGAGCGGTCCGCTGGGAACCAACGCAC  >  minE/824616‑824744

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: