Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 222,997 T→A 35.9% intergenic (+52/‑128) lpcA → / → prfH D‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE222,9970TA35.9% 73.9 / 46.7 50intergenic (+52/‑128)lpcA/prfHD‑sedoheptulose 7‑phosphate isomerase/predicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base T (18/11);  new base A (6/11);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

AAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAA  >  minE/222927‑223066
                                                                      |                                                                     
aaaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                       <  1:1474638/71‑1 (MQ=255)
 aaaGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGct                                                                       >  1:3409786/1‑70 (MQ=255)
       aTGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctg                                                                      <  1:1277879/65‑1 (MQ=255)
        tGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcc                                                                 >  1:836347/1‑69 (MQ=255)
         ggTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgat                                                                    >  1:38340/1‑65 (MQ=255)
           ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctga                                                              <  1:2551613/69‑1 (MQ=255)
           ttAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatg                                                                   >  1:2776481/1‑64 (MQ=255)
                  agtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                                                     >  1:3113530/1‑71 (MQ=255)
                   gtCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTg                                                    <  1:1454332/71‑1 (MQ=255)
                       tccGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAAcgc                                                >  1:1188511/3‑71 (MQ=14)
                        gCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGAtgcg                                               >  1:67823/1‑71 (MQ=255)
                          gTAGGCATGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgatg                                                            >  1:1449724/1‑56 (MQ=21)
                            aGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCt                                                     >  1:3040764/1‑61 (MQ=255)
                              gCCGGACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATCa                                       <  1:3162524/71‑1 (MQ=9)
                               ccGGATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATc                                        <  1:1109263/71‑1 (MQ=14)
                                 ggACAAGACGTTTA‑‑TCGCATCCGGCAGTCGTGTGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGc                                    <  1:2257632/71‑1 (MQ=2)
                                  gATAAGGCGTTTACGCCGCATCCGGCAGTCGTGCGcagatgcctgatgc                                                           <  1:3630928/49‑1 (MQ=11)
                                   aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGctgatgcctgat                                                             <  1:3661025/46‑1 (MQ=255)
                                   aTAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGc                                    >  1:368178/1‑71 (MQ=255)
                                      aggcGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTa                                 >  1:3539000/4‑71 (MQ=11)
                                      aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                  <  1:1566096/67‑1 (MQ=11)
                                      aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                  <  1:1592322/67‑1 (MQ=11)
                                      aggcGTTTACGCCGCATCCGGCAGTTGTACGCAGGTGCCTGATGCGACGCTGGCGCGTCTTATCATGCCt                                  <  1:330049/67‑1 (MQ=11)
                                      aggcGTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTa                                 >  1:2336579/4‑71 (MQ=11)
                                       gACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTa                                 >  1:1374841/1‑70 (MQ=255)
                                        aCGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCTTGCCTACa                               <  1:1527663/71‑1 (MQ=255)
                                         cGTTTATGCCGCATCCGGCAGTTATGCGCAGATGCCTGATGCGACGCTGGCGCGTCTTATCATGCCTACaa                              <  1:2412485/71‑1 (MQ=14)
                                         cGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACa                               >  1:1412261/1‑70 (MQ=255)
                                          gTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACaaa                             >  1:770164/1‑71 (MQ=255)
                                          gTTTACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCag                                      <  1:203916/62‑2 (MQ=14)
                                           tttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt                            <  1:2860960/71‑1 (MQ=255)
                                            ttACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATc                           >  1:2275068/1‑71 (MQ=255)
                                            tcaCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                          >  1:1025131/3‑56 (MQ=14)
                                              aCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAttcc                         >  1:2556266/1‑68 (MQ=11)
                                              aCGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGACCTACAAAttct                         <  1:3662595/71‑4 (MQ=11)
                                                gccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTa                       >  1:323003/1‑71 (MQ=255)
                                                 ccgcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAc                      <  1:3931188/71‑1 (MQ=255)
                                                  tgcATCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCTa                                 <  1:442082/57‑1 (MQ=14)
                                                   gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc                    >  1:2146681/1‑71 (MQ=255)
                                                   gcATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTAcgc                    >  1:1371132/1‑71 (MQ=255)
                                                     aTCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGa                  <  1:483843/71‑1 (MQ=255)
                                                      tCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACAAAt                            >  1:747356/1‑60 (MQ=11)
                                                      tCCGGCACTGTGTGCCGATGCCTGATGCGACGCTGACGCGTTTTAt                                         <  1:430686/46‑1 (MQ=2)
                                                       ccGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAAt                            >  1:457429/1‑59 (MQ=255)
                                                          gCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCAtcc                                    >  1:2705439/1‑46 (MQ=255)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg              >  1:2067074/1‑69 (MQ=255)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCg              >  1:1868962/1‑69 (MQ=255)
                                                           cATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTa            <  1:2612615/71‑1 (MQ=255)
                                                            aTTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGt             <  1:381638/69‑1 (MQ=255)
                                                            aTTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAg           >  1:1276661/1‑71 (MQ=255)
                                                             actgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                  <  1:3002515/45‑1 (MQ=14)
                                                             actgtgtgCCGATGCCTGATGCGACGCTGACGCGTTTTATCATGCCt                                  <  1:225542/45‑1 (MQ=14)
                                                              ttGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGc         >  1:768215/1‑71 (MQ=255)
                                                              tcgtgCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACa                               <  1:905556/47‑1 (MQ=11)
                                                               cgtgCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTACa                               <  1:430042/47‑1 (MQ=11)
                                                                 tgCGCTGATGCCTGATGCGACGCTGACGCGTCTTa                                          <  1:1956420/35‑1 (MQ=255)
                                                                 tgCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCg       <  1:2390605/70‑1 (MQ=255)
                                                                  gcgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaa     >  1:2012727/1‑71 (MQ=255)
                                                                   cgcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaat    <  1:3570070/71‑1 (MQ=255)
                                                                    gcTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaata   <  1:2282270/71‑1 (MQ=255)
                                                                     ctgatgcctgatgcGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGaataa  >  1:3428121/1‑71 (MQ=255)
                                                                      |                                                                     
AAAAGAGATGGTTAAGTAAGTCTGGCGTAGGCCGGATAAGACGTTTACGCCGCATCCGGCATTTGTGCGCTGATGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACAAATCTGTACGCGAACCGTAGGCCGAATAA  >  minE/222927‑223066

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: