Predicted mutation
evidence seq id position mutation annotation gene description
RA minE 2,349,606 A→G V681A (GTG→GCG)  polA ← fused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,349,6060AG99.0% 363.5 / ‑5.1 98V681A (GTG→GCG) polAfused DNA polymerase I 5'‑>3' exonuclease, 3'‑>5' polymerase and 3'‑>5' exonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (43/54);  minor base T (0/1);  total (43/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CTGAGACAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCACCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGTGA  >  minE/2349538‑2349671
                                                                    |                                                                 
cTGAGACAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGcc                                                                 <  1:174099/71‑1 (MQ=255)
 tGAGACAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCg                                                                >  1:2093793/1‑71 (MQ=255)
  gagaCAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCg                                                                >  1:4083120/1‑70 (MQ=255)
  gagaCAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCg                                                                >  1:1838103/1‑70 (MQ=255)
   agaCAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGCCCTTCTTCGTTACGCGCCgg                                                               >  1:2404618/1‑70 (MQ=255)
   agaCAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCgg                                                               <  1:3023012/70‑1 (MQ=255)
      cAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATg                                                           <  1:1735836/71‑1 (MQ=255)
       aaTCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGt                                                          <  1:3012005/71‑1 (MQ=255)
        aTCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGtt                                                         <  1:1351169/71‑1 (MQ=255)
        aTCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGtt                                                         <  1:2277402/71‑1 (MQ=255)
         tCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGttt                                                        >  1:771981/1‑71 (MQ=255)
          cacaTAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGtttt                                                       >  1:803833/1‑71 (MQ=255)
           acaTAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTg                                                      <  1:3069364/71‑1 (MQ=255)
             aTAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCa                                                    <  1:1297233/71‑1 (MQ=255)
               aaTCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAgg                                                  >  1:3039969/1‑71 (MQ=255)
                aTCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGt                                                 <  1:356985/71‑1 (MQ=255)
                aTCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGt                                                 <  1:378526/71‑1 (MQ=255)
                   ctctGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAg                                              <  1:3279525/71‑1 (MQ=255)
                   ctctGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAg                                              <  1:281207/71‑1 (MQ=255)
                   ctctGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAg                                              <  1:366270/71‑1 (MQ=255)
                     ctGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGtt                                                         >  1:939625/1‑58 (MQ=255)
                     ctGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGtt                                                         >  1:1893664/1‑58 (MQ=255)
                        gcgcAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGAt                                           <  1:4031059/69‑1 (MQ=255)
                           cAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGtt                                      <  1:2853623/71‑1 (MQ=255)
                            aaTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCgg                                        >  1:2857255/1‑68 (MQ=255)
                            aaTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGtt                                      >  1:845551/1‑70 (MQ=255)
                            aaTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGtt                                      >  1:1442526/1‑70 (MQ=255)
                            aaTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTg                                     <  1:910742/71‑1 (MQ=255)
                            aaTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTg                                     <  1:2059245/71‑1 (MQ=255)
                            aaTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTg                                     <  1:79344/71‑1 (MQ=255)
                             aTAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGtt                                                <  1:1467880/59‑1 (MQ=255)
                              tAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGAc                                   >  1:906947/1‑71 (MQ=255)
                              tAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGAc                                   >  1:3661278/1‑71 (MQ=255)
                                aaCGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTa                                               >  1:1447860/1‑57 (MQ=255)
                                aaCGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACg                                  <  1:2728252/70‑1 (MQ=255)
                                aaCGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACg                                  <  1:3634801/70‑1 (MQ=255)
                                aaCGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGa                                 >  1:1191845/1‑71 (MQ=255)
                                aaCGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGa                                 <  1:3624397/71‑1 (MQ=255)
                                aaCGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGa                                 <  1:2272169/71‑1 (MQ=255)
                                 aCGCCTAGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGAt                                <  1:1842924/71‑1 (MQ=255)
                                  cGCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATa                               <  1:864484/71‑1 (MQ=255)
                                   gCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATa                               >  1:461086/1‑70 (MQ=255)
                                   gCCTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATa                               <  1:62174/70‑1 (MQ=255)
                                     cTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATaaa                             >  1:673832/1‑70 (MQ=255)
                                     cTGGCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAAc                            >  1:2540842/1‑71 (MQ=255)
                                       ggCGGATACGACGACCTTCTTCGTTACGCGCCGGGATGTTTTGCAGGTTAGGAt                                           >  1:36408/1‑54 (MQ=39)
                                       ggCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGAt                                           >  1:4050374/1‑54 (MQ=255)
                                       ggCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAAc                            >  1:99491/1‑69 (MQ=255)
                                       ggCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGt                          >  1:3309690/1‑71 (MQ=255)
                                        gCGGATACGACGCCCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTc                         >  1:1052747/1‑71 (MQ=255)
                                        gCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATaa                              <  1:1168039/66‑1 (MQ=255)
                                        gCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACg                           <  1:200780/69‑1 (MQ=255)
                                        gCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACg                           <  1:2474006/69‑1 (MQ=255)
                                        gCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTc                         <  1:1458717/71‑1 (MQ=255)
                                        gCGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTc                         >  1:2140502/1‑71 (MQ=255)
                                         cGGATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTcc                        >  1:1448794/1‑71 (MQ=255)
                                          ggATACGACGACCTTCTTCGTTACGCGCCGGAATGttt                                                        >  1:3427301/1‑38 (MQ=255)
                                          ggATACGACGACCTTCTTCGTTACGCGCCGGAATGttt                                                        >  1:3753186/1‑38 (MQ=255)
                                          ggATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGt                                                 <  1:3110253/45‑1 (MQ=255)
                                           gATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCg                      <  1:2147472/71‑1 (MQ=255)
                                           gATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCg                      <  1:2056336/71‑1 (MQ=255)
                                           gATACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCg                      <  1:2983761/71‑1 (MQ=255)
                                             tACGACGACCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGtt                    <  1:3945946/71‑1 (MQ=255)
                                              acgacgACCTTCTTCGTTACGCGCCGGAATGTTTTGGAGGTTAGGATCGGTTGACGATAAACGTCCCGTTg                   >  1:837650/1‑71 (MQ=255)
                                               cgacgaCCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTg                   <  1:1268852/70‑1 (MQ=255)
                                               cgacgaCCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTg                   <  1:377654/70‑1 (MQ=255)
                                               cgacgaCCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGc                  <  1:2069812/71‑1 (MQ=255)
                                                gacgacCTTCTTCGTTACGCTCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGAt                                <  1:3457239/56‑1 (MQ=255)
                                                gacgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCg                                         <  1:2242754/47‑1 (MQ=255)
                                                gacgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTg                   >  1:291859/1‑69 (MQ=255)
                                                gacgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCa                 <  1:2383162/71‑1 (MQ=255)
                                                 acgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAAc                            >  1:1885003/1‑59 (MQ=255)
                                                 acgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAAc                            >  1:341595/1‑59 (MQ=255)
                                                  cgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAg                >  1:3544794/1‑70 (MQ=255)
                                                  cgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGt               >  1:2909852/1‑71 (MQ=255)
                                                  cgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGt               >  1:455152/1‑71 (MQ=255)
                                                  cgacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGt               >  1:3640486/1‑71 (MQ=255)
                                                   gacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGt                     <  1:3680178/64‑1 (MQ=255)
                                                   gacCTTCTTCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGtt              <  1:2196373/71‑1 (MQ=255)
                                                      cttcttCGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACt           >  1:2046155/1‑71 (MQ=255)
                                                       ttcttcGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTg          >  1:4032738/1‑71 (MQ=255)
                                                       ttcttcGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTg          >  1:151426/1‑71 (MQ=255)
                                                         cttcGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTg          <  1:3853243/69‑1 (MQ=255)
                                                         cttcGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGcc        >  1:3478680/1‑71 (MQ=255)
                                                          ttcGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGtt                                      <  1:1322514/40‑1 (MQ=255)
                                                          ttcGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGa                                    >  1:3755353/1‑42 (MQ=255)
                                                            cGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGcc        <  1:3039001/68‑1 (MQ=255)
                                                            cGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGcc        <  1:4048005/68‑1 (MQ=255)
                                                            cGTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCt       <  1:1312538/69‑1 (MQ=255)
                                                             gTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGcc        <  1:205866/67‑1 (MQ=255)
                                                             gTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGcc        <  1:3127172/67‑1 (MQ=255)
                                                             gTTACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGt    <  1:3938813/71‑1 (MQ=255)
                                                              ttACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATaa                              <  1:447652/44‑1 (MQ=255)
                                                              ttACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATaa                              <  1:3060805/44‑1 (MQ=255)
                                                              ttACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGTg   <  1:3518342/71‑1 (MQ=255)
                                                              ttACGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGTg   <  1:854743/71‑1 (MQ=255)
                                                                aCGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGTGa  >  1:2946954/1‑70 (MQ=255)
                                                                aCGCGCCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGTGa  >  1:2160204/1‑70 (MQ=255)
                                                                    |                                                                 
CTGAGACAATCACATAATCCTCTGGCGCAATAAACGCCTGGCGGATACGACGACCTTCTTCGTTACGCACCGGAATGTTTTGCAGGTTAGGATCGGTTGACGATAAACGTCCCGTTGCAGTTACTGCCTGGTGA  >  minE/2349538‑2349671

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: