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.7% 223.8 / ‑0.5 61intergenic (+100/‑95)dnaG/rpoDDNA primase/RNA polymerase sigma factor
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (24/35);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

ATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATG  >  REL606/3145696‑3145841
                                                                         |                                                                        
aTCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTtggtgg                                                                        >  2:807213/1‑76 (MQ=255)
 tCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTggtggt                                                                       <  2:1229439/76‑1 (MQ=255)
  cGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTa                                                                      <  1:702554/76‑1 (MQ=255)
   gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaa                                                                     <  1:1438909/76‑1 (MQ=255)
   gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaa                                                                     <  1:116099/76‑1 (MQ=255)
   gggAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTaa                                                                     <  1:1103/76‑1 (MQ=255)
     gAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaac                                                                   <  1:1659343/76‑1 (MQ=255)
      aaGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaaca                                                                  <  2:583530/76‑1 (MQ=255)
      aaGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaaca                                                                  >  2:1188067/1‑76 (MQ=255)
      aaGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaaca                                                                  <  2:1518497/76‑1 (MQ=255)
       aGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaacaa                                                                 <  1:934514/76‑1 (MQ=255)
        gCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAaacaac                                                                >  2:1741113/1‑76 (MQ=255)
           cccGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGt                                                             >  1:1122453/1‑76 (MQ=255)
             cGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTg                                                           <  1:1236960/76‑1 (MQ=255)
             cGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTg                                                           >  1:35694/1‑76 (MQ=255)
                cAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAt                                                        >  2:859844/1‑76 (MQ=255)
                cAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAt                                                        >  2:527283/1‑76 (MQ=255)
                cAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAt                                                        >  2:420154/1‑76 (MQ=255)
                 aGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAtt                                                       <  1:474285/76‑1 (MQ=255)
                  gCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAttt                                                      <  2:100033/76‑1 (MQ=255)
                  gCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAttt                                                      <  2:1667216/76‑1 (MQ=255)
                  gCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAttt                                                      <  1:1736781/76‑1 (MQ=255)
                  gCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGAttt                                                      <  1:967762/76‑1 (MQ=255)
                   ccGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTc                                                     <  2:452894/76‑1 (MQ=255)
                    cGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCa                                                    <  2:1031382/76‑1 (MQ=255)
                    cGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCa                                                    <  2:1615341/76‑1 (MQ=255)
                    cGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCa                                                    <  1:613316/76‑1 (MQ=255)
                     gCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAg                                                   <  1:626987/76‑1 (MQ=255)
                        cTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGt                                                <  1:386243/76‑1 (MQ=255)
                          gagaGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTa                                              <  1:1081816/76‑1 (MQ=255)
                           tgagGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTaa                                             >  2:1453723/2‑76 (MQ=255)
                            gagGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAAc                                            <  2:799498/76‑1 (MQ=255)
                                cAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCt                                        >  1:1093107/1‑76 (MQ=255)
                                cAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCt                                        >  1:1243042/1‑76 (MQ=255)
                                 aGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGcca                                       >  1:889763/1‑74 (MQ=255)
                                 aGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTg                                       >  1:573483/1‑76 (MQ=255)
                                 aGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTg                                       >  1:1415031/1‑76 (MQ=255)
                                 aGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTg                                       <  1:93729/76‑1 (MQ=255)
                                       aaaTATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGAc                                 <  1:848991/76‑1 (MQ=255)
                                       aaaTATATAAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGAc                                 >  1:3087/1‑76 (MQ=255)
                                        aaTATATAAGTACGCCCTCGTACTTATCGTTGGTGGTAAACAACCTTTGGATTTCAGCGTTAACGGCTAAAAGaaa                                >  2:1481550/1‑74 (MQ=255)
                                          tatataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATc                              <  2:209195/76‑1 (MQ=255)
                                          tatataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATc                              <  1:648096/76‑1 (MQ=255)
                                          tatataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATc                              <  2:1507039/76‑1 (MQ=255)
                                          tatataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATc                              <  1:950381/76‑1 (MQ=255)
                                          tatataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATc                              <  2:946337/76‑1 (MQ=255)
                                           atataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCg                             >  1:873074/1‑76 (MQ=255)
                                            tataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCgg                            <  1:449133/76‑1 (MQ=255)
                                            tataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCgg                            <  2:776904/76‑1 (MQ=255)
                                             ataAGTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCggg                           >  2:1661384/1‑76 (MQ=255)
                                                aGTACGCCCTCGTAATTATCGTTGGTGTTAAACAACCGTTGGATTTCAGCGTTAACGGCTGACGGACATCGGGTCa                        <  2:1024206/76‑1 (MQ=255)
                                                 gTACGCCCTCGTAATTATCGTTGGTGGTAAACAACCGTTGGGTTTCAGCGTTAACGGCTGCAGGACATCGAGGCaa                       >  2:712215/1‑76 (MQ=255)
                                                          cGTAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAac              >  2:1385306/1‑76 (MQ=255)
                                                           gTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAaca             <  1:1622594/76‑1 (MQ=255)
                                                            tAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAacac            >  2:344145/1‑76 (MQ=255)
                                                            tAATTATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAacac            >  2:956806/1‑76 (MQ=255)
                                                               ttATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCaa         <  2:1219763/76‑1 (MQ=255)
                                                                tATCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAAc        >  1:1197762/1‑76 (MQ=255)
                                                                 aTCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAAcc       >  2:192868/1‑76 (MQ=255)
                                                                 aTCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAAcc       >  1:983705/1‑76 (MQ=255)
                                                                  tCGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCt      <  1:1115955/76‑1 (MQ=255)
                                                                   cGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTc     <  1:96137/76‑1 (MQ=255)
                                                                   cGTTGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTc     <  2:843380/76‑1 (MQ=255)
                                                                     ttGGTGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCAt   >  2:1125398/1‑76 (MQ=255)
                                                                      tggtggTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATg  <  2:212490/76‑1 (MQ=255)
                                                                      tggtggTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATg  <  2:1665869/76‑1 (MQ=255)
                                                                         |                                                                        
ATCGGGAAGCCCCCGACAGCCGCACTGAGAGGCAGCGGCAAATATATAAGTACGCCCTCGTAATTATCGTTGGCGGTAAACAACCGTTGGATTTCAGCGTTAACGGCTGAAGGACATCGGGTCAATCGCCCAACACCAACCTCATG  >  REL606/3145696‑3145841

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