Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 705,295 A→G 18.4% G32G (GGT→GGC ymcB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE705,2950AG18.4% 220.0 / 43.6 98G32G (GGT→GGCymcBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/45);  new base G (7/11);  total (42/56)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.83e-01

Only 100 of 104 total aligned reads displayed.

CGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTG  >  minE/705225‑705361
                                                                      |                                                                  
cgcgCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCa                                                                    >  1:2746375/1‑71 (MQ=255)
       ttGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCa                                                             <  1:3472183/71‑1 (MQ=255)
        tGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCa                                                             <  1:1203374/70‑1 (MQ=255)
             ttGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGAt                                                          <  1:2274008/68‑1 (MQ=255)
             ttGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATaaa                                                       <  1:1282454/71‑1 (MQ=255)
              tgtgTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATaaa                                                       >  1:2645319/1‑70 (MQ=255)
               gtgtCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCAcc                                                                  >  1:1241686/1‑58 (MQ=255)
               gtgtCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATaaa                                                       >  1:3528927/1‑69 (MQ=255)
                tgtCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATa                                                     <  1:783976/70‑1 (MQ=255)
                 gtCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGt                                                   >  1:872385/1‑71 (MQ=255)
                  tCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGt                                                   >  1:1941369/1‑70 (MQ=255)
                  tCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGt                                                   <  1:3597261/70‑1 (MQ=255)
                    aCCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCa                                                 >  1:1382098/1‑70 (MQ=255)
                    aCCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCAc                                                >  1:3686727/1‑71 (MQ=255)
                    aCCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGt                                                   >  1:1456012/1‑68 (MQ=255)
                        aCTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGAt                                                          <  1:4073570/57‑1 (MQ=255)
                          tGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCa                                                                    >  1:1482620/1‑45 (MQ=255)
                          tGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCa                                                                    >  1:1531475/1‑45 (MQ=255)
                          tGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCa                                                             <  1:2504169/52‑1 (MQ=255)
                           ggACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCa                                                                    >  1:3513257/1‑44 (MQ=255)
                           ggACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAg                                                    >  1:1749820/1‑60 (MQ=255)
                           ggACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAg                                                    >  1:1895063/1‑60 (MQ=255)
                           ggACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCt                                         >  1:169064/1‑71 (MQ=255)
                            gACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCtt                                        >  1:2090653/1‑71 (MQ=255)
                             acaacaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCAc                                                                   <  1:672483/43‑1 (MQ=255)
                              caacaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTc                                           <  1:1615723/66‑1 (MQ=255)
                              caacaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgc                                      <  1:3359878/71‑1 (MQ=255)
                              caacaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgc                                      <  1:4028444/71‑1 (MQ=255)
                                acaTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATa                                                         <  1:3074285/50‑1 (MQ=255)
                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                   <  1:755283/71‑1 (MQ=255)
                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                   <  1:298600/71‑1 (MQ=255)
                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTCCGCa                                   <  1:2123715/71‑1 (MQ=255)
                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgcgc                                    >  1:2574141/1‑70 (MQ=255)
                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgcgc                                    >  1:3126149/1‑70 (MQ=255)
                                 caTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgcgc                                    >  1:1125215/1‑70 (MQ=255)
                                  aTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGa                                                           <  1:898051/46‑1 (MQ=255)
                                  aTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTgcgc                                    >  1:1017828/1‑69 (MQ=255)
                                  aTTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                   >  1:1012953/1‑70 (MQ=255)
                                   ttttCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATaaa                                                       <  1:1108265/49‑1 (MQ=255)
                                    tttCCACGGGGCCAACGGAAAGTGTCAGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAg                                >  1:684618/1‑71 (MQ=255)
                                     ttCCACGGGGCCAACGGAAAGTGTCTGTTGTTCAc                                                                   >  1:3346676/1‑35 (MQ=255)
                                      tCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATa                                                     >  1:346573/1‑48 (MQ=255)
                                      tCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATa                                                     >  1:3220996/1‑48 (MQ=255)
                                      tCCACGGGGCCAACGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAgg                               <  1:728139/70‑1 (MQ=255)
                                      tCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAgg                               <  1:2232568/70‑1 (MQ=255)
                                       ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATaa                                                        <  1:690886/44‑1 (MQ=255)
                                       ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAgg                               <  1:1736904/69‑1 (MQ=255)
                                       ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAgg                               <  1:528228/69‑1 (MQ=255)
                                       ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                              <  1:4033203/70‑1 (MQ=255)
                                       ccACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                              <  1:245805/70‑1 (MQ=255)
                                             ggCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGgt                       <  1:3229876/71‑1 (MQ=255)
                                              gCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAAt                                                      >  1:1264203/1‑39 (MQ=255)
                                                  aCGTAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCtt                                        <  1:444375/49‑1 (MQ=255)
                                                  aCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCtt                                        <  1:1255738/49‑1 (MQ=255)
                                                  aCGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCAt                  <  1:152846/71‑1 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATg                           <  1:546343/61‑1 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATg                           <  1:443937/61‑1 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                   <  1:2438729/53‑1 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCa                             <  1:2648695/59‑1 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCa                             <  1:81604/59‑1 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCAt                  >  1:3084002/1‑70 (MQ=255)
                                                   cGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCAt                  >  1:4024565/1‑70 (MQ=255)
                                                    ggAAAGTGTCTGTTGTTGACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTa                >  1:3603196/1‑71 (MQ=255)
                                                    ggAAAGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCa                             >  1:471474/1‑58 (MQ=255)
                                                    ggAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTa                >  1:3684149/1‑71 (MQ=255)
                                                     gAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCtt                                        <  1:527600/46‑1 (MQ=255)
                                                      aaaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                              >  1:1347806/1‑55 (MQ=255)
                                                      aaaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGc                              >  1:1325549/1‑55 (MQ=255)
                                                      aaaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATggg                         <  1:2125473/60‑1 (MQ=255)
                                                       aaGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACg              >  1:3246406/1‑70 (MQ=255)
                                                        aGTGTCTGTTGTTCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                   <  1:3751287/48‑1 (MQ=255)
                                                        aGTGTCTGTTGTTCACCAGGCATATAAATAGTCACCGTTCCTTGCGCAAAGGCAt                            <  1:2296055/55‑1 (MQ=255)
                                                        aGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCAt                            <  1:1882768/55‑1 (MQ=255)
                                                        aGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCAt                            <  1:2713888/55‑1 (MQ=255)
                                                        aGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCAt                            <  1:2420925/55‑1 (MQ=255)
                                                        aGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTa            <  1:133461/71‑1 (MQ=255)
                                                          tgtCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAgg                               <  1:1947606/50‑1 (MQ=255)
                                                           gtCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTa            >  1:3794354/1‑68 (MQ=255)
                                                           gtCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTa            >  1:1330468/1‑68 (MQ=255)
                                                            tCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCAt                            <  1:2009659/51‑1 (MQ=255)
                                                            tCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTaaa          <  1:432375/69‑1 (MQ=255)
                                                            tCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAg         <  1:3048419/70‑1 (MQ=255)
                                                            tCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAg         <  1:1117459/70‑1 (MQ=255)
                                                            tCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAg         <  1:2579639/70‑1 (MQ=255)
                                                              tgttgtTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAg                                >  1:67971/1‑45 (MQ=255)
                                                                ttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTa            <  1:4055154/63‑1 (MQ=255)
                                                                ttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTa            <  1:1514486/63‑1 (MQ=255)
                                                                ttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTa       <  1:2638033/68‑1 (MQ=255)
                                                                ttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTa       <  1:68710/68‑1 (MQ=255)
                                                                ttgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGc    >  1:1769310/1‑71 (MQ=255)
                                                                 tgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATgggg                        >  1:3394658/1‑50 (MQ=255)
                                                                 tgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATgggg                        >  1:2051746/1‑50 (MQ=255)
                                                                 tgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATgggg                        >  1:160171/1‑50 (MQ=255)
                                                                 tgttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCt   >  1:2367881/1‑71 (MQ=255)
                                                                  gttCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTg  >  1:3989922/1‑71 (MQ=255)
                                                                  gttCGCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTg  >  1:3347121/1‑71 (MQ=255)
                                                                   ttCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCa                                   <  1:3083997/37‑1 (MQ=255)
                                                                     cACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCAt                            >  1:664525/1‑42 (MQ=255)
                                                                      aCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGgtgt                     <  1:3997921/48‑1 (MQ=255)
                                                                      aCCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTg  <  1:4029494/67‑1 (MQ=255)
                                                                      |                                                                  
CGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAACGGAAAGTGTCTGTTGTTCACCAGGCAGATAAATAGTCACCGTTCCTTGCGCAAAGGCATGGGGTGTCATTACGTAAAGTACGCTG  >  minE/705225‑705361

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: