Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,254,785 A→G 15.9% D751D (GAT→GAC rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,254,7850AG15.9% 150.0 / 19.9 63D751D (GAT→GACrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (23/30);  new base G (6/4);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

AGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAC  >  Exported/2254715‑2254852
                                                                      |                                                                   
agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa                                                                     <  1:3638102/71‑1 (MQ=255)
 gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                                    <  1:3370327/71‑1 (MQ=255)
 gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                                    <  1:1619030/71‑1 (MQ=255)
     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                                                >  1:3747987/1‑71 (MQ=255)
      aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                                                               <  1:3297121/71‑1 (MQ=255)
      aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                                                               <  1:2367583/71‑1 (MQ=255)
            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCg                                                           >  1:3662170/1‑69 (MQ=255)
            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                                          >  1:2802122/1‑70 (MQ=255)
            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                                          >  1:4010258/1‑70 (MQ=255)
            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                                          >  1:3499390/1‑70 (MQ=255)
                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                                                <  1:4690708/60‑1 (MQ=255)
                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                                                <  1:3306006/60‑1 (MQ=255)
                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                                                <  1:3641668/60‑1 (MQ=255)
                 aCGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                                                <  1:1160143/59‑1 (MQ=255)
                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                                              >  1:1581500/1‑60 (MQ=255)
                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                                              >  1:1737997/1‑60 (MQ=255)
                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                                              >  1:4047195/1‑60 (MQ=255)
                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                                              >  1:294702/1‑60 (MQ=255)
                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                                              >  1:2760828/1‑60 (MQ=255)
                  agTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                                              >  1:461485/2‑60 (MQ=255)
                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAcc                                                 <  1:1159017/71‑1 (MQ=255)
                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc                                                  >  1:1283729/1‑70 (MQ=255)
                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc                                                  >  1:4213329/1‑70 (MQ=255)
                      cAGAGCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                                               >  1:1040516/1‑71 (MQ=255)
                      cAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                                               >  1:3658936/1‑71 (MQ=255)
                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                                               <  1:4016088/70‑1 (MQ=255)
                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                                               <  1:3303827/70‑1 (MQ=255)
                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACg                                              <  1:3205649/71‑1 (MQ=255)
                       aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATACGGCTTCGCCATCAGACCAc                                               <  1:3090693/70‑1 (MQ=255)
                           cTTCACGGAAGTTCGCGGTGATTGGCGTTTCGAAGATGGAGCCATCCgg                                                                <  1:2967508/49‑1 (MQ=255)
                            ttCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                                                <  1:4973791/48‑1 (MQ=255)
                             tCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa                                            <  1:886499/67‑1 (MQ=255)
                              cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCa                                            <  1:4818177/66‑1 (MQ=255)
                              cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                                         >  1:3823471/1‑69 (MQ=255)
                              cACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATAc                                         >  1:1726566/1‑69 (MQ=255)
                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                                     <  1:4048458/71‑1 (MQ=255)
                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                                     <  1:1202514/71‑1 (MQ=255)
                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                                     <  1:7006/71‑1 (MQ=255)
                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                                     <  1:4027435/71‑1 (MQ=255)
                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                                     <  1:2769271/71‑1 (MQ=255)
                                 ggAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                                     >  1:1557053/1‑70 (MQ=255)
                                    aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                                                <  1:2906056/40‑1 (MQ=255)
                                    aGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGc                                 >  1:3509438/1‑71 (MQ=255)
                                          cgGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGa                           <  1:1402842/71‑1 (MQ=255)
                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt                         <  1:5072192/71‑1 (MQ=255)
                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt                         <  1:4459727/71‑1 (MQ=255)
                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt                         >  1:4995036/1‑71 (MQ=255)
                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt                         <  1:4762503/71‑1 (MQ=255)
                                            gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt                         >  1:3248259/1‑71 (MQ=255)
                                             tGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATc                        >  1:1870126/1‑71 (MQ=255)
                                                   gCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGc                   <  1:5027828/70‑1 (MQ=255)
                                                      tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt                                <  1:1519824/54‑1 (MQ=255)
                                                      tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt                                <  1:4115351/54‑1 (MQ=255)
                                                      tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt                                <  1:2461947/54‑1 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:2925638/1‑61 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:2466236/1‑61 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:5184857/1‑61 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:1859121/1‑61 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:4603414/1‑61 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:433424/1‑61 (MQ=255)
                                                        tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt                       >  1:1575597/1‑61 (MQ=255)
                                                                   gCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAc  <  1:621537/71‑1 (MQ=255)
                                                                   gCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAc  <  1:3800115/71‑1 (MQ=255)
                                                                    ccATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAc  <  1:1569903/70‑1 (MQ=255)
                                                                      |                                                                   
AGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGCCGCAGAACCACGCGCAC  >  Exported/2254715‑2254852

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: