Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,151 C→T 11.8% D352D (GAC→GAT)  yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,1510CT11.8% 370.5 / 30.7 127D352D (GAC→GAT) yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (65/47);  new base T (10/5);  total (75/52)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 129 total aligned reads displayed.

GGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTT  >  minE/2571085‑2571217
                                                                  |                                                                  
ggACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                                                >  1:1804343/1‑71 (MQ=255)
 gACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                                               >  1:2962501/1‑71 (MQ=255)
   cccTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                                                >  1:2413395/1‑68 (MQ=255)
     cTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAg                                                            <  1:157861/70‑1 (MQ=255)
      tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:2109264/1‑70 (MQ=255)
      tAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:1472993/1‑70 (MQ=255)
       aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAttt                                                               >  1:3609083/1‑65 (MQ=255)
       aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGcc                                                          >  1:3531010/1‑70 (MQ=255)
       aaCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCa                                                         >  1:1344110/1‑71 (MQ=255)
        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGcc                                                          <  1:1697132/69‑1 (MQ=255)
        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:2648577/1‑71 (MQ=255)
        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:605397/1‑71 (MQ=255)
        aCGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:2941390/1‑71 (MQ=255)
         cGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGAc                                                                    <  1:1458731/58‑1 (MQ=255)
           atgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGtt                                                     <  1:46342/71‑1 (MQ=255)
            tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                    >  1:3773897/1‑71 (MQ=255)
            tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGtt                                                     <  1:1428399/70‑1 (MQ=255)
             gatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                   >  1:3033650/1‑71 (MQ=255)
              atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                   <  1:3890893/70‑1 (MQ=255)
              atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                                                  <  1:2177103/71‑1 (MQ=255)
               tatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                   >  1:3011480/1‑69 (MQ=255)
                 tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGcc                                                          <  1:2051528/60‑1 (MQ=255)
                 tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                    <  1:1774778/66‑1 (MQ=255)
                   gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              <  1:1986810/70‑1 (MQ=255)
                   gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                                             <  1:1574837/71‑1 (MQ=255)
                     tGCCTTACCCGACGGCATCTACGCACAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                           >  1:3582736/1‑71 (MQ=255)
                      gctttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGtctc                                                   <  1:2136055/59‑1 (MQ=255)
                       ccTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCg                                          >  1:3995500/1‑70 (MQ=255)
                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAt                                                                 <  1:202761/46‑1 (MQ=255)
                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCgg                                         >  1:2943204/1‑70 (MQ=255)
                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:604/71‑1 (MQ=255)
                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:126297/71‑1 (MQ=255)
                        cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:2756688/71‑1 (MQ=255)
                         ttACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        >  1:491004/1‑70 (MQ=255)
                         ttACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        >  1:2457986/1‑70 (MQ=255)
                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:2861070/1‑50 (MQ=255)
                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:3310591/1‑50 (MQ=255)
                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 <  1:105244/60‑1 (MQ=255)
                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 <  1:95550/60‑1 (MQ=255)
                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 <  1:128929/60‑1 (MQ=255)
                          tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        >  1:3770779/1‑69 (MQ=255)
                               gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGc                                                           >  1:194186/1‑45 (MQ=255)
                                aCGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                >  1:230951/1‑70 (MQ=255)
                                aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCAtt                                                                <  1:74586/39‑1 (MQ=255)
                                aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 <  1:3999344/54‑1 (MQ=255)
                                aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 <  1:2820720/54‑1 (MQ=255)
                                aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:717041/63‑1 (MQ=255)
                                 cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                               <  1:1012803/71‑3 (MQ=255)
                                  ggCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTc                                                              >  1:3264489/1‑39 (MQ=255)
                                      tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                       <  1:1554937/58‑1 (MQ=255)
                                      tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:416130/57‑1 (MQ=255)
                                      tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:1714441/57‑1 (MQ=255)
                                      tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                           >  1:1434260/1‑70 (MQ=255)
                                       cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                         >  1:2402462/1‑71 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:2197641/1‑39 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:4102977/1‑39 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCAt                                                        >  1:1464632/1‑39 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                   >  1:3959796/1‑44 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                   >  1:3264540/1‑44 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 >  1:1949026/1‑46 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                 >  1:1605564/1‑46 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              >  1:3163798/1‑49 (MQ=255)
                                        tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                      >  1:119096/1‑57 (MQ=255)
                                         aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                   >  1:1845410/1‑43 (MQ=255)
                                         aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                <  1:292240/46‑1 (MQ=255)
                                          cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                        <  1:4080024/69‑1 (MQ=255)
                                          cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                      <  1:3668258/71‑1 (MQ=255)
                                             ccAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCa                   <  1:80103/71‑1 (MQ=255)
                                                ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                <  1:3194704/39‑1 (MQ=255)
                                                ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                            <  1:2728731/59‑1 (MQ=255)
                                                 gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg               >  1:1095014/1‑71 (MQ=255)
                                                 gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg               >  1:210843/1‑71 (MQ=255)
                                                 gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                        <  1:2594655/46‑1 (MQ=255)
                                                 gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtg                                   <  1:3857489/51‑1 (MQ=255)
                                                 gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg               >  1:3945463/1‑71 (MQ=255)
                                                   ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCgg                                         >  1:938061/1‑43 (MQ=255)
                                                   ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa             >  1:545315/1‑71 (MQ=255)
                                                    gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                        >  1:1617535/1‑59 (MQ=255)
                                                    gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa             <  1:952344/70‑1 (MQ=255)
                                                     cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg           <  1:239614/71‑1 (MQ=255)
                                                     cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg           >  1:2814993/1‑71 (MQ=255)
                                                      ggTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc          <  1:3009714/71‑1 (MQ=255)
                                                      ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                              >  1:2255146/1‑35 (MQ=255)
                                                       gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCa         <  1:3919056/71‑1 (MQ=255)
                                                        tttACTCTGATCATTTCAGCCATGTCTCTGTGATGCgc                                         <  1:2470087/38‑2 (MQ=19)
                                                        tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                            <  1:1068004/51‑1 (MQ=255)
                                                         ttACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaaa                                      >  1:3343647/1‑40 (MQ=255)
                                                         ttACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                >  1:1349466/1‑45 (MQ=255)
                                                         ttACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                >  1:480131/1‑45 (MQ=255)
                                                         ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                      >  1:3380717/1‑56 (MQ=255)
                                                         ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt        >  1:3169414/1‑70 (MQ=255)
                                                         ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt        >  1:2338357/1‑70 (MQ=255)
                                                          tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt       <  1:230781/70‑1 (MQ=255)
                                                          tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa      >  1:1855529/1‑71 (MQ=255)
                                                           aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa      <  1:248070/70‑1 (MQ=255)
                                                           aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa      <  1:253048/70‑1 (MQ=255)
                                                               cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                           >  1:2545533/1‑45 (MQ=255)
                                                               cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                           >  1:2717872/1‑45 (MQ=255)
                                                               cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                      >  1:2593160/1‑50 (MQ=255)
                                                               cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                >  1:3403223/1‑56 (MQ=255)
                                                                  |                                                                  
GGACCCTAACGATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTT  >  minE/2571085‑2571217

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