Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,145,769 C→T intergenic (+100/‑95) dnaG → / → rpoD DNA primase/RNA polymerase sigma factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,145,7690CT96.0% 267.2 / 1.5 74intergenic (+100/‑95)dnaG/rpoDDNA primase/RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (37/34);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

TCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCAT  >  REL606/3145697‑3145840
                                                                        |                                                                       
tCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTggtggt                                                                      <  2:401913/76‑1 (MQ=255)
tCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTggtggt                                                                      <  2:780751/76‑1 (MQ=255)
 cGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTa                                                                     <  1:851896/76‑1 (MQ=255)
  gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaa                                                                    <  2:1428270/76‑1 (MQ=255)
  gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaa                                                                    <  2:187781/76‑1 (MQ=255)
  gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaa                                                                    <  1:223667/76‑1 (MQ=255)
  gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTaa                                                                    >  1:1173067/1‑76 (MQ=255)
   tgAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaaa                                                                   >  2:925684/2‑76 (MQ=255)
   ggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaaa                                                                   <  1:1413988/76‑1 (MQ=255)
    gAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaac                                                                  <  2:1373087/76‑1 (MQ=255)
      aGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaacaa                                                                <  2:1915939/76‑1 (MQ=255)
       gCCCCCGACAGCCGCACTGAGAGACAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaacaac                                                               >  2:344676/1‑76 (MQ=255)
        cccccGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAAcc                                                              >  2:1675016/1‑76 (MQ=255)
        cccccGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAAcc                                                              <  1:1580877/76‑1 (MQ=255)
            cGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACACCCTCGTAATTATCGTTGGTGGTAAACAACCGTTg                                                          >  1:2053835/1‑76 (MQ=255)
                aGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAtt                                                      <  2:858003/76‑1 (MQ=255)
                  ccGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTc                                                    <  2:1929335/76‑1 (MQ=255)
                  ccGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTc                                                    <  1:2071920/76‑1 (MQ=255)
                  ccGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTc                                                    <  2:582485/76‑1 (MQ=255)
                  ccGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTc                                                    <  2:571805/76‑1 (MQ=255)
                    gCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAg                                                  <  1:1651491/76‑1 (MQ=255)
                    gCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAg                                                  <  1:1651482/76‑1 (MQ=255)
                     cACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGc                                                 <  1:185672/76‑1 (MQ=255)
                        tGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGtt                                              >  1:335762/1‑76 (MQ=255)
                         gagaGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTa                                             >  2:1947221/1‑76 (MQ=255)
                         gagaGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTa                                             >  2:1730243/1‑76 (MQ=255)
                          agagGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTaa                                            >  2:1727245/1‑76 (MQ=255)
                           gagGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAAc                                           <  1:1012140/76‑1 (MQ=255)
                            agGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACg                                          >  2:447431/1‑76 (MQ=255)
                                  cGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGaa                                    <  1:554839/76‑1 (MQ=255)
                                   ggCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAg                                   >  2:556531/1‑76 (MQ=255)
                                   ggCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAg                                   <  1:884715/76‑1 (MQ=255)
                                   ggCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAg                                   >  1:717046/1‑76 (MQ=255)
                                     cAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGa                                 >  2:1003604/1‑76 (MQ=255)
                                     cAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGa                                 >  2:192196/1‑76 (MQ=255)
                                      aaaTATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGAc                                >  2:1923284/1‑76 (MQ=255)
                                        atatatAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACAt                              >  1:1027746/1‑76 (MQ=255)
                                         tatataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATc                             <  2:102296/76‑1 (MQ=255)
                                          atataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTACCGGCTGAAGGACATCg                            >  1:1646548/1‑76 (MQ=255)
                                          atataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCg                            >  2:827730/1‑76 (MQ=255)
                                           tataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCgg                           <  1:551447/76‑1 (MQ=255)
                                           tataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCgg                           <  1:169560/76‑1 (MQ=255)
                                            ataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCggg                          >  1:757418/1‑76 (MQ=255)
                                            ataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCggg                          <  2:1447220/76‑1 (MQ=255)
                                            ataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCggg                          >  2:831830/1‑76 (MQ=255)
                                             taAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGt                         >  1:502150/1‑76 (MQ=255)
                                               aGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCa                       <  1:766829/76‑1 (MQ=255)
                                               aGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCa                       <  2:1242282/76‑1 (MQ=255)
                                               aGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCa                       <  2:227689/76‑1 (MQ=255)
                                               aGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCa                       >  2:333142/1‑76 (MQ=255)
                                                  aCGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATc                    >  1:571020/1‑76 (MQ=255)
                                                   cGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCg                   <  1:925684/76‑1 (MQ=255)
                                                   cGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCg                   <  1:834660/76‑1 (MQ=255)
                                                    gCCCCCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGc                  >  1:73092/1‑76 (MQ=255)
                                                     cccTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGcc                 >  2:180581/1‑76 (MQ=255)
                                                       cTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCa               >  1:2049592/1‑76 (MQ=255)
                                                        tCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCaa              >  2:1074155/1‑76 (MQ=255)
                                                          gTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAaca            >  2:164172/1‑76 (MQ=255)
                                                          gTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAaca            >  1:1729906/1‑76 (MQ=255)
                                                          gTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAaca            <  2:1244965/76‑1 (MQ=255)
                                                          gTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAaca            >  1:1428385/1‑76 (MQ=255)
                                                           tAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAacac           >  1:1554567/1‑76 (MQ=255)
                                                           tAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAacac           >  2:1088198/1‑76 (MQ=255)
                                                           tAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAacac           >  1:1414013/1‑76 (MQ=255)
                                                             aTTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCa         >  1:1020687/1‑76 (MQ=255)
                                                             aTTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCa         <  2:1985427/76‑1 (MQ=255)
                                                             aTTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCa         <  2:1605602/76‑1 (MQ=255)
                                                             aTTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCa         <  2:1515887/76‑1 (MQ=255)
                                                              ttATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCaa        >  2:999068/1‑76 (MQ=255)
                                                              ttATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCac        >  2:48865/1‑75 (MQ=255)
                                                               tATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAAACGCCCAACACCAAc       <  2:2053835/76‑1 (MQ=255)
                                                                aTCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAAcc      <  2:1033088/76‑1 (MQ=255)
                                                                 tCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCt     >  2:568858/1‑76 (MQ=255)
                                                                  cGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTc    >  2:80732/1‑76 (MQ=255)
                                                                  cGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTc    >  1:1782092/1‑76 (MQ=255)
                                                                  cGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTc    >  2:942351/1‑76 (MQ=255)
                                                                    ttGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCAt  <  1:630006/76‑1 (MQ=255)
                                                                        |                                                                       
TCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCAT  >  REL606/3145697‑3145840

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