Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 609,333 G→A 16.2% D307D (GAC→GAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE609,3330GA16.2% 155.7 / 19.6 68D307D (GAC→GATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base A (3/8);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 9.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCA  >  minE/609264‑609403
                                                                     |                                                                      
gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                                                       >  1:3271539/1‑71 (MQ=255)
gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                                                                       >  1:5303577/1‑71 (MQ=255)
  cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                                     <  1:5031927/71‑1 (MQ=255)
  cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                                     <  1:278814/71‑1 (MQ=255)
    gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                                    >  1:2931471/1‑69 (MQ=255)
    gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                                    >  1:878281/1‑69 (MQ=255)
      gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                                                                 <  1:1885619/71‑1 (MQ=255)
        aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                                               >  1:1719359/1‑71 (MQ=255)
        aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                                               >  1:4089161/1‑71 (MQ=255)
         cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                               >  1:2999037/1‑70 (MQ=255)
              cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                                                         >  1:3684893/1‑69 (MQ=255)
              cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                         >  1:1981534/1‑71 (MQ=255)
              cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                         >  1:3693472/1‑71 (MQ=255)
                 caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                                                       >  1:1560493/4‑70 (MQ=255)
                  gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                                                     >  1:4107251/1‑71 (MQ=255)
                  aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                                                      <  1:939224/68‑2 (MQ=255)
                     ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTcta                                                                    >  1:958715/1‑51 (MQ=255)
                     ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                                                  >  1:5283728/1‑71 (MQ=255)
                      tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                               <  1:5318468/57‑1 (MQ=255)
                      tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                               <  1:4064226/57‑1 (MQ=255)
                       ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                                                                    >  1:1607554/1‑51 (MQ=255)
                        cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                               <  1:236672/71‑1 (MQ=255)
                        cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                               <  1:4578514/71‑1 (MQ=255)
                          gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                          >  1:1729511/3‑58 (MQ=255)
                          gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                          >  1:4909761/3‑58 (MQ=255)
                          gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                                     <  1:3155037/47‑1 (MQ=255)
                          gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                                              <  1:4530147/54‑1 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:339874/60‑3 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:3122764/60‑3 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:4567928/60‑3 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:4251922/60‑3 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:1205687/60‑3 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:2205378/60‑3 (MQ=255)
                             tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                     <  1:1301603/60‑3 (MQ=255)
                               cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                             >  1:299515/1‑50 (MQ=255)
                               cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                             >  1:1787213/1‑50 (MQ=255)
                                  aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                             >  1:2102882/1‑47 (MQ=255)
                                  aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                             >  1:4362599/1‑47 (MQ=255)
                                  aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                             >  1:1682521/1‑47 (MQ=255)
                                   ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                               <  1:242077/44‑1 (MQ=255)
                                   ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                               <  1:811417/44‑1 (MQ=255)
                                   ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                    >  1:299446/1‑71 (MQ=255)
                                   ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                    >  1:1523225/1‑71 (MQ=255)
                                   ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                    <  1:3169799/71‑1 (MQ=255)
                                    gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt                                      >  1:4128862/1‑68 (MQ=255)
                                    gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                   >  1:5601422/1‑71 (MQ=255)
                                      aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                    >  1:5283690/1‑68 (MQ=255)
                                      aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                   >  1:2070553/1‑69 (MQ=255)
                                         gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc                                       <  1:4481714/62‑1 (MQ=255)
                                         gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                              >  1:990197/1‑71 (MQ=255)
                                         gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                              <  1:1970934/71‑1 (MQ=255)
                                            ccATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                             >  1:3614193/1‑37 (MQ=255)
                                              aTACGCAGTTTCTTGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                               <  1:30087/49‑1 (MQ=255)
                                                 cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                      <  1:4365751/71‑1 (MQ=255)
                                                 cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                      <  1:297722/71‑1 (MQ=255)
                                                     gTTTCATGACCAGGTGGTCTACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGt                  <  1:5111201/71‑1 (MQ=255)
                                                      tttCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                            <  1:960013/60‑1 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                     >  1:422164/1‑49 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                     >  1:4065374/1‑49 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                <  1:4254341/70‑1 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                <  1:1026504/70‑1 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                <  1:1665705/70‑1 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                >  1:728656/1‑70 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                <  1:2058404/70‑1 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                >  1:3027650/1‑70 (MQ=255)
                                                        tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                <  1:333195/70‑1 (MQ=255)
                                                           tGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATg            >  1:178792/1‑71 (MQ=255)
                                                               cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                          >  1:2169634/1‑53 (MQ=255)
                                                               cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                         >  1:2332258/1‑54 (MQ=255)
                                                               cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                         >  1:1121064/1‑54 (MQ=255)
                                                                 ggtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                          >  1:680557/1‑51 (MQ=255)
                                                                  gtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACAt     >  1:5601343/1‑71 (MQ=255)
                                                                    ggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCa                                 >  1:2663644/1‑41 (MQ=255)
                                                                     gtCAACCATTTCCTGCTCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATcc   <  1:1769834/70‑1 (MQ=255)
                                                                     gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc    >  1:2417452/1‑69 (MQ=255)
                                                                     gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc    >  1:4725418/1‑69 (MQ=255)
                                                                     gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCa  <  1:4826755/71‑1 (MQ=255)
                                                                     |                                                                      
GCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCA  >  minE/609264‑609403

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: