Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,045,818 C→T 10.4% intergenic (+580/+215) yqiK → / ← rfaE conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,045,8180CT10.4% 143.2 / 5.7 48intergenic (+580/+215)yqiK/rfaEconserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base C (15/28);  new base T (2/3);  total (17/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01

AGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGT  >  minE/2045752‑2045887
                                                                  |                                                                      
aGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTcc                                                                    >  1:637868/1‑71 (MQ=255)
  cATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTccc                                                                   <  1:2102473/70‑1 (MQ=255)
   aTTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCga                                                                 <  1:3295006/71‑1 (MQ=255)
         cgGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGaagaag                                                            <  1:3649451/70‑1 (MQ=255)
          ggggTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAg‑gg                                                         <  1:3065948/71‑1 (MQ=255)
           gggTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAg‑gg                                                         <  1:904446/70‑1 (MQ=255)
              gccGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCt                                                     <  1:753594/70‑1 (MQ=255)
                   cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                <  1:2460908/71‑1 (MQ=255)
                   cgcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGta                                                 >  1:934496/1‑70 (MQ=255)
                    gcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtata                                               <  1:1240967/71‑1 (MQ=255)
                    gcTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtata                                               <  1:893500/71‑1 (MQ=255)
                      tcatcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                <  1:3810546/68‑1 (MQ=255)
                        atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGta                                                 >  1:299217/1‑65 (MQ=255)
                        atcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAgg                                            >  1:1848227/1‑70 (MQ=255)
                         tcaACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGc                                                       <  1:3191778/58‑1 (MQ=255)
                          caacaaTCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGa                                           >  1:383422/1‑69 (MQ=255)
                                   gggggCAGCAAGGGGCGGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCtt                                                    <  1:2321566/51‑1 (MQ=255)
                                    ggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAg‑ggg                                                        <  1:10232/46‑1 (MQ=255)
                                    ggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGtat                                                <  1:379648/54‑1 (MQ=255)
                                    ggggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatgat                               <  1:2822319/71‑1 (MQ=255)
                                      ggCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTAtgatga                             <  1:2250737/71‑1 (MQ=255)
                                       gcagcaAGGGGCTGAAACGGGAAAGCCCCTCCCgaaga                                                              <  1:424134/38‑1 (MQ=255)
                                          gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCt                         >  1:3826505/1‑71 (MQ=255)
                                           caAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAg                           >  1:1959571/1‑68 (MQ=255)
                                            aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaaa                                         >  1:2838940/1‑53 (MQ=255)
                                            aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                       >  1:1376259/1‑71 (MQ=255)
                                               acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGt                      >  1:1708430/3‑70 (MQ=0)
                                                ggCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                   <  1:1936945/71‑1 (MQ=255)
                                                 gCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTatgatg                              <  1:3100332/59‑1 (MQ=255)
                                                   tGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                       <  1:748236/64‑1 (MQ=255)
                                                    gAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAgg                                       >  1:3395419/1‑47 (MQ=255)
                                                     aaaCGTGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                <  1:3550215/69‑1 (MQ=255)
                                                     aaaCGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTAtgatga                             >  1:116466/1‑57 (MQ=0)
                                                     aaaCGGGAAAGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcat                <  1:3128594/69‑1 (MQ=11)
                                                     aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGaaa                                         <  1:2098056/44‑1 (MQ=255)
                                                     aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTat                                  <  1:3155389/51‑1 (MQ=255)
                                                     aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                <  1:2923315/69‑1 (MQ=255)
                                                     aaaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc              >  1:673150/1‑71 (MQ=255)
                                                      aaCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc              <  1:2633942/70‑1 (MQ=255)
                                                       aCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTTTATAAg                                             >  1:3420149/1‑38 (MQ=255)
                                                       aCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                        >  1:3282549/1‑59 (MQ=255)
                                                       aCGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                        >  1:2261239/1‑59 (MQ=255)
                                                         gggAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTgg           <  1:1049747/71‑1 (MQ=11)
                                                         gggAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTgg           <  1:3376571/71‑1 (MQ=11)
                                                          ggAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTa                                   <  1:3541859/46‑1 (MQ=0)
                                                           gAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg        >  1:2125140/1‑71 (MQ=255)
                                                            aaaGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgt       <  1:484967/71‑1 (MQ=255)
                                                              aGCCCCTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatc               <  1:1334196/61‑2 (MQ=9)
                                                                 cccTCCCGAGGAAG‑GGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatcg              <  1:994508/59‑3 (MQ=255)
                                                                 cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAg                           <  1:273649/46‑1 (MQ=255)
                                                                 cccTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt  <  1:3639705/71‑1 (MQ=255)
                                                                  |                                                                      
AGCATTGCGCGGGGTCCGGCGCTCATCAACAATCGGGGGGCAGCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAG‑GGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGT  >  minE/2045752‑2045887

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: ‑