Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 608,142 C→T 5.7% Q704Q (CAG→CAA)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE608,1420CT 5.7% 301.6 / 5.0 88Q704Q (CAG→CAA) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (57/26);  new base T (2/3);  total (59/29)
Fisher's exact test for biased strand distribution p-value = 3.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGA  >  minE/608073‑608211
                                                                     |                                                                     
cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                                      <  1:2993827/71‑1 (MQ=255)
 ggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTg                                                                     <  1:808871/71‑1 (MQ=12)
  gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                                                    >  1:87146/1‑71 (MQ=255)
   atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                                                   >  1:2458871/1‑71 (MQ=255)
        tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                                                      >  1:3168975/1‑63 (MQ=255)
        tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                                               >  1:747842/1‑70 (MQ=255)
           ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                                            <  1:2698739/70‑1 (MQ=255)
           ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                                            <  1:1257222/70‑1 (MQ=255)
            cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacc                                                                  <  1:1602284/63‑1 (MQ=255)
              ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                        <  1:2205812/71‑1 (MQ=255)
              ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                        <  1:757415/71‑1 (MQ=255)
              ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                        <  1:662534/71‑1 (MQ=255)
               gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                        >  1:1728911/1‑70 (MQ=255)
                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                                    >  1:1585647/1‑71 (MQ=255)
                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                                    >  1:319676/1‑71 (MQ=255)
                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                                                     >  1:3293199/1‑70 (MQ=255)
                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                        <  1:1791464/67‑1 (MQ=255)
                   tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                                                   >  1:994847/1‑71 (MQ=255)
                    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                                                  >  1:1545759/1‑71 (MQ=255)
                         gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                                                          >  1:1281253/1‑58 (MQ=255)
                         gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                                             >  1:789870/1‑71 (MQ=255)
                           gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                           >  1:3071801/1‑71 (MQ=255)
                              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                      >  1:2311773/1‑57 (MQ=255)
                              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                                             >  1:1126238/1‑66 (MQ=255)
                              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                         >  1:1333066/1‑70 (MQ=255)
                              gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg                                         >  1:3085286/1‑70 (MQ=255)
                                  agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa                                      >  1:591115/1‑68 (MQ=12)
                                   gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                      >  1:856523/1‑52 (MQ=255)
                                   gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                      >  1:58018/1‑52 (MQ=255)
                                   gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                      >  1:65986/1‑52 (MQ=255)
                                   gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa                                   <  1:91753/71‑1 (MQ=255)
                                      tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                                            <  1:2349467/43‑1 (MQ=255)
                                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                                                 >  1:511709/1‑36 (MQ=255)
                                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                                                 >  1:1106761/1‑36 (MQ=255)
                                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                                    >  1:3546612/1‑49 (MQ=255)
                                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                                                    >  1:3355235/1‑49 (MQ=255)
                                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                           >  1:1086844/1‑58 (MQ=255)
                                        tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                              >  1:1453375/1‑71 (MQ=255)
                                         cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                              >  1:924633/1‑70 (MQ=255)
                                           cGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                           <  1:2671877/71‑1 (MQ=255)
                                            ggctcatcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                           <  1:3003433/67‑1 (MQ=255)
                                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                                                      >  1:581999/1‑43 (MQ=255)
                                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                           >  1:1601783/1‑70 (MQ=255)
                                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                           >  1:1273188/1‑70 (MQ=255)
                                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                           >  1:354624/1‑70 (MQ=255)
                                            ggTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatccat                           >  1:2027532/1‑70 (MQ=255)
                                             gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                         <  1:585938/71‑1 (MQ=255)
                                             gTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                         <  1:1962647/71‑1 (MQ=255)
                                               tcatcaCGTTAACGTTCAGGTGCTGACCACCTTCgatg                                                        >  1:2193590/1‑38 (MQ=255)
                                               tcatcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc                              >  1:1791889/1‑64 (MQ=255)
                                               tcatcaCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCCTCGTGGTGGAAGTAAc                                  >  1:3146625/1‑60 (MQ=255)
                                                catcaCGTTAAAGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGa                       <  1:2986580/70‑1 (MQ=255)
                                                  ttacgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccg                    <  1:1481784/69‑2 (MQ=255)
                                                   caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGaa                                       >  1:1169785/1‑51 (MQ=255)
                                                   caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAcc                                 >  1:400568/1‑57 (MQ=255)
                                                   caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                   >  1:2253541/1‑71 (MQ=255)
                                                   caCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                   >  1:4059787/1‑71 (MQ=255)
                                                    acgttgacgttgAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGccga                   >  1:825321/1‑68 (MQ=9)
                                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                         >  1:3804569/1‑64 (MQ=255)
                                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                  >  1:304540/1‑71 (MQ=255)
                                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagc                  >  1:3863955/1‑71 (MQ=255)
                                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAc                      >  1:126087/1‑67 (MQ=255)
                                                    aCGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATAAGAccagc                  >  1:1568640/1‑71 (MQ=255)
                                                    aCGTTAACGTTCAGGTGCTGAACACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccag                   <  1:570405/70‑1 (MQ=255)
                                                      gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                <  1:3344541/71‑1 (MQ=255)
                                                      gTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagcca                <  1:2048259/71‑1 (MQ=255)
                                                       ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                           >  1:3926271/1‑43 (MQ=255)
                                                       ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                           >  1:3071311/1‑43 (MQ=255)
                                                       ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                           >  1:595676/1‑43 (MQ=255)
                                                       ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt                                           >  1:806294/1‑43 (MQ=255)
                                                       ttAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGAccagccag               <  1:1042484/71‑1 (MQ=255)
                                                        tAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAgg              >  1:560480/1‑71 (MQ=255)
                                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                         >  1:3458225/1‑59 (MQ=255)
                                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCa                         >  1:1012167/1‑59 (MQ=255)
                                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt             >  1:4110437/1‑71 (MQ=255)
                                                         aaCGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt             >  1:2299336/1‑71 (MQ=255)
                                                           cGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTg           >  1:2453846/1‑71 (MQ=255)
                                                            gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGa                       >  1:290658/1‑58 (MQ=255)
                                                            gTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTgg          <  1:2637070/71‑1 (MQ=255)
                                                             ttgaGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGAAATACCCATCCAGCAGGCCGACAAGGtt            <  1:1091981/65‑1 (MQ=255)
                                                             ttCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGt         <  1:3293469/71‑1 (MQ=255)
                                                               cAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt       >  1:1531091/1‑71 (MQ=255)
                                                                aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCtt      <  1:88109/71‑1 (MQ=255)
                                                                aGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCt       <  1:612462/70‑1 (MQ=255)
                                                                 ggTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTa     <  1:1226099/71‑1 (MQ=255)
                                                                  gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGt             <  1:1007309/62‑1 (MQ=255)
                                                                  gTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTAc    >  1:1617533/1‑71 (MQ=255)
                                                                   tGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACg   <  1:1012123/71‑1 (MQ=255)
                                                                     cTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGa  <  1:3227277/70‑1 (MQ=255)
                                                                     |                                                                     
CGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATCCATCAGACCAGCCAGGTTGGTCTTACGA  >  minE/608073‑608211

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: ‑