Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 275,782 G→A 13.9% H131H (CAC→CAT cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE275,7820GA13.9% 184.4 / 18.6 72H131H (CAC→CATcyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base A (4/6);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAG  >  minE/275715‑275845
                                                                   |                                                               
ccGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTg                                                               <  1:4296316/70‑1 (MQ=255)
            cGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTAc                                                  <  1:3142322/71‑1 (MQ=255)
               tCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGttt                                                         <  1:3986820/61‑1 (MQ=255)
               tCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGttt                                                         <  1:2782101/61‑1 (MQ=255)
                 aCCGGCGGCGCAGCGCCGGAGAGCGCTCCAATCAACGTGCCGTAGACAGAGTGGCGtt                                                          >  1:4628987/1‑58 (MQ=255)
                 aCCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGtttg                                                        >  1:4588113/1‑59 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGtt                                                          >  1:1979808/1‑57 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                   >  1:1755552/1‑64 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                            <  1:647212/71‑1 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                            <  1:3920803/71‑1 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                            >  1:3985082/1‑71 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                            >  1:1748151/1‑71 (MQ=255)
                  ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                            <  1:587487/71‑1 (MQ=255)
                    ggcggcGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                          <  1:4844492/71‑1 (MQ=255)
                    ggcggcGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                          <  1:231100/71‑1 (MQ=255)
                    ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                          <  1:2368056/68‑1 (MQ=255)
                    ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                          <  1:3221242/68‑1 (MQ=255)
                    ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                          <  1:3986980/68‑1 (MQ=255)
                    ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                          <  1:3544004/68‑1 (MQ=255)
                     gcggcgCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                          >  1:4447083/1‑70 (MQ=255)
                     gaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                          <  1:2728495/68‑1 (MQ=255)
                     gaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                          <  1:1735597/68‑1 (MQ=255)
                        gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACATAGTGGCGTTTCATGTa                                                   <  1:2284072/58‑1 (MQ=255)
                        gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                   <  1:4120523/58‑1 (MQ=255)
                        gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                   <  1:4242085/58‑1 (MQ=255)
                        gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                   <  1:4422122/58‑1 (MQ=255)
                        gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                   <  1:919926/58‑1 (MQ=255)
                        gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGATTCATGTa                                                   <  1:4527918/58‑1 (MQ=255)
                         cgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGc                                     <  1:332874/71‑1 (MQ=255)
                                    agagCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTAGAGGCTATAAACGCCGCCATAAAcc                          >  1:4686549/1‑71 (MQ=255)
                                     gagagaACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAc                           >  1:2182092/1‑69 (MQ=255)
                                      agCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTc                                                        <  1:114723/39‑1 (MQ=255)
                                      agCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTc                                                        <  1:5159024/39‑1 (MQ=255)
                                      agCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTc                                                        <  1:5122471/39‑1 (MQ=255)
                                       gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGtt                                                          >  1:2352134/1‑36 (MQ=255)
                                       gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGttc                                                  >  1:3957869/1‑42 (MQ=255)
                                       gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATaaa                                        >  1:1751736/1‑54 (MQ=255)
                                       gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATaaa                                        >  1:1627046/1‑54 (MQ=255)
                                       gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTAAAGGCTATaaa                                        >  1:2226716/1‑54 (MQ=255)
                                        cGAACCAATCAACGTGCCGTAGACAGAGTGGCGtt                                                          >  1:4293477/1‑35 (MQ=255)
                                        cGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCAt                                                      >  1:1919448/1‑39 (MQ=255)
                                        cGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCAt                                                      >  1:1357461/1‑39 (MQ=255)
                                        cGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCAt                                                      >  1:265585/1‑39 (MQ=255)
                                         gAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                              <  1:4113407/62‑1 (MQ=255)
                                         gAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                              <  1:20914/62‑1 (MQ=255)
                                         gAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                              <  1:1359386/62‑1 (MQ=255)
                                          aaCCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAg                    <  1:3591443/71‑1 (MQ=255)
                                              aaTCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCa                <  1:3948864/71‑1 (MQ=255)
                                               aTCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCTCAAAGCCCAt               <  1:4237202/71‑1 (MQ=255)
                                               aTCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAt               <  1:830229/71‑1 (MQ=255)
                                               aTCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAt               <  1:1974583/71‑1 (MQ=255)
                                                 cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATc              >  1:2980946/1‑70 (MQ=255)
                                                 cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCa             >  1:897753/1‑71 (MQ=255)
                                                 cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCa             >  1:2975155/1‑71 (MQ=255)
                                                 cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCa             >  1:2259594/1‑71 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAgg                                               >  1:2543747/1‑36 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGc                                              >  1:2405654/1‑37 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGc                                              >  1:5026328/1‑37 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGc                                              >  1:889031/1‑37 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAtggg            >  1:4314018/1‑68 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAtggg            >  1:908327/1‑68 (MQ=255)
                                                  aaCGTGCCGTAGACAGAGTGGCGATTCATGTACAg                                                >  1:5167827/1‑35 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:1319475/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:4714381/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:1208184/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:4302455/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:523196/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:120653/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:2232473/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                          <  1:3946307/56‑1 (MQ=255)
                                                   aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCAcc           >  1:1147603/1‑71 (MQ=255)
                                                             gACAGAGTGGCGTTTCATGTAAAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCccagccag  >  1:4073825/1‑70 (MQ=255)
                                                                   |                                                               
CCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAG  >  minE/275715‑275845

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: