Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,255,098 G→A 7.7% P647L (CCG→CTG)  rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,255,0980GA 7.7% 205.3 / 6.4 65P647L (CCG→CTG) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (25/35);  new base A (3/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACATG  >  Exported/2255029‑2255166
                                                                     |                                                                    
gacTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCgg                                                                     >  1:3811/1‑71 (MQ=255)
   tGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGAt                                                                  >  1:2154696/1‑71 (MQ=255)
   tGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGAt                                                                  >  1:752846/1‑71 (MQ=255)
   tGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGAt                                                                  >  1:1407357/1‑71 (MQ=255)
      aaCTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCgg                                                                     <  1:1936310/65‑1 (MQ=255)
       aCTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGAcc                                                              <  1:226418/71‑1 (MQ=255)
       aCTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGAc                                                               >  1:2396318/1‑70 (MQ=255)
           cTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGt                                                          <  1:1643979/71‑1 (MQ=255)
              cTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTc                                                         <  1:1244229/69‑1 (MQ=255)
              cTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTGTCTTCTCCGGGATGACCATGTc                                                         <  1:2200537/69‑1 (MQ=255)
               tGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCgg                                                                     <  1:2720776/56‑1 (MQ=255)
                  aTTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGAt                                                   <  1:542872/71‑1 (MQ=255)
                      cAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATATCATCGATAcc                                                <  1:1309580/70‑1 (MQ=255)
                       aGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCa                                               <  1:2413213/70‑1 (MQ=255)
                        gCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAAc                                             <  1:1720516/71‑1 (MQ=255)
                          aaCTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:1586035/62‑1 (MQ=255)
                          aaCTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:3392659/62‑1 (MQ=255)
                          aaCTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:1768693/62‑1 (MQ=255)
                          aaCTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:1917864/62‑1 (MQ=255)
                          aaCTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCaa                                              >  1:1225643/1‑68 (MQ=255)
                               ctgcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGc                                       <  1:2591086/70‑1 (MQ=255)
                               ctgcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCa                                      >  1:689517/1‑71 (MQ=255)
                               ctgcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCa                                      >  1:2206753/1‑71 (MQ=255)
                               ctgcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCa                                      >  1:2612955/1‑71 (MQ=255)
                               ctgcttctgcCTCGGAGATGATTTCGTCTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCa                                      >  1:1340997/1‑71 (MQ=255)
                                tgcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGa                                                                   >  1:2916863/1‑41 (MQ=255)
                                 gcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCAc                                     >  1:1024110/1‑70 (MQ=255)
                                 gcttctgcCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCAc                                     >  1:1210356/1‑70 (MQ=255)
                                        ccTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAAcgc                             <  1:300712/71‑1 (MQ=255)
                                            ggAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGa                                          >  1:2542537/1‑54 (MQ=255)
                                            ggAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGa                                          >  1:1180507/1‑54 (MQ=255)
                                            ggAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTg                         <  1:3339319/71‑1 (MQ=255)
                                             gagaTGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATc                                                      <  1:211057/41‑1 (MQ=255)
                                             gagaTGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATc                                                      <  1:984707/41‑1 (MQ=255)
                                               gatgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCa                       >  1:2880833/1‑70 (MQ=255)
                                               gatgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCa                       >  1:3347902/1‑70 (MQ=255)
                                               gatgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCAt                      <  1:3239876/71‑1 (MQ=255)
                                               gatgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCAt                      <  1:590662/71‑1 (MQ=255)
                                                atgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCaa                                              <  1:3326271/46‑1 (MQ=255)
                                                atgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCaa                                              <  1:994329/46‑1 (MQ=255)
                                                atgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCaa                                              <  1:1636332/46‑1 (MQ=255)
                                                atgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCt                          >  1:400302/1‑66 (MQ=255)
                                                atgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCAt                      <  1:1152975/70‑1 (MQ=255)
                                                atgatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATa                     <  1:2215439/71‑1 (MQ=255)
                                                 tgatTTCGTGTTTCTTCTCCAGGATGACCATGTCATCGATACCAACAGATGCACCAGa                                 <  1:920031/58‑1 (MQ=255)
                                                 tgatTTCGTGTTTCTTCTCCAGGATGACCATGTCATCGATACCAACAGATGCACCAGa                                 <  1:575001/58‑1 (MQ=255)
                                                  gatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATAc                                                 >  1:2616099/1‑41 (MQ=255)
                                                  gatTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAgg                   >  1:2903479/1‑71 (MQ=255)
                                                   atTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGc                  >  1:1055658/1‑71 (MQ=255)
                                                   atTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGc                  >  1:2913860/1‑71 (MQ=255)
                                                   atTTCGTGTTTCTTCTCCAGGATGACCATGTCATCGATACCAACAGATGCACCAGAAc                               >  1:275265/1‑58 (MQ=255)
                                                   atTTCGTGTTTCTTCTCCAGGATGACCATGTCATCGATACCAACAGATGCACCAGAAc                               >  1:2005206/1‑58 (MQ=255)
                                                   atTTCGTGTTTCTTCTCCAGGATGACCATGTCATCGATACCAACAGATGCACCAGAAc                               >  1:3310040/1‑58 (MQ=255)
                                                    tttCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCg                 <  1:1226607/71‑1 (MQ=255)
                                                     ttCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:1726775/35‑1 (MQ=255)
                                                     ttCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:1491470/35‑1 (MQ=255)
                                                     ttCGTGTTTCTTCTCCGGGATGACCATGTCATCGa                                                    <  1:2676414/35‑1 (MQ=255)
                                                      tCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGAt                                         >  1:2623742/1‑45 (MQ=255)
                                                      tCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGa                <  1:2212466/70‑1 (MQ=255)
                                                       cGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGaa               <  1:1765327/70‑1 (MQ=255)
                                                           tttCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCg           <  1:179049/70‑1 (MQ=255)
                                                           tttCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCg           <  1:308584/70‑1 (MQ=255)
                                                                 ctcCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTa      <  1:1450738/69‑1 (MQ=255)
                                                                 ctcCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACa    >  1:1668692/1‑71 (MQ=255)
                                                                 ctcCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACa    >  1:1208614/1‑71 (MQ=255)
                                                                 ctcCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACa    >  1:1171809/1‑71 (MQ=255)
                                                                    cGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACATg  <  1:572838/70‑1 (MQ=255)
                                                                    cGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACATg  <  1:2882091/70‑1 (MQ=255)
                                                                     |                                                                    
GACTGGAACTGCTCCTGAATTTCAGCAACTTCTGCTTCTGCCTCGGAGATGATTTCGTGTTTCTTCTCCGGGATGACCATGTCATCGATACCAACAGATGCACCAGAACGCGCTGCATAGGCGAAGCCGGTGTACATG  >  Exported/2255029‑2255166

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: