Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,206 A→G 7.4% M371V (ATG→GTG)  yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,2060AG 7.4% 262.5 / 11.0 94M371V (ATG→GTG) yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/48);  new base G (4/3);  total (43/51)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

Only 100 of 101 total aligned reads displayed.

GGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAG  >  minE/2571136‑2571271
                                                                      |                                                                 
ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   >  1:545315/1‑71 (MQ=255)
 gCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   <  1:952344/70‑1 (MQ=255)
  cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                                                 >  1:2814993/1‑71 (MQ=255)
  cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                                                 <  1:239614/71‑1 (MQ=255)
   ggTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc                                                                <  1:3009714/71‑1 (MQ=255)
    gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCa                                                               <  1:3919056/71‑1 (MQ=255)
      ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                                              >  1:3169414/1‑70 (MQ=255)
      ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                                              >  1:2338357/1‑70 (MQ=255)
       tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                                             <  1:230781/70‑1 (MQ=255)
       tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                                              >  1:3680293/1‑69 (MQ=255)
       tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                                            >  1:1855529/1‑71 (MQ=255)
        aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                                            <  1:253048/70‑1 (MQ=255)
        aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                                            <  1:248070/70‑1 (MQ=255)
            cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCCGCTAAGAATGACCAGCTGGATGCATTATCtt                                                        >  1:2423530/1‑70 (MQ=255)
                cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                                                               >  1:3276909/1‑59 (MQ=255)
                cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                                                               >  1:930626/1‑59 (MQ=255)
                 aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                                                      <  1:1667984/67‑1 (MQ=255)
                  tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                                 <  1:679027/71‑1 (MQ=255)
                    tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGa                                                <  1:2342399/70‑1 (MQ=255)
                    tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               >  1:552544/1‑71 (MQ=255)
                    tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               >  1:746082/1‑71 (MQ=255)
                      aGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                             <  1:2335195/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                                                      <  1:457302/60‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:423590/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:3093849/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:925931/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:3976246/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:795084/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:44487/71‑1 (MQ=255)
                        ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:630160/71‑1 (MQ=255)
                         cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                                                        >  1:4082641/1‑57 (MQ=255)
                         cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                          <  1:775427/71‑1 (MQ=255)
                         cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                          <  1:3087011/71‑1 (MQ=255)
                         cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                          <  1:3683127/71‑1 (MQ=255)
                         cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                          <  1:176674/71‑1 (MQ=255)
                          aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           >  1:3613364/1‑69 (MQ=255)
                          aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           >  1:2265557/1‑69 (MQ=255)
                           tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                             <  1:824364/66‑1 (MQ=255)
                           tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                                            <  1:2066690/67‑1 (MQ=255)
                           tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                                       <  1:2911709/56‑1 (MQ=255)
                           tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                                       <  1:2338841/56‑1 (MQ=255)
                             tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   >  1:2596392/1‑42 (MQ=255)
                             tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   >  1:148076/1‑42 (MQ=255)
                             tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                                                   <  1:2604808/58‑1 (MQ=255)
                             tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg                                      >  1:4022860/1‑71 (MQ=255)
                             tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                                                     >  1:1880178/3‑56 (MQ=255)
                             tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                                                     >  1:2445816/3‑56 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   >  1:1281653/1‑40 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   >  1:899206/1‑40 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                                            >  1:3041325/1‑47 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                                            >  1:1650799/1‑47 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                                            >  1:3582595/1‑47 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           >  1:3999649/1‑64 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                                            >  1:3111456/1‑63 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                     <  1:3462416/70‑1 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                     <  1:1448316/70‑1 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                     <  1:1567541/70‑1 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                     <  1:594621/70‑1 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                     <  1:3833739/70‑1 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                                     <  1:3936350/70‑1 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg                                      >  1:1463137/1‑69 (MQ=255)
                               tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                                         >  1:1228556/1‑66 (MQ=255)
                                cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                                    <  1:1839896/70‑1 (MQ=255)
                                 ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                                   >  1:3888278/1‑38 (MQ=255)
                                   tGATGCGGAAAGTGCTGCTAAGAATGACCATCTGGATGCATTATCTTTATTTTTATCa                                             <  1:2245708/58‑1 (MQ=255)
                                    gATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagggg                               >  1:1510288/1‑71 (MQ=255)
                                      tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                             <  1:2250476/55‑1 (MQ=255)
                                      tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                                        <  1:409267/60‑1 (MQ=255)
                                       gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATc                                                          >  1:1592872/1‑41 (MQ=255)
                                       gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                                             >  1:3914078/1‑54 (MQ=255)
                                          gAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                                             <  1:3005181/35‑1 (MQ=255)
                                          gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt                         >  1:295657/1‑71 (MQ=255)
                                          gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt                         >  1:3659699/1‑71 (MQ=255)
                                           aaaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt                         >  1:2594483/1‑70 (MQ=255)
                                            aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                                           <  1:1692103/51‑1 (MQ=255)
                                               tgctgcTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                                              <  1:1206929/45‑1 (MQ=255)
                                               tgctgcTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATg                     >  1:263586/1‑70 (MQ=255)
                                                gctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                                    >  1:695984/1‑54 (MQ=255)
                                                 ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa                  >  1:1196068/1‑71 (MQ=255)
                                                 ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGaa                  >  1:3857652/1‑71 (MQ=255)
                                                  tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                          <  1:3705773/46‑1 (MQ=255)
                                                     tAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                                              <  1:785489/39‑1 (MQ=255)
                                                      aaGAATGACCAGCTGGATGCATTATCTTTATTTGGATc                                              >  1:3502510/1‑38 (MQ=255)
                                                      aaGAATGACCAGCT‑‑ATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAc                 >  1:2990063/1‑65 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               <  1:2691761/35‑1 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               <  1:1603974/35‑1 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               <  1:2799994/35‑1 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               <  1:2799712/35‑1 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                               <  1:2464111/35‑1 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCt           >  1:1467265/1‑71 (MQ=255)
                                                        gAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCt           >  1:3198189/1‑71 (MQ=255)
                                                         aaTGACCAGCTGGATGCATTATCTTTATTTGGATc                                              <  1:652628/35‑1 (MQ=255)
                                                         aaTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATg                     >  1:3270006/1‑60 (MQ=255)
                                                          aTGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACtt               >  1:3955087/1‑65 (MQ=255)
                                                            gACCAGCTGGATGCATTATCTTTATTTGGATCATTg                                          <  1:826268/36‑1 (MQ=255)
                                                             cccAGCTGGATGCATTTTCTTTATTTGGATCATTgaggggg                                    >  1:3639601/2‑41 (MQ=39)
                                                             aCCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTg      >  1:680785/1‑71 (MQ=255)
                                                               cAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGTGGGAGATTATGGAACTTTGCTCCTTGTa    <  1:3827956/71‑1 (MQ=255)
                                                                  cTGGATGCATTATCTTTATTTGGATCATTGagggggaggggg                              >  1:878995/1‑42 (MQ=255)
                                                                  cTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAg  >  1:4040246/1‑70 (MQ=255)
                                                                      |                                                                 
GGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAACTTTGCTCCTTGTAAG  >  minE/2571136‑2571271

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: