Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 607,997 T→A 6.9% T753S (ACT→TCT) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE607,9970TA 6.9% 205.7 / 6.0 72T753S (ACT→TCT) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (44/23);  new base A (2/3);  total (46/26)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

CCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGG  >  minE/607927‑608065
                                                                      |                                                                    
ccTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGt                                                                      <  1:3035322/71‑1 (MQ=255)
 cTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGt                                                                      >  1:4589035/1‑70 (MQ=255)
    tATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGtaat                                                                   >  1:3601005/1‑70 (MQ=255)
    tATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGtaat                                                                   <  1:5612729/70‑1 (MQ=255)
      ttGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGtaata                                                                  >  1:3284179/1‑69 (MQ=255)
           cGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGt                                                              <  1:258612/68‑1 (MQ=255)
           cGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCct                                                           >  1:1132147/1‑71 (MQ=255)
             cTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgc                                                         >  1:5226863/1‑71 (MQ=255)
              tttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGt                                                              >  1:4540037/1‑65 (MQ=255)
               ttttACTGTACGATTTCAGTCAAATCTAATTACCTAGATTGAGTGAAGGTACGAGtaat                                                                   <  1:2383788/59‑1 (MQ=255)
               ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGtaata                                                                  <  1:3605694/60‑1 (MQ=255)
               ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGtaat                                                                   <  1:2282079/59‑1 (MQ=255)
               ttttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                       >  1:2236666/1‑71 (MQ=255)
                tttACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgc                                                      >  1:4449022/1‑71 (MQ=255)
                  tACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgctg                                                    >  1:3002028/1‑71 (MQ=255)
                     tGTACTATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctg                                                       <  1:970631/65‑1 (MQ=255)
                       tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGt                                                                      <  1:3169384/48‑1 (MQ=255)
                       tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGtt                                                  >  1:3642488/1‑68 (MQ=255)
                       tACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCtt                                               >  1:2897703/1‑71 (MQ=255)
                          gATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAg                                            <  1:4346096/71‑1 (MQ=255)
                           aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGt                                                   >  1:1604128/1‑63 (MQ=255)
                           aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGt                                                   >  1:1914988/1‑63 (MQ=255)
                           aTTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGt                                                   >  1:1146529/1‑63 (MQ=255)
                            tttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTc                                          <  1:355191/71‑1 (MQ=255)
                             ttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct                                                     >  1:3999315/1‑59 (MQ=255)
                             ttCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCctgctgct                                                     >  1:2868653/1‑59 (MQ=255)
                                aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGc                                       <  1:2413541/70‑1 (MQ=255)
                                aGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGc                                       <  1:2762980/70‑1 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGtaataa                                                                 >  1:531398/1‑40 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGtaataa                                                                 >  1:3156303/1‑40 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGtaataa                                                                 >  1:2425201/1‑40 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGtaataa                                                                 >  1:1404656/1‑40 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGtaataa                                                                 >  1:3706145/1‑40 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCct                                                           >  1:1269137/1‑46 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                  >  1:2591745/1‑71 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                  >  1:4352840/1‑71 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                  >  1:973397/1‑71 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                  >  1:229397/1‑71 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                  <  1:645509/71‑1 (MQ=255)
                                    aaaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGtt                                  <  1:2151974/71‑1 (MQ=255)
                                     aaTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCtt                                               <  1:894144/57‑1 (MQ=255)
                                      aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTc                                                             >  1:5653643/1‑42 (MQ=255)
                                      aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTc                                                             >  1:380619/1‑42 (MQ=255)
                                      aTCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTc                                                             >  1:919959/1‑42 (MQ=255)
                                      aTCTAATTACATAGATTGACTGAAGGTACGAGTAATAACGt                                                              >  1:1233973/1‑41 (MQ=255)
                                       tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAAc                                                                >  1:964899/1‑38 (MQ=255)
                                       tCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGa                                     <  1:1309261/65‑1 (MQ=255)
                                        cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                              >  1:4477486/1‑71 (MQ=255)
                                        cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                              >  1:982748/1‑71 (MQ=255)
                                        cTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                              >  1:2170650/1‑71 (MQ=255)
                                           aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGaaa                              <  1:4419698/68‑1 (MQ=255)
                                           aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                            >  1:2298041/1‑70 (MQ=255)
                                           aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                            <  1:1962975/70‑1 (MQ=255)
                                           aTTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                           <  1:3026259/71‑1 (MQ=255)
                                            ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                           >  1:4396549/1‑70 (MQ=255)
                                            ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTa                          >  1:2300988/1‑71 (MQ=255)
                                            ttACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGATGCTGTTCTTTAGTCAGCGAGTTGAAACGt                           >  1:3494523/1‑70 (MQ=255)
                                               cATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACt                        <  1:4815956/70‑1 (MQ=255)
                                                aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTa                                             <  1:3725039/48‑1 (MQ=255)
                                                aTAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTa                                             <  1:756690/48‑1 (MQ=255)
                                                aTAGATTGAGTGAAGGTACGAGTAATAAAGTCCTGCTGCTGTTCTTTa                                             <  1:4201361/48‑1 (MQ=255)
                                                 tAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGGCAGCGAGTTGAAACGTACTGc                      <  1:3749529/70‑1 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:71531/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:2338224/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:2561674/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:2244981/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:3197352/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:210154/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGt                                   >  1:176854/1‑56 (MQ=255)
                                                  aGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCg                     >  1:3536577/1‑70 (MQ=255)
                                                     ttGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGc                 <  1:2617571/71‑1 (MQ=255)
                                                      tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                            >  1:1953076/1‑59 (MQ=255)
                                                      tGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                            >  1:1801633/1‑59 (MQ=255)
                                                       gAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGAGAGTTGAAACGTACTGCGTAGcc                <  1:83408/70‑1 (MQ=255)
                                                         gtaaaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGa             <  1:350587/68‑1 (MQ=255)
                                                          tGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGt                           <  1:1294776/56‑1 (MQ=255)
                                                            aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga           <  1:2969174/70‑3 (MQ=255)
                                                            aaGGTACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGaga           <  1:404675/70‑3 (MQ=255)
                                                            aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa           >  1:227011/1‑70 (MQ=255)
                                                            aaGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATa           >  1:5315911/1‑70 (MQ=255)
                                                             aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                            >  1:5010636/1‑52 (MQ=255)
                                                             aGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACg                            >  1:805418/1‑52 (MQ=255)
                                                                tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGacac        >  1:3279435/1‑69 (MQ=255)
                                                                tACGTGAAATAACATCCTGTTGCTGTTCACGGGTCAGTGCGTTGAAGCGCACGGCGTAGCCAGAGACACg       >  1:398507/1‑70 (MQ=255)
                                                                    aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACTGATgg  <  1:87811/71‑1 (MQ=255)
                                                                    aGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATgg  <  1:3814349/71‑1 (MQ=255)
                                                                      taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg      <  1:5158368/65‑1 (MQ=255)
                                                                      taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACgg      <  1:3326983/65‑1 (MQ=255)
                                                                      taataaCGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACg       <  1:366473/64‑1 (MQ=255)
                                                                      |                                                                    
CCTTTATTGTACGCTTTTTACTGTACGATTTCAGTCAAATCTAATTACATAGATTGAGTGAAGGTACGAGTAATAACGTCCTGCTGCTGTTCTTTAGTCAGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGG  >  minE/607927‑608065

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