Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,435,202 C→A 40.6% I273I (ATC→ATA gpp → guanosine pentaphosphatase/exopolyphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,435,2020CA40.6% 54.8 / 78.8 91I273I (ATC→ATAgppguanosine pentaphosphatase/exopolyphosphatase
Reads supporting (aligned to +/- strand):  ref base C (11/43);  new base A (13/24);  total (24/67)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-01

GCTGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGAAGGC  >  minE/2435151‑2435265
                                                   |                                                               
gCTGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGaa                                                  >  1:2097687/1‑67 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:90493/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:458483/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:3364216/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:3281738/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2114437/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2212132/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2901282/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2750174/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2265379/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2533991/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:2363664/69‑1 (MQ=255)
  tGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGa                                              <  1:1702594/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:274315/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:1966011/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:2407040/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:2498260/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:110935/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:2953313/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:2962792/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:3075162/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:999470/69‑1 (MQ=255)
   gACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:1170271/69‑1 (MQ=255)
    aCGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGaa                                             <  1:3237225/68‑1 (MQ=255)
        tGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTc                                        <  1:3022588/69‑1 (MQ=255)
         ggAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCa                                       <  1:104897/69‑1 (MQ=255)
         ggAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCa                                       <  1:2755386/69‑1 (MQ=255)
            aCGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:1915090/69‑1 (MQ=255)
            aCGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:3425799/69‑1 (MQ=255)
            aCGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:1876605/69‑1 (MQ=255)
            aCGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:602025/69‑1 (MQ=255)
            aCGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:224257/69‑1 (MQ=255)
             cGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:2515295/68‑1 (MQ=255)
                   ttAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGAcc                             >  1:1828452/1‑69 (MQ=255)
                     aGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCAtttt                                                         >  1:2391960/1‑39 (MQ=255)
                     aGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCAtt                                                           >  1:2363312/1‑37 (MQ=255)
                     aGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGAccc                            >  1:1623242/1‑68 (MQ=255)
                     aGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGAccc                            >  1:1785928/1‑68 (MQ=255)
                     aGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCt                           >  1:2361211/1‑69 (MQ=255)
                           cccGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:658520/54‑1 (MQ=255)
                           cccGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:118079/54‑1 (MQ=255)
                           cccGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:553485/54‑1 (MQ=255)
                           cccGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:1328188/54‑1 (MQ=255)
                           cccGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:421939/54‑1 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:732092/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:1950760/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:276783/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:2916028/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:1855142/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:2748743/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:1460382/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:1334699/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:1177968/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:1988956/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:3085863/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggcgggcg                   >  1:3049413/1‑69 (MQ=255)
                             cGAGTGGTCTGGCGATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTggc                        >  1:2744338/1‑64 (MQ=255)
                               aGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTc                                        <  1:2598850/46‑1 (MQ=255)
                                  ggTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTgcg              <  1:2522156/69‑1 (MQ=255)
                                     cTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTg           <  1:1905445/69‑1 (MQ=255)
                                          gATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTTCGCTGCGTg       <  1:1820154/68‑1 (MQ=255)
                                          gATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTg       <  1:1941641/68‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:355185/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:436003/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:462266/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:494927/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:55670/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:3340298/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:903323/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1188506/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1453556/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1712059/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1857639/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1884836/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:2099863/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:215377/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1760188/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1756518/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1547956/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:2539072/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:2618071/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:2672268/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:2702187/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:2751870/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1388690/69‑1 (MQ=255)
                                          gATCCTGATAGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGa      <  1:1215525/69‑1 (MQ=255)
                                              cTGATCGCCATTTTTACCGAACTGAATATTCAgtg                                    <  1:3222146/35‑1 (MQ=255)
                                              cTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGAAGGc  >  1:1742453/1‑69 (MQ=255)
                                              cTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGAAGGc  >  1:2785994/1‑69 (MQ=255)
                                              cTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGAAGGc  >  1:981233/1‑69 (MQ=255)
                                              cTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGAAGGc  >  1:997579/1‑69 (MQ=255)
                                                   |                                                               
GCTGACGCTGGAACGTGCGTTAGTGTTCCCGAGTGGTCTGGCGATCCTGATCGCCATTTTTACCGAACTGAATATTCAGTGTATGACCCTGGCGGGCGGTGCGCTGCGTGAAGGC  >  minE/2435151‑2435265

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: