Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,166 T→C 13.8% V357V (GTT→GTC yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,1660TC13.8% 336.3 / 42.0 130V357V (GTT→GTCyibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (47/65);  new base C (12/6);  total (59/71)
Fisher's exact test for biased strand distribution p-value = 7.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.00e-01

Only 100 of 140 total aligned reads displayed.

ATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGG  >  minE/2571096‑2571235
                                                                      |                                                                     
atgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGtt                                                                       <  1:46342/71‑1 (MQ=255)
   atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     <  1:3890893/70‑1 (MQ=255)
   atatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                                                                    <  1:2177103/71‑1 (MQ=255)
    tatCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     >  1:3011480/1‑69 (MQ=255)
      tCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                                      <  1:1774778/66‑1 (MQ=255)
        gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                                                <  1:1986810/70‑1 (MQ=255)
          tGCCTTACCCGACGGCATCTACGCACAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                                             >  1:3582736/1‑71 (MQ=255)
           gctttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGtctc                                                                     <  1:2136055/59‑1 (MQ=255)
            ccTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCg                                                            >  1:3995500/1‑70 (MQ=255)
             cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCgg                                                           >  1:2943204/1‑70 (MQ=255)
             cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          <  1:604/71‑1 (MQ=255)
             cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          <  1:2756688/71‑1 (MQ=255)
             cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          <  1:1553612/71‑1 (MQ=255)
              ttACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          >  1:491004/1‑70 (MQ=255)
              ttACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          >  1:2457986/1‑70 (MQ=255)
               tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   <  1:95550/60‑1 (MQ=255)
               tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   <  1:105244/60‑1 (MQ=255)
               tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   <  1:128929/60‑1 (MQ=255)
               tACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          >  1:3770779/1‑69 (MQ=255)
                aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                         <  1:3116470/69‑1 (MQ=255)
                     aCGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                                  >  1:230951/1‑70 (MQ=255)
                     aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   <  1:3999344/54‑1 (MQ=255)
                     aCGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          <  1:717041/63‑1 (MQ=255)
                      cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                                 <  1:1012803/71‑3 (MQ=255)
                           tttACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaaa                                                        >  1:248902/3‑59 (MQ=255)
                           tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     <  1:2887505/46‑1 (MQ=255)
                           tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                         <  1:1554937/58‑1 (MQ=255)
                           tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                          <  1:1714441/57‑1 (MQ=255)
                           tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                             >  1:1434260/1‑70 (MQ=255)
                            cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                           >  1:2402462/1‑71 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     >  1:3264540/1‑44 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     >  1:3959796/1‑44 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   >  1:1451582/1‑46 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   >  1:1949026/1‑46 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCgg                                                                   >  1:1605564/1‑46 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                  >  1:1932371/1‑47 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                                        >  1:119096/1‑57 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                          >  1:133312/1‑71 (MQ=255)
                             tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                          >  1:111777/1‑71 (MQ=255)
                              aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     >  1:1845410/1‑43 (MQ=255)
                              aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                  <  1:292240/46‑1 (MQ=255)
                               cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                        <  1:3668258/71‑1 (MQ=255)
                                   cAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTc                                                                     <  1:1743392/38‑1 (MQ=255)
                                     ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                  <  1:3194704/39‑1 (MQ=255)
                                     ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                              <  1:2728731/59‑1 (MQ=255)
                                      gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                                 >  1:210843/1‑71 (MQ=255)
                                      gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtg                                                     <  1:3857489/51‑1 (MQ=255)
                                        ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                               >  1:545315/1‑71 (MQ=255)
                                         gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                               <  1:952344/70‑1 (MQ=255)
                                          cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                             >  1:2814993/1‑71 (MQ=255)
                                           ggTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc                            <  1:3009714/71‑1 (MQ=255)
                                           ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                                                >  1:2255146/1‑35 (MQ=255)
                                           ggTTTACTCCGACCATTTCAGCCATGTTACGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                 <  1:3312450/66‑1 (MQ=255)
                                            gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCa                           <  1:3919056/71‑1 (MQ=255)
                                             tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                              <  1:1068004/51‑1 (MQ=255)
                                             tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                              <  1:4142337/51‑1 (MQ=255)
                                              ttACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaaa                                                        >  1:1607759/1‑40 (MQ=255)
                                              ttACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                                  >  1:480131/1‑45 (MQ=255)
                                              ttACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                                  >  1:1349466/1‑45 (MQ=255)
                                              ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                        >  1:3380717/1‑56 (MQ=255)
                                              ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                          >  1:3169414/1‑70 (MQ=255)
                                               tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGc                                   >  1:2947690/1‑60 (MQ=255)
                                               tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                         <  1:230781/70‑1 (MQ=255)
                                                aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                        <  1:253048/70‑1 (MQ=255)
                                                aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                        <  1:248070/70‑1 (MQ=255)
                                                    cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                             >  1:2545533/1‑45 (MQ=255)
                                                    cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                        >  1:2593160/1‑50 (MQ=255)
                                                    cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                  >  1:3403223/1‑56 (MQ=255)
                                                    cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCCGCTAAGAATGACCAGCTGGATGCATTATCtt                    >  1:2423530/1‑70 (MQ=255)
                                                        cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                           >  1:930626/1‑59 (MQ=255)
                                                        cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                           >  1:3276909/1‑59 (MQ=255)
                                                         aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                  <  1:1667984/67‑1 (MQ=255)
                                                          tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                             <  1:2863215/39‑1 (MQ=255)
                                                          tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg             <  1:679027/71‑1 (MQ=255)
                                                            tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGa            <  1:2342399/70‑1 (MQ=255)
                                                              aGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa         <  1:2335195/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:44487/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:630160/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:3093849/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:795084/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:423590/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:925931/71‑1 (MQ=255)
                                                                ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       <  1:3976246/71‑1 (MQ=255)
                                                                 cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                    >  1:4082641/1‑57 (MQ=255)
                                                                 cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg      <  1:3087011/71‑1 (MQ=255)
                                                                 cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg      <  1:775427/71‑1 (MQ=255)
                                                                 cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg      <  1:3683127/71‑1 (MQ=255)
                                                                  aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                 <  1:2191325/43‑1 (MQ=255)
                                                                  aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt       >  1:3613364/1‑69 (MQ=255)
                                                                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                 <  1:2223298/42‑1 (MQ=255)
                                                                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                 <  1:1544684/42‑1 (MQ=255)
                                                                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                 <  1:419792/42‑1 (MQ=255)
                                                                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa         <  1:824364/66‑1 (MQ=255)
                                                                   tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt        <  1:2066690/67‑1 (MQ=255)
                                                                   tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                   <  1:2338841/56‑1 (MQ=255)
                                                                    gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                   <  1:623439/39‑1 (MQ=255)
                                                                     tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                  >  1:485237/1‑39 (MQ=255)
                                                                     tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                               >  1:2596392/1‑42 (MQ=255)
                                                                     tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt               <  1:2604808/58‑1 (MQ=255)
                                                                     tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg  >  1:4022860/1‑71 (MQ=255)
                                                                      |                                                                     
ATGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGG  >  minE/2571096‑2571235

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: