Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,221,878 A→G 5.9% Q151Q (CAA→CAG yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,221,8780AG 5.9% 293.6 / 8.8 101Q151Q (CAA→CAGyebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (55/40);  new base G (2/4);  total (57/44)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01

Only 100 of 107 total aligned reads displayed.

GGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAAC  >  minE/1221810‑1221947
                                                                    |                                                                     
gggAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAAcc                                                                     >  1:1015567/1‑71 (MQ=255)
   aaCTATATCTGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGa                                                                   >  1:3709790/1‑70 (MQ=255)
   aaCTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaa                                                                  <  1:3085703/71‑1 (MQ=255)
     cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                <  1:351015/71‑1 (MQ=255)
     cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                <  1:3931379/71‑1 (MQ=255)
     cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                <  1:2595770/71‑1 (MQ=255)
        tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAatat                                                              >  1:2924238/1‑70 (MQ=255)
        tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATc                                                             <  1:2798219/71‑1 (MQ=255)
        tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATc                                                             >  1:1528017/1‑71 (MQ=255)
          tCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCgg                                                           >  1:2140211/1‑71 (MQ=255)
             gCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                        >  1:70409/1‑71 (MQ=255)
               atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                      <  1:2802799/71‑2 (MQ=38)
               atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                      <  1:3955505/71‑2 (MQ=38)
               atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                      <  1:2652691/71‑2 (MQ=38)
               atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                      <  1:250044/71‑2 (MQ=38)
                tgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCaa                                                                       <  1:2098125/53‑1 (MQ=255)
                tgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCaa                                                                       <  1:3212438/53‑1 (MQ=255)
                tgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                      >  1:2012182/1‑70 (MQ=255)
                 gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                    >  1:906849/1‑68 (MQ=37)
                  atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:2725512/1‑71 (MQ=255)
                  atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:802471/1‑71 (MQ=255)
                  atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:3852980/1‑71 (MQ=255)
                  atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:1029857/1‑71 (MQ=255)
                         gtaaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                            >  1:4177585/1‑71 (MQ=255)
                          taaaggtaaagAGCAGTATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                              <  1:4320184/68‑1 (MQ=255)
                          taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                <  1:2924015/50‑1 (MQ=255)
                          taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                              <  1:3412172/68‑1 (MQ=255)
                          taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                              <  1:3559872/68‑1 (MQ=255)
                          taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATc                                           >  1:2849333/1‑71 (MQ=255)
                            aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                         >  1:3292777/1‑55 (MQ=255)
                            aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                         >  1:3794254/1‑55 (MQ=255)
                            aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                         >  1:4073706/1‑55 (MQ=255)
                            aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCt                                          <  1:4350618/70‑1 (MQ=255)
                             aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                        >  1:3657989/1‑71 (MQ=255)
                             aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                        >  1:3089539/1‑71 (MQ=255)
                             aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                        >  1:3518189/1‑71 (MQ=255)
                              ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                      >  1:2567144/1‑56 (MQ=255)
                              ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                      >  1:1809656/1‑56 (MQ=255)
                              ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                       <  1:3418502/71‑1 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaaa                                                                 >  1:440773/1‑44 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaaa                                                                 >  1:1159695/1‑44 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                        <  1:2656646/53‑1 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                        <  1:859358/53‑1 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                        <  1:4046307/53‑1 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:2414956/1‑58 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgatg                                      <  1:1226216/71‑1 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgatg                                      <  1:452570/71‑1 (MQ=255)
                               gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                       >  1:4003162/1‑70 (MQ=255)
                                taaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGAGATCtgatg                                      >  1:2503251/1‑70 (MQ=255)
                                 aaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                            >  1:3519554/1‑63 (MQ=255)
                                  aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                            >  1:3959734/1‑62 (MQ=255)
                                  aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                            >  1:3108247/1‑62 (MQ=255)
                                  aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATc                                   >  1:3429980/1‑71 (MQ=255)
                                   agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                            <  1:2750500/45‑1 (MQ=255)
                                   agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                            <  1:868985/45‑1 (MQ=255)
                                   agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                            <  1:1423413/45‑1 (MQ=255)
                                      gCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCAc                                <  1:2187188/70‑1 (MQ=255)
                                       cAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCt                              <  1:1151494/71‑1 (MQ=255)
                                       cAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCt                              <  1:1546748/71‑1 (MQ=255)
                                        aGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTg                             >  1:3115701/1‑71 (MQ=255)
                                         ggATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGc                            <  1:3742953/71‑1 (MQ=255)
                                           aTCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGc                            <  1:2464709/69‑1 (MQ=255)
                                                tttGTCGCACTCGATACCCAACCGAAATATCGGCTGGAc                                                     <  1:1780079/39‑1 (MQ=255)
                                                tttGTCGCACTCGATACCCAACCGAAATATCGGCTGGAc                                                     <  1:806404/39‑1 (MQ=255)
                                                  tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:1275829/1‑39 (MQ=255)
                                                  tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                   >  1:3476280/1‑39 (MQ=255)
                                                  tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                            >  1:1969534/1‑46 (MQ=255)
                                                  tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGccc                      <  1:1301739/68‑1 (MQ=255)
                                                      gCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct               <  1:1955496/71‑1 (MQ=255)
                                                       cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc              >  1:191557/1‑71 (MQ=255)
                                                       cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc              >  1:3299012/1‑71 (MQ=255)
                                                       cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc              <  1:1100457/71‑1 (MQ=255)
                                                         cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg            >  1:77109/1‑71 (MQ=255)
                                                         cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg            <  1:2811734/71‑1 (MQ=255)
                                                         cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg            >  1:2097516/1‑71 (MQ=255)
                                                          tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                       <  1:3488793/43‑1 (MQ=255)
                                                          tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                       <  1:3575275/43‑1 (MQ=255)
                                                           cGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCAc                                >  1:2631649/1‑49 (MQ=255)
                                                           cGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGtt          >  1:2608283/1‑71 (MQ=255)
                                                            gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAg                         >  1:2740496/1‑55 (MQ=255)
                                                            gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc         >  1:4136526/1‑71 (MQ=255)
                                                            gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc         <  1:2952791/71‑1 (MQ=255)
                                                            gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc         >  1:3863571/1‑71 (MQ=255)
                                                            gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc         >  1:1348553/1‑71 (MQ=255)
                                                            gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc         >  1:1940407/1‑71 (MQ=255)
                                                              tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                    <  1:1265361/42‑1 (MQ=255)
                                                              tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                    <  1:3089195/42‑1 (MQ=255)
                                                              tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                           >  1:741907/1‑51 (MQ=255)
                                                              tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                           >  1:1649786/1‑51 (MQ=255)
                                                                cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCt      >  1:614977/1‑70 (MQ=255)
                                                                cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTg     <  1:4404829/71‑1 (MQ=255)
                                                                 ccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct               <  1:2667694/60‑1 (MQ=255)
                                                                 ccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct               <  1:1911608/60‑1 (MQ=255)
                                                                  cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa   >  1:2990037/1‑71 (MQ=255)
                                                                  cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa   >  1:2272879/1‑71 (MQ=255)
                                                                  cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa   >  1:1774182/1‑71 (MQ=255)
                                                                  cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa   >  1:3159878/1‑71 (MQ=255)
                                                                   aaCCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAAc  <  1:2911884/71‑1 (MQ=255)
                                                                    aCCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa   >  1:3417214/1‑69 (MQ=255)
                                                                    aCCGAAATATAGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa   >  1:1017712/1‑69 (MQ=255)
                                                                    |                                                                     
GGGAACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAAC  >  minE/1221810‑1221947

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: