Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,178 G→C 14.5% R361R (CGG→CGC yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,1780GC14.5% 331.3 / 41.8 131R361R (CGG→CGCyibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (49/63);  new base C (14/5);  total (63/68)
Fisher's exact test for biased strand distribution p-value = 2.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.80e-01

Only 100 of 137 total aligned reads displayed.

CTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAT  >  minE/2571109‑2571248
                                                                     |                                                                      
cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:126297/71‑1 (MQ=255)
cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:604/71‑1 (MQ=255)
cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:1553612/71‑1 (MQ=255)
 ttACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       >  1:2457986/1‑70 (MQ=255)
 ttACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       >  1:491004/1‑70 (MQ=255)
   aCCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                                      <  1:3116470/69‑1 (MQ=255)
        aCGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgcta                                                               >  1:230951/1‑70 (MQ=255)
         cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                                              <  1:1012803/71‑3 (MQ=255)
              tttACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCaaa                                                                     >  1:248902/3‑59 (MQ=255)
              tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                                      <  1:1554937/58‑1 (MQ=255)
              tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:1714441/57‑1 (MQ=255)
              tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:416130/57‑1 (MQ=255)
              tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          >  1:1434260/1‑70 (MQ=255)
               cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                        >  1:2402462/1‑71 (MQ=255)
               cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                                         >  1:3994567/1‑70 (MQ=255)
                tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaaa                                                                     >  1:119096/1‑57 (MQ=255)
                tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                       >  1:133312/1‑71 (MQ=255)
                tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                       >  1:3220899/1‑71 (MQ=255)
                  cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                       <  1:4080024/69‑1 (MQ=255)
                  cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                     <  1:3668258/71‑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)
                          tGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                                       <  1:2116866/45‑1 (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:2814993/1‑71 (MQ=255)
                             cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATg                                          <  1:239614/71‑1 (MQ=255)
                              ggTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGc                                         <  1:3009714/71‑1 (MQ=255)
                              ggTTTACTCCGACCATTTCAGCCATGTTACGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:3312450/66‑1 (MQ=255)
                               gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCa                                        <  1:3919056/71‑1 (MQ=255)
                                tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTaa                                                           <  1:1068004/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:2338357/1‑70 (MQ=255)
                                  tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGc                                                >  1:2947690/1‑60 (MQ=255)
                                  tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                       >  1:3680293/1‑69 (MQ=255)
                                  tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                     >  1:1855529/1‑71 (MQ=255)
                                   aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                     <  1:253048/70‑1 (MQ=255)
                                       cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          >  1:2717872/1‑45 (MQ=255)
                                       cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                     >  1:2593160/1‑50 (MQ=255)
                                       cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                               >  1:3403223/1‑56 (MQ=255)
                                           cATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                        >  1:251483/1‑43 (MQ=255)
                                           cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCa                                        >  1:3276909/1‑59 (MQ=255)
                                            aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTa                               <  1:1667984/67‑1 (MQ=255)
                                             tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                          <  1:2863215/39‑1 (MQ=255)
                                             tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATg                                                      <  1:2207404/43‑1 (MQ=255)
                                             tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                          <  1:679027/71‑1 (MQ=255)
                                               tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGa                         <  1:2342399/70‑1 (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:795084/71‑1 (MQ=255)
                                                   ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    <  1:3976246/71‑1 (MQ=255)
                                                   ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    <  1:630160/71‑1 (MQ=255)
                                                   ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    <  1:423590/71‑1 (MQ=255)
                                                   ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    <  1:925931/71‑1 (MQ=255)
                                                    cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                                 >  1:4082641/1‑57 (MQ=255)
                                                    cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                   <  1:3087011/71‑1 (MQ=255)
                                                    cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                   <  1:3683127/71‑1 (MQ=255)
                                                    cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg                   <  1:775427/71‑1 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    >  1:3613364/1‑69 (MQ=255)
                                                     aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    >  1:2265557/1‑69 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:2223298/42‑1 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:419792/42‑1 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:295266/42‑1 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                              <  1:1544684/42‑1 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                      <  1:824364/66‑1 (MQ=255)
                                                      tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                     <  1:2066690/67‑1 (MQ=255)
                                                      tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCttt                                <  1:2911709/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)
                                                        tttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                            >  1:148076/1‑42 (MQ=255)
                                                        tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                              >  1:2445816/3‑56 (MQ=255)
                                                        tctCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAt                              >  1:1880178/3‑56 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                            >  1:899206/1‑40 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                     >  1:3582595/1‑47 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTa                                     >  1:3041325/1‑47 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                    >  1:3999649/1‑64 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt                     >  1:3111456/1‑63 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg              <  1:3833739/70‑1 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg              <  1:594621/70‑1 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg              <  1:3936350/70‑1 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg              <  1:3462416/70‑1 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggg               >  1:1463137/1‑69 (MQ=255)
                                                          tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                  >  1:1228556/1‑66 (MQ=255)
                                                          tCGGAGATGCGGAAAGTGCTGCTAAGAATGACCAg                                                 >  1:2288865/1‑35 (MQ=255)
                                                           cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg             <  1:1839896/70‑1 (MQ=255)
                                                            ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                            >  1:3888278/1‑38 (MQ=255)
                                                              tGATGCGGAAAGTGCTGCTAAGAATGACCATCTGGATGCATTATCTTTATTTTTATCa                      <  1:2245708/58‑1 (MQ=255)
                                                                 tGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgag                 <  1:409267/60‑1 (MQ=255)
                                                                  gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATc                                   >  1:1592872/1‑41 (MQ=255)
                                                                     gAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt  >  1:295657/1‑71 (MQ=255)
                                                                     |                                                                      
CTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAT  >  minE/2571109‑2571248

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: