Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 275,785 G→A 13.5% R130R (CGC→CGT cyoE ← protoheme IX farnesyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported275,7850GA13.5% 192.9 / 17.9 74R130R (CGC→CGTcyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (27/37);  new base A (4/6);  total (31/43)
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.83e-01

CGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGC  >  Exported/275727‑275854
                                                          |                                                                     
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:1748151/1‑71 (MQ=255)
      ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                                     <  1:3920803/71‑1 (MQ=255)
      ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                                     >  1:3985082/1‑71 (MQ=255)
      ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                                     <  1:587487/71‑1 (MQ=255)
      ccGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCta                                                     <  1:647212/71‑1 (MQ=255)
        ggcggcGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                                   <  1:231100/71‑1 (MQ=255)
        ggcggcGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                                   <  1:4844492/71‑1 (MQ=255)
        ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                                   <  1:3986980/68‑1 (MQ=255)
        ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                                   <  1:3221242/68‑1 (MQ=255)
        ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                                   <  1:3544004/68‑1 (MQ=255)
        ggaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                                   <  1:2368056/68‑1 (MQ=255)
         gcggcgCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                                   >  1:4447083/1‑70 (MQ=255)
         gaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                                   <  1:1735597/68‑1 (MQ=255)
         gaggcgCAGCGCCGGAGAGCGAACCAATCAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCtata                                                   <  1:2728495/68‑1 (MQ=255)
            gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACATAGTGGCGTTTCATGTa                                                            <  1:2284072/58‑1 (MQ=255)
            gcgcAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                            <  1:919926/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)
            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:1627046/1‑54 (MQ=255)
                           gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATaaa                                                 >  1:1751736/1‑54 (MQ=255)
                           gCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTAAAGGCTATaaa                                                 >  1:2226716/1‑54 (MQ=255)
                            cGAACCAATCAACGTGCCGTAGACAGAGTGGCGtt                                                                   >  1:4293477/1‑35 (MQ=255)
                            cGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCAt                                                               >  1:1357461/1‑39 (MQ=255)
                            cGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCAt                                                               >  1:265585/1‑39 (MQ=255)
                            cGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCAt                                                               >  1:1919448/1‑39 (MQ=255)
                             gAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                                       <  1:1359386/62‑1 (MQ=255)
                             gAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                                       <  1:20914/62‑1 (MQ=255)
                             gAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                                       <  1:4113407/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:1974583/71‑1 (MQ=255)
                                   aTCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAt                        <  1:830229/71‑1 (MQ=255)
                                     cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATc                       >  1:2980946/1‑70 (MQ=255)
                                     cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCa                      >  1:897753/1‑71 (MQ=255)
                                     cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCa                      >  1:2259594/1‑71 (MQ=255)
                                     cAAGGTGCCGTAGACAGAATGACGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCa                      >  1:2975155/1‑71 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAgg                                                        >  1:2543747/1‑36 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGc                                                       >  1:889031/1‑37 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGc                                                       >  1:2405654/1‑37 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGc                                                       >  1:5026328/1‑37 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAtggg                     >  1:908327/1‑68 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCAtggg                     >  1:4314018/1‑68 (MQ=255)
                                      aaCGTGCCGTAGACAGAGTGGCGATTCATGTACAg                                                         >  1:5167827/1‑35 (MQ=255)
                                       aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                                   <  1:4714381/56‑1 (MQ=255)
                                       aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                                   <  1:120653/56‑1 (MQ=255)
                                       aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                                   <  1:523196/56‑1 (MQ=255)
                                       aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                                   <  1:1208184/56‑1 (MQ=255)
                                       aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                                   <  1:1319475/56‑1 (MQ=255)
                                       aCGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAcc                                   <  1:4302455/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)
                                                        tGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCAcccc                  <  1:2754027/56‑1 (MQ=255)
                                                        tGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCAcccc                  <  1:4079880/56‑1 (MQ=255)
                                                        tGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCAcccc                  <  1:3020647/56‑1 (MQ=255)
                                                         ggCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGc  >  1:3216627/1‑71 (MQ=255)
                                                          |                                                                     
CGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGC  >  Exported/275727‑275854

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: