Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,254,770 T→C 15.1% E756E (GAA→GAG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,254,7700TC15.1% 158.9 / 19.7 66E756E (GAA→GAGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (23/33);  new base C (6/4);  total (29/37)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

CACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGC  >  Exported/2254703‑2254835
                                                                   |                                                                 
cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:1036823/71‑1 (MQ=255)
cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:3674579/71‑1 (MQ=255)
cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:3680419/71‑1 (MQ=255)
cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:2490041/71‑1 (MQ=255)
cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                                <  1:866263/71‑1 (MQ=255)
  ccGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgat                                                              <  1:3583096/71‑1 (MQ=255)
    gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:2753853/1‑71 (MQ=255)
    gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:1140902/1‑71 (MQ=255)
    gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:4459948/1‑71 (MQ=255)
    gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            <  1:2328737/71‑1 (MQ=255)
    gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                            >  1:2141417/1‑71 (MQ=255)
     tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                                           >  1:3934677/1‑71 (MQ=255)
     tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                                           >  1:4466266/1‑71 (MQ=255)
       ggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGatgatg                                                          <  1:250970/70‑1 (MQ=255)
           gagaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGc                                                      >  1:1973328/1‑70 (MQ=255)
            agaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCa                                                    <  1:3638102/71‑1 (MQ=255)
             gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                   <  1:1619030/71‑1 (MQ=255)
             gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                                   <  1:3370327/71‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:2658214/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:4158060/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:4338288/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:593299/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:5078207/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:740292/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                                  <  1:2083262/52‑1 (MQ=255)
                 aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCgg                                               >  1:3747987/1‑71 (MQ=255)
                  aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                                              <  1:2367583/71‑1 (MQ=255)
                  aGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGc                                              <  1:3297121/71‑1 (MQ=255)
                        ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCg                                          >  1:3662170/1‑69 (MQ=255)
                        ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:3499390/1‑70 (MQ=255)
                        ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:2802122/1‑70 (MQ=255)
                        ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:4010258/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:4047195/1‑60 (MQ=255)
                              cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:294702/1‑60 (MQ=255)
                              cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:2760828/1‑60 (MQ=255)
                              cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:1581500/1‑60 (MQ=255)
                              cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                             >  1:1737997/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:3303827/70‑1 (MQ=255)
                                   aGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCAc                              <  1:4016088/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:1202514/71‑1 (MQ=255)
                                            cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:4048458/71‑1 (MQ=255)
                                            cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:4027435/71‑1 (MQ=255)
                                            cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:7006/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:4762503/71‑1 (MQ=255)
                                                        gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        >  1:4995036/1‑71 (MQ=255)
                                                        gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:5072192/71‑1 (MQ=255)
                                                        gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:4459727/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:2461947/54‑1 (MQ=255)
                                                                  tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt               <  1:1519824/54‑1 (MQ=255)
                                                                  tttCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCt               <  1:4115351/54‑1 (MQ=255)
                                                                   |                                                                 
CACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGC  >  Exported/2254703‑2254835

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: