Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 394,838 G→A 15.7% P63P (CCG→CCA ahpF → alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE394,8380GA15.7% 121.7 / 13.9 51P63P (CCG→CCAahpFalkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base G (22/21);  new base A (6/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.69e-01

TGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGG  >  minE/394768‑394900
                                                                      |                                                              
tGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCg                                                                >  1:1171283/1‑71 (MQ=255)
 ggCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCg                                                                <  1:735979/70‑1 (MQ=255)
   cTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCt                                                             >  1:2543070/1‑71 (MQ=255)
   cTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCt                                                             >  1:1098656/1‑71 (MQ=255)
    tGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCtt                                                            >  1:3189008/1‑71 (MQ=255)
      aaaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCtt                                                            <  1:1180880/69‑1 (MQ=255)
       aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTcc                                                         <  1:931339/71‑1 (MQ=255)
       aaTCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTcc                                                         <  1:2142631/71‑1 (MQ=255)
             aGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCa                                                   <  1:2176445/71‑1 (MQ=255)
               aaCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt                                                     <  1:353972/67‑1 (MQ=255)
               aaCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCAc                                                  >  1:1416119/1‑70 (MQ=255)
                aCTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCa                                                <  1:3943057/71‑1 (MQ=255)
                 cTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt                                                     >  1:2009776/1‑65 (MQ=255)
                 cTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGAt                                                     >  1:2130067/1‑65 (MQ=255)
                       gACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGt                                                               >  1:10810/1‑49 (MQ=255)
                        aCAAATTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGt                                        >  1:1651079/1‑71 (MQ=255)
                        aCAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGt                                        >  1:2574983/1‑71 (MQ=255)
                         cAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGtt                                       >  1:1352170/1‑71 (MQ=255)
                          aaaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAAccc                                            <  1:299581/65‑1 (MQ=255)
                          aaaGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTc                                      <  1:3524168/71‑1 (MQ=255)
                            aGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCt                                                        >  1:3827875/1‑51 (MQ=255)
                                 ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa                                <  1:3951080/70‑1 (MQ=255)
                                 ccTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               <  1:507782/71‑1 (MQ=255)
                                  cTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTg                                                       >  1:1908924/1‑46 (MQ=255)
                                  cTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               <  1:1255016/70‑1 (MQ=255)
                                    ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCaa                                               >  1:2816917/1‑52 (MQ=255)
                                    ttAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAgggg                            >  1:2664229/1‑71 (MQ=255)
                                     tAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGc                           >  1:201199/1‑71 (MQ=255)
                                            gaTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAAc                                  <  1:686568/57‑1 (MQ=255)
                                             aTAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGttt                    <  1:2591426/70‑1 (MQ=255)
                                              tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGTGCCACGTTTTg                  >  1:2968532/1‑71 (MQ=255)
                                              tAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  <  1:3230035/71‑1 (MQ=255)
                                               aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCAcc                                                 <  1:1479219/39‑1 (MQ=255)
                                               aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACg                       <  1:2657327/65‑1 (MQ=255)
                                               aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  >  1:2758197/1‑70 (MQ=255)
                                               aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  >  1:1291326/1‑70 (MQ=255)
                                               aaCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  >  1:3490388/1‑70 (MQ=255)
                                                aCAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGc                 <  1:1776191/70‑1 (MQ=255)
                                                aCAGCTTGCCGGTGCGTAAGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGc                 <  1:2854796/70‑1 (MQ=255)
                                                 cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCa                                >  1:3297002/1‑54 (MQ=255)
                                                 cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAgg                              >  1:1136182/1‑56 (MQ=255)
                                                 cAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGCGCCACGTTTTGCa                >  1:925430/1‑70 (MQ=255)
                                                    cTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCt            >  1:3939194/1‑71 (MQ=255)
                                                       gCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAgg                              <  1:2464742/50‑1 (MQ=255)
                                                           gTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCt     <  1:472128/71‑1 (MQ=255)
                                                             gCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg   >  1:2367742/1‑71 (MQ=255)
                                                              cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTggg  >  1:3784870/1‑71 (MQ=255)
                                                              cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg   >  1:3280324/1‑70 (MQ=255)
                                                              cGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTgg   <  1:734275/70‑1 (MQ=255)
                                                                tAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAg                               <  1:697844/40‑1 (MQ=255)
                                                                  aGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTg    >  1:2270697/1‑65 (MQ=255)
                                                                  aGCCATCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTg    >  1:319668/1‑65 (MQ=255)
                                                                     cGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTcccc        <  1:2045783/58‑1 (MQ=255)
                                                                      gTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTg                  <  1:181817/47‑1 (MQ=255)
                                                                      |                                                              
TGGCTGAAATCGCAGAACTGTCAGACAAAGTCACCTTTAAAGAAGATAACAGCTTGCCGGTGCGTAAGCCGTCTTTCCTGATCACCAACCCAGGTTCCAACCAGGGGCCACGTTTTGCAGGCTCCCCGCTGGG  >  minE/394768‑394900

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: