Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,168,904 C→T A1315V (GCA→GTA)  rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,168,9040CT97.6% 149.6 / NA 41A1315V (GCA→GTA) rpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (21/19);  minor base G (1/0);  total (22/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCG  >  REL606/4168831‑4168974
                                                                         |                                                                      
gAATTGGAAGCGAACGGCAAANTGGGTGCAACTTACTCTCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTaa                                                                      <  2:269756/76‑1 (MQ=255)
gAATTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTaa                                                                      <  2:710485/76‑1 (MQ=255)
gAATTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTaa                                                                      <  1:36506/76‑1 (MQ=255)
 aaTTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAAc                                                                     <  1:1546846/76‑1 (MQ=255)
  aTTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAAcc                                                                    >  2:998464/1‑76 (MQ=255)
  aTTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAAcc                                                                    >  2:1119511/1‑76 (MQ=255)
   ttGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCg                                                                   <  1:322338/76‑1 (MQ=255)
   ttGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTATGGTAACCg                                                                   >  2:720967/1‑76 (MQ=255)
     ggAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAg                                                                 <  2:382350/76‑1 (MQ=255)
      gAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGt                                                                <  2:890631/76‑1 (MQ=255)
          cGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCtt                                                            <  2:22649/76‑1 (MQ=255)
            aaCGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCa                                                          <  1:40822/76‑1 (MQ=255)
                 cAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTcc                                                     <  1:1126395/76‑1 (MQ=255)
                 cAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTcc                                                     <  2:1281388/76‑1 (MQ=255)
                  aaaGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCcg                                                    <  1:1570846/76‑1 (MQ=255)
                  aaaGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCcg                                                    >  2:1565442/1‑76 (MQ=255)
                      tGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGc                                                >  1:223933/1‑76 (MQ=255)
                             aaCTTACTCCCGCGATCTGCTGGGTATCACTAAAGCGTCTCTGGGAACCGAGTCGTTCATCTCCGTGGCATCGGTc                                         >  1:1524968/1‑76 (MQ=255)
                             aaCTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTc                                         >  2:1202060/1‑76 (MQ=255)
                                 tACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAgg                                     <  1:987292/76‑1 (MQ=255)
                                   cTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGgag                                   >  1:852634/1‑76 (MQ=255)
                                     cccGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGAc                                 >  1:847569/1‑76 (MQ=255)
                                     cccGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGAc                                 >  2:266107/1‑76 (MQ=255)
                                      ccGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGAcc                                >  1:1570978/1‑76 (MQ=255)
                                       cgcgATCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCa                               >  2:1660973/1‑76 (MQ=255)
                                           aTCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTcg                           >  1:1522054/1‑76 (MQ=255)
                                           aTCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTcg                           >  2:623135/1‑76 (MQ=255)
                                           aTCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTcg                           <  2:411647/76‑1 (MQ=255)
                                           aTCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTcg                           <  2:1320002/76‑1 (MQ=255)
                                            tCTGCTGGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTcgc                          >  1:709292/1‑76 (MQ=255)
                                                 tgGGTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCt                     >  2:1612187/1‑76 (MQ=255)
                                                    gTATCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGAc                  <  1:1367553/76‑1 (MQ=255)
                                                      aTCACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCg                >  2:834623/1‑76 (MQ=255)
                                                       tCACCATAGCGTCTCTGGTAAACGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGa               <  2:702623/76‑1 (MQ=255)
                                                        cACCAAAGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGaa              >  1:702818/1‑76 (MQ=255)
                                                            aaaGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAagcag          <  2:1396216/76‑1 (MQ=255)
                                                              aGCGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGcc        >  1:1056989/1‑76 (MQ=255)
                                                                cGTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGt      >  2:116404/1‑76 (MQ=255)
                                                                 gTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGtt     <  2:356449/76‑1 (MQ=255)
                                                                 gTCTCTGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGtt     <  2:1389564/76‑1 (MQ=255)
                                                                   ctctGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGc   >  1:308661/1‑76 (MQ=255)
                                                                   ctctGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGc   <  1:379466/76‑1 (MQ=255)
                                                                    tctGGTAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCg  >  2:468122/1‑76 (MQ=255)
                                                                         |                                                                      
GAATTGGAAGCGAACGGCAAAGTGGGTGCAACTTACTCCCGCGATCTGCTGGGTATCACCAAAGCGTCTCTGGCAACCGAGTCCTTCATCTCCGCGGCATCGTTCCAGGAGACCACTCGCGTGCTGACCGAAGCAGCCGTTGCG  >  REL606/4168831‑4168974

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