Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,254,779 G→C 16.1% S753S (TCC→TCG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,254,7790GC16.1% 151.3 / 21.2 62S753S (TCC→TCGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (24/28);  new base C (6/4);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGC  >  Exported/2254710‑2254835
                                                                     |                                                        
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)
          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:4010258/1‑70 (MQ=255)
                 ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:2802122/1‑70 (MQ=255)
                 ggAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGc                                         >  1:3499390/1‑70 (MQ=255)
                     tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:3306006/60‑1 (MQ=255)
                     tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:3641668/60‑1 (MQ=255)
                     tACGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:4690708/60‑1 (MQ=255)
                      aCGTTCAGGCCTTCCCGGAAGTTCGCGGTGATCGGGGTCTCGATGATCGAGCCGTCCgg                                               <  1:1160143/59‑1 (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:1581500/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: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:2769271/71‑1 (MQ=255)
                                     cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:7006/71‑1 (MQ=255)
                                     cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:4027435/71‑1 (MQ=255)
                                     cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:4048458/71‑1 (MQ=255)
                                     cGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGc                    <  1:1202514/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:3248259/1‑71 (MQ=255)
                                                 gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:5072192/71‑1 (MQ=255)
                                                 gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        >  1:4995036/1‑71 (MQ=255)
                                                 gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:4762503/71‑1 (MQ=255)
                                                 gTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAAt        <  1:4459727/71‑1 (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:4603414/1‑61 (MQ=255)
                                                             tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt      >  1:433424/1‑61 (MQ=255)
                                                             tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt      >  1:1575597/1‑61 (MQ=255)
                                                             tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt      >  1:2925638/1‑61 (MQ=255)
                                                             tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt      >  1:1859121/1‑61 (MQ=255)
                                                             tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt      >  1:5184857/1‑61 (MQ=255)
                                                             tCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCt      >  1:2466236/1‑61 (MQ=255)
                                                                     |                                                        
GGTGGAGATGAAGTACTGGAGTACGTTCAGACCTTCACGGAAGTTCGCGGTGATTGGCGTTTCGATGATGGAGCCATCCGGCTTCGCCATCAGACCACGCATACCAGCAAGCTGACGAATCTGTGC  >  Exported/2254710‑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: