Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,735,602 T→C 11.0% E157E (GAA→GAG era ← membrane‑associated, 16S rRNA‑binding GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,735,6020TC11.0% 227.5 / 17.7 82E157E (GAA→GAGeramembrane‑associated, 16S rRNA‑binding GTPase
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base C (5/4);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

GGAAGTGATGAGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGGCAGCAG  >  minE/1735532‑1735670
                                                                      |                                                                    
ggAAGTGATGAGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGtt                                                                      >  1:2260312/1‑71 (MQ=255)
ggAAGTGATGAGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGtt                                                                      >  1:4040201/1‑71 (MQ=255)
 gAAGTGATGAGTCGCTTCAGGTAGATGCTTACGCCCGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     >  1:1436114/1‑71 (MQ=255)
    gTGATGAGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCgg                                                                  >  1:542882/1‑71 (MQ=255)
       atgaGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGc                                                                 <  1:286712/69‑1 (MQ=255)
       atgaGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTCTCGGc                                                                 <  1:4414510/69‑1 (MQ=255)
        tgaGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGc                                                                 <  1:1767592/68‑1 (MQ=255)
             cGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTCTCGGCAGAGATTg                                                         >  1:2309235/1‑71 (MQ=255)
             cGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTCTCGGCAGAGATTg                                                         >  1:1961905/1‑71 (MQ=255)
              gCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTgg                                                        <  1:3054077/71‑1 (MQ=255)
              gCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTgg                                                        <  1:1564464/71‑1 (MQ=255)
               cTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGc                                                       <  1:1535428/71‑1 (MQ=255)
               cTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGc                                                       <  1:3376615/71‑1 (MQ=255)
                  cAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCa                                                                >  1:2569818/1‑59 (MQ=255)
                  cAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCa                                                                >  1:3835602/1‑59 (MQ=255)
                  cAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCa                                                                >  1:4706895/1‑59 (MQ=255)
                  cAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCa                                                                >  1:4687029/1‑59 (MQ=255)
                      tAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGatat                                                <  1:4307365/71‑1 (MQ=255)
                        gATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCATAGATTGGCACGATATc                                               <  1:3383294/70‑1 (MQ=255)
                        gATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATc                                               <  1:3763256/70‑1 (MQ=255)
                         aTGCTTACGGACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:3146755/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:3352437/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:3813019/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:3109318/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:2664069/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:3060940/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:2061434/47‑1 (MQ=255)
                         aTGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGttt                                                                     <  1:785108/47‑1 (MQ=255)
                           gCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAg                                            >  1:1217208/1‑70 (MQ=255)
                           gCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAg                                            >  1:3627026/1‑70 (MQ=255)
                           gCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAg                                            >  1:3826865/1‑70 (MQ=255)
                           gCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAg                                            >  1:3451761/1‑70 (MQ=255)
                           gCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAg                                            >  1:2659620/1‑70 (MQ=255)
                             ttACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGa                                             <  1:1005742/67‑1 (MQ=255)
                             ttACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAgg                                           <  1:3558688/69‑1 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGtt                                      >  1:1408170/1‑71 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGtt                                      >  1:2573264/1‑71 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGt                                       >  1:3943847/1‑70 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGt                                       >  1:3761683/1‑70 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTCTCGGCAgaga                                                            >  1:4653903/1‑49 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTCTCGGCAgaga                                                            >  1:4269339/1‑49 (MQ=255)
                                cGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTCTCGGCAgaga                                                            >  1:3282150/1‑49 (MQ=255)
                                  cACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTc                                     >  1:2640955/1‑70 (MQ=255)
                                  cACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCa                                    <  1:3969772/71‑1 (MQ=255)
                                  cACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCa                                    <  1:2152006/71‑1 (MQ=255)
                                  cACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCa                                    <  1:588615/71‑1 (MQ=255)
                                  cACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCa                                    <  1:2468675/71‑1 (MQ=255)
                                   aCGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCa                                    <  1:374947/70‑1 (MQ=255)
                                    cGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATc                                  >  1:2228539/1‑71 (MQ=255)
                                     gATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 <  1:1248830/71‑1 (MQ=255)
                                     gATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 <  1:3814548/71‑1 (MQ=255)
                                     gATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 <  1:3327859/71‑1 (MQ=255)
                                     gATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 >  1:194349/1‑71 (MQ=255)
                                     gATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 <  1:2515718/71‑1 (MQ=255)
                                     gATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 <  1:2841964/71‑1 (MQ=255)
                                      aTTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCAt                                   >  1:3269167/1‑68 (MQ=255)
                                      aTTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCAt                                   >  1:3031455/1‑68 (MQ=255)
                                      aTTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTg                                >  1:1443294/1‑71 (MQ=255)
                                          ccGCAATAGTGTCAACATTCAGCCCGGTTTCGGCtgag                                                             <  1:3771919/38‑1 (MQ=255)
                                            gCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGAt                                                           >  1:3716053/1‑38 (MQ=255)
                                            gCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGc                          <  1:4698178/71‑1 (MQ=255)
                                            gCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGc                          <  1:2284834/71‑1 (MQ=255)
                                            gCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGc                          <  1:433706/71‑1 (MQ=255)
                                            gCAATAGTGTCAACATTCAGCCCGGTCTCGGCAGAGATTGGCACGATATCGAgg                                           <  1:3322933/54‑1 (MQ=255)
                                                tAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAgg                      <  1:803499/71‑1 (MQ=255)
                                                   tgtCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTg                           <  1:363288/63‑1 (MQ=255)
                                                   tgtCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTg                           <  1:4601650/63‑1 (MQ=255)
                                                   tgtCAACATTCAGCCCGGTCTCGGCAGAGATTGGCACGATATCGAGGAAGTTc                                     <  1:4545974/53‑1 (MQ=255)
                                                   tgtCAACATTCAGCCCGGTCTCGGCAGAGATTGGCACGATATCGAGGAAGTTc                                     <  1:2809755/53‑1 (MQ=255)
                                                     tCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGcc                         >  1:814491/1‑63 (MQ=255)
                                                       aaCATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTg                 >  1:751647/1‑69 (MQ=255)
                                                       aaCATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTg                 >  1:1519449/1‑69 (MQ=255)
                                                       aaCATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTg                 >  1:515470/1‑69 (MQ=255)
                                                       aaCATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTg                 >  1:2101028/1‑69 (MQ=255)
                                                       aaCATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGc                <  1:4241739/70‑1 (MQ=255)
                                                          aTTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGaa                                         >  1:1907088/1‑42 (MQ=255)
                                                          aTTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGaa                                         >  1:3098239/1‑42 (MQ=255)
                                                          aTTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGaa                                         >  1:2573847/1‑42 (MQ=255)
                                                           ttCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAgg                                           <  1:1793153/39‑1 (MQ=255)
                                                           ttCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAgg                                           <  1:2835376/39‑1 (MQ=255)
                                                           ttCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTg           >  1:3408008/1‑71 (MQ=255)
                                                             cAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGaa                                         >  1:4454215/1‑39 (MQ=255)
                                                              aGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGaa                    >  1:829221/1‑59 (MQ=255)
                                                              aGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGaa                    >  1:3710521/1‑59 (MQ=255)
                                                              aGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCgg        <  1:1045645/71‑1 (MQ=255)
                                                                 ccGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGgcag     <  1:4226239/71‑1 (MQ=255)
                                                                 ccGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGgcag     <  1:3473044/71‑1 (MQ=255)
                                                                 ccGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGgcag     <  1:1938032/71‑1 (MQ=255)
                                                                 ccGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGgcag     <  1:112932/71‑1 (MQ=255)
                                                                  cGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGc                <  1:255252/59‑1 (MQ=255)
                                                                   ggTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCt                                 <  1:3961440/41‑1 (MQ=255)
                                                                   ggTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGgcag     >  1:1604827/1‑69 (MQ=255)
                                                                    gTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGGcagcag  <  1:2789403/71‑1 (MQ=255)
                                                                     tttCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGgcag     >  1:499856/1‑67 (MQ=255)
                                                                      |                                                                    
GGAAGTGATGAGTCGCTTCAGGTAGATGCTTACGCACGATTGCCGCAATAGTGTCAACATTCAGCCCGGTTTCGGCAGAGATTGGCACGATATCGAGGAAGTTCATCTGGCTTGCCAGGAACTGCAGGTGCGGCAGCAG  >  minE/1735532‑1735670

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: