Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,523,433 G→T 24.9% A78S (GCC→TCC)  flk → predicted flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,523,4330GT24.9% 91.5 / 12.6 36A78S (GCC→TCC) flkpredicted flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base G (17/10);  new base T (7/2);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CTGTGGGCGGGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCTTT  >  minE/1523364‑1523499
                                                                     |                                                                  
cTGTGGGCGGGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGc                                                                   >  1:101162/1‑71 (MQ=255)
           tATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCAt                                                         >  1:2197560/1‑70 (MQ=255)
             tCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATa                                                        <  1:878762/69‑1 (MQ=255)
               aGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCa                                                          <  1:1450519/65‑1 (MQ=255)
                   aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGgcg                                                >  1:792730/1‑71 (MQ=255)
                   aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATATGgcg                                                >  1:1894131/1‑71 (MQ=255)
                     gATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCa                                              <  1:1077721/71‑1 (MQ=255)
                       tatTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAAc                                            <  1:781086/71‑1 (MQ=255)
                          tGGTCTGTCGGGAGATTCACCGGTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTc                                         >  1:1331746/1‑71 (MQ=255)
                           ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:470285/1‑71 (MQ=255)
                            gTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTctg                                       <  1:855600/71‑1 (MQ=255)
                               tGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTctgctg                                    >  1:1281654/1‑71 (MQ=255)
                               tGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTctgctg                                    >  1:1141098/1‑71 (MQ=255)
                               tGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAAt                                                      >  1:1316640/1‑53 (MQ=255)
                                   gggAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCcg                                >  1:987513/1‑71 (MQ=255)
                                       gaTTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATa                                                        <  1:1525236/43‑1 (MQ=255)
                                          tCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaa                          >  1:1684358/1‑70 (MQ=255)
                                            aCCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaaaa                        <  1:1485817/70‑1 (MQ=255)
                                              cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaaaaa                       >  1:1401554/1‑69 (MQ=255)
                                              cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaaaaa                       >  1:279447/1‑69 (MQ=255)
                                              cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc                     >  1:1156866/1‑71 (MQ=255)
                                               gctgctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAg                      >  1:997944/1‑69 (MQ=255)
                                               gctgctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc                     >  1:1991758/1‑70 (MQ=255)
                                                   ctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGAAACGTCTGCTGGCCGCGCAAAAAAGCCAt                  <  1:2084260/69‑1 (MQ=255)
                                                     ttCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTg              >  1:1202674/1‑71 (MQ=255)
                                                          gTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGc         >  1:1474428/1‑71 (MQ=255)
                                                          gTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGc         >  1:969321/1‑71 (MQ=255)
                                                           tCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGc         >  1:1472562/1‑70 (MQ=255)
                                                           tCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGcc        <  1:552271/71‑1 (MQ=255)
                                                             aCTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAg      >  1:553459/1‑71 (MQ=255)
                                                             aCTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAg      >  1:1971287/1‑71 (MQ=255)
                                                             aCTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAg      >  1:1265149/1‑71 (MQ=255)
                                                               ttCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCaaa                          <  1:850902/49‑1 (MQ=255)
                                                                 cccTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCttt  >  1:7928/1‑71 (MQ=255)
                                                                   cTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc                     <  1:538945/50‑1 (MQ=255)
                                                                   cTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc                     <  1:1218433/50‑1 (MQ=255)
                                                                     |                                                                  
CTGTGGGCGGGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCTTT  >  minE/1523364‑1523499

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: