Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,523,433 G→T 29.2% 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,4330GT29.2% 85.8 / 23.2 41A78S (GCC→TCC) flkpredicted flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base G (20/9);  new base T (9/3);  total (29/12)
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.22e-01

GGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCTTT  >  minE/1523373‑1523499
                                                            |                                                                  
ggTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCa                                                          <  1:4072489/71‑1 (MQ=255)
  tATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCTTCACTTCCCTGCCGCTGAGCATa                                                        <  1:3734627/71‑1 (MQ=255)
    tCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAAt                                                      <  1:1121344/71‑1 (MQ=255)
         cACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGgc                                                 >  1:800784/1‑71 (MQ=255)
          aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGgcg                                                >  1:606130/1‑71 (MQ=255)
          aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGgcg                                                >  1:27827/1‑71 (MQ=255)
          aCGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGgcg                                                >  1:2377015/1‑71 (MQ=255)
               atTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACg                                           >  1:311033/1‑71 (MQ=255)
               atTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACg                                           >  1:765224/1‑71 (MQ=255)
               atTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACg                                           >  1:3197688/1‑71 (MQ=255)
               atTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACg                                           >  1:1065587/1‑71 (MQ=255)
                  ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:1982076/1‑71 (MQ=255)
                  ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:1546216/1‑71 (MQ=255)
                  ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:3292627/1‑71 (MQ=255)
                  ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:1438711/1‑71 (MQ=255)
                  ggTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:3840319/1‑71 (MQ=255)
                   gTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct                                        <  1:2254851/70‑1 (MQ=255)
                   gTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct                                        <  1:3777420/70‑1 (MQ=255)
                        tCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTctgct                                     <  1:2845640/68‑1 (MQ=255)
                        tCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTctgct                                     <  1:3321152/68‑1 (MQ=255)
                           ggAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCa                                                          <  1:2820415/44‑1 (MQ=255)
                            gagaTTCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCa                                              >  1:470470/1‑55 (MQ=255)
                                ttCACCGCTGCTTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACg                                           <  1:1615768/54‑1 (MQ=255)
                                     cGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGc                     <  1:4076140/71‑1 (MQ=255)
                                         gctTTCGCGTCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCAtt                 >  1:2779977/1‑71 (MQ=255)
                                            ttCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTctg                                       >  1:3928848/1‑46 (MQ=255)
                                            ttCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTct                                        >  1:2108983/1‑45 (MQ=255)
                                             tCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCaa                                             >  1:1961125/1‑39 (MQ=255)
                                                  tCACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCAt                  <  1:2974848/61‑1 (MQ=255)
                                                   cACTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCAt                  <  1:1438788/60‑1 (MQ=255)
                                                    aCTTCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTATGCCCGCCAg      >  1:553499/1‑71 (MQ=255)
                                                      ttCCCTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCt    >  1:755643/1‑71 (MQ=255)
                                                      ttCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTc                >  1:1267845/1‑59 (MQ=255)
                                                       tCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTgg                                   >  1:3574750/1‑39 (MQ=255)
                                                       tCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTgg                                   >  1:4018458/1‑39 (MQ=255)
                                                       tCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTgg                                   >  1:3225826/1‑39 (MQ=255)
                                                       tCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTgg                                   >  1:1214576/1‑39 (MQ=255)
                                                       tCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCTTTc                >  1:546662/1‑58 (MQ=255)
                                                       tCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTc                >  1:2881717/1‑58 (MQ=255)
                                                         ccTTCCGGTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCttt  >  1:815182/1‑70 (MQ=255)
                                                         ccTTCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCttt  >  1:163029/1‑70 (MQ=255)
                                                            |                                                                  
GGTATCAGGCACGATATTGGTCTGTCGGGAGATTCACCGCTGCTTTCGCGTCACTTCCCTGCCGCTGAGCATAATCTGGCGCAACGTCTGCTGGCCGCGCAAAAAAGCCATTCTGCCCGCCAGCTTT  >  minE/1523373‑1523499

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