Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,254,764 T→C 13.7% P758P (CCA→CCG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,254,7640TC13.7% 184.0 / 19.3 73P758P (CCA→CCGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (23/40);  new base C (6/4);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 1.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATC  >  Exported/2254699‑2254830
                                                                 |                                                                  
cGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGt                                                               <  1:3602964/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                           <  1:1036823/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                           <  1:2490041/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                           <  1:3680419/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                           <  1:3674579/71‑1 (MQ=255)
    cACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCg                                                           <  1:866263/71‑1 (MQ=255)
      ccGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgat                                                         <  1:3583096/71‑1 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                       >  1:4459948/1‑71 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                       >  1:2753853/1‑71 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                       >  1:1140902/1‑71 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                       <  1:2328737/71‑1 (MQ=255)
        gtgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatga                                                       >  1:2141417/1‑71 (MQ=255)
         tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                                      >  1:4466266/1‑71 (MQ=255)
         tgggtggAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCgatgat                                                      >  1:3934677/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:3370327/71‑1 (MQ=255)
                 gaTGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCAt                                              <  1:1619030/71‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                             <  1:2083262/52‑1 (MQ=255)
                     aaGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGttt                                                             <  1:5078207/52‑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:740292/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:2802122/1‑70 (MQ=255)
                            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                    >  1:4010258/1‑70 (MQ=255)
                            ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                    >  1:3499390/1‑70 (MQ=255)
                                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                          <  1:3306006/60‑1 (MQ=255)
                                tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                          <  1:4690708/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:2760828/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                        >  1:1737997/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                        >  1:1581500/1‑60 (MQ=255)
                                  cGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                        >  1:294702/1‑60 (MQ=255)
                                  agTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCGGCt                                        >  1:461485/2‑60 (MQ=255)
                                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAcc                           <  1:1159017/71‑1 (MQ=255)
                                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc                            >  1:4213329/1‑70 (MQ=255)
                                    ttCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGAc                            >  1:1283729/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:2769271/71‑1 (MQ=255)
                                                cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc               <  1:4027435/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)
                                                                 |                                                                  
CGAGCACCGTGGGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATC  >  Exported/2254699‑2254830

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: