Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,253,051 A→G 41.7% T1329T (ACT→ACC rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,253,0510AG41.7% 29.4 / 69.4 60T1329T (ACT→ACCrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (13/22);  new base G (17/8);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 3.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

GTTCTCTTTCAGGCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATAC  >  Exported/2252988‑2253115
                                                               |                                                                
gTTCTCTTTCAGACCACGCAGGAAGTCGCGCTTGCCGGTGACCGCCGCCTCGGTGAGGACGCGGGTGGtct                                                           >  1:3843741/1‑71 (MQ=255)
gTTCTCTTTCAGACCACGCAGGAAGTCGCGCTTGCCGGTGACCGCCGCCTCGGTGAGGACGCGGGTGGtct                                                           >  1:1357261/1‑71 (MQ=255)
gTTCTCTTTCAGACCACGCAGGAAGTCGCGCTTGCCGGTGACCGCCGCCTCGGTGAGGACGCGGGTGGtct                                                           >  1:3086410/1‑71 (MQ=255)
gTTCTCTTTCAGACCACGCAGGAAGTCGCGCTTGCCGGTGACCGCCGCCTCGGTGAGGACGCGGGTGGtct                                                           >  1:2904359/1‑71 (MQ=255)
     ctTTCAGGCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTg                                                       >  1:118890/1‑70 (MQ=255)
     ctTTCAGGCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTg                                                       >  1:4902793/1‑70 (MQ=255)
            gCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATg                                               <  1:4869075/71‑1 (MQ=255)
            gCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATg                                               <  1:3488794/71‑1 (MQ=255)
            gCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATg                                               <  1:3113567/71‑1 (MQ=255)
            gCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATg                                               <  1:1438947/71‑1 (MQ=255)
            gCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATg                                               <  1:153693/71‑1 (MQ=255)
            gCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATg                                               <  1:793870/71‑1 (MQ=255)
             ccGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGc                                              <  1:99208/71‑1 (MQ=255)
                cgcAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCcgc                                           <  1:4151091/71‑1 (MQ=255)
                cgcAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCcgc                                           <  1:3930343/71‑1 (MQ=255)
                      tcgtcgCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGAt                                     >  1:446644/1‑71 (MQ=255)
                        gtcgCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGa                                   <  1:4677691/71‑1 (MQ=255)
                        gtcgCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGa                                   <  1:3792987/71‑1 (MQ=255)
                                tggccGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTGTCCTGGAACGATGCCGCGGAGATGAAGGACTc                            <  1:4065730/67‑1 (MQ=255)
                                 gCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCg                           <  1:1987193/70‑1 (MQ=255)
                                   ccGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGtt                        >  1:4323849/1‑71 (MQ=255)
                                   ccGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGtt                        >  1:1041079/1‑71 (MQ=255)
                                    cGCAACGGCTGCTTCGGTCAGCACGCGAGTGGGCTCCTGGAACGATg                                               <  1:3660481/47‑1 (MQ=255)
                                                 tCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCttt          >  1:3491352/1‑71 (MQ=255)
                                                  cGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAggt                               >  1:3801863/1‑48 (MQ=255)
                                                  cGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCtt           <  1:1634887/69‑1 (MQ=255)
                                                   ggTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTgg        <  1:1251468/71‑1 (MQ=255)
                                                    gTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGt       >  1:2425062/1‑71 (MQ=255)
                                                       aGCACGGGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:2737307/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:4396902/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:1224086/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:1256645/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:1253965/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:2295718/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:2681639/42‑2 (MQ=255)
                                                       aGCACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGaac                                 <  1:4302292/42‑2 (MQ=255)
                                                       aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGa     >  1:2187203/1‑70 (MQ=255)
                                                       aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGa     >  1:4706076/1‑70 (MQ=255)
                                                       aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGa     >  1:3678840/1‑70 (MQ=255)
                                                       aGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAACGACTCGGTTGCCa                    >  1:2363949/1‑55 (MQ=255)
                                                        gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa   <  1:638836/71‑1 (MQ=255)
                                                        gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa   <  1:3874895/71‑1 (MQ=255)
                                                        gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa   <  1:845628/71‑1 (MQ=255)
                                                        gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa   <  1:1362392/71‑1 (MQ=255)
                                                        gCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa   <  1:3052781/71‑1 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:838693/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:357429/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:4950113/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:58809/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:860139/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:3801056/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:3768434/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:2832008/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:2244842/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:1985820/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:1504887/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGATGCGGCCGAGATGAACGa                               >  1:1059048/1‑42 (MQ=255)
                                                         cACGCGGGTGGTTTCCTGGAACGACGCGGCCGAGATGAACGa                               >  1:4894285/1‑42 (MQ=16)
                                                         cACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATa   <  1:2065545/70‑1 (MQ=255)
                                                         cACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATAc  >  1:4366853/1‑71 (MQ=255)
                                                               |                                                                
GTTCTCTTTCAGGCCGCGCAGTTCGTCGCGTTTGCCCGCAACGGCTGCTTCGGTCAGCACGCGAGTGGTCTCCTGGAACGATGCCGCGGAGATGAAGGACTCGGTTGCCAGAGACGCTTTGGTGATAC  >  Exported/2252988‑2253115

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: