Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,552,005 G→A 14.3% intergenic (+55/+126) pydA → / ← rutR NADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,552,0050GA14.3% 61.4 / 4.3 28intergenic (+55/+126)pydA/rutRNADP‑dependent dihydropyrimidine dehydrogenase subunit PreA/transcriptional regulator, TetR family
Reads supporting (aligned to +/- strand):  ref base G (10/14);  new base A (2/2);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

TGCATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCC  >  NC_002947/4551864‑4552150
                                                                                                                                             |                                                                                                                                                 
tGCATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCggg                                                                                                                                           >  1:389411/1‑150 (MQ=255)
  cATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGcc                                                                                                                                         <  2:760569/150‑1 (MQ=255)
   aTCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCg                                                                                                                                        >  2:451364/1‑150 (MQ=255)
         aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc                                                                                                                                  <  1:405279/150‑1 (MQ=255)
         aTGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGc                                                                                                                                  <  1:47584/150‑1 (MQ=255)
                cacaGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCca                                                                                                                           <  1:726917/150‑1 (MQ=255)
                   aGGATACCGGCAAGCCGTACCTGAACGGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGTGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCcagca                                                                                                                        <  2:360646/150‑1 (MQ=255)
                     gATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGt                                                                                                                      <  2:659939/150‑1 (MQ=255)
                      aTACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTc                                                                                                                     <  2:389411/150‑1 (MQ=255)
                            gCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGcc                                                                                                               >  2:570091/1‑150 (MQ=255)
                                    tACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACa                                                                                                       <  1:570091/150‑1 (MQ=255)
                                                           gcgTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCgg                                                                                <  1:451364/150‑1 (MQ=255)
                                                              caaCCCCTACCGTGAGGCCAGTTGAAGCTGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGccc                                                                             <  2:694265/149‑1 (MQ=255)
                                                                  cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCg                                                                                        <  1:548985/135‑1 (MQ=255)
                                                                  cccTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCg                                                                                        >  2:548985/1‑135 (MQ=255)
                                                                         gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTcg                                                                                                                                                   <  1:538663/69‑1 (MQ=255)
                                                                         gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTcg                                                                                                                                                   >  2:538663/1‑69 (MQ=255)
                                                                         gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCTCTTCGCGGCCCGATCGCGac                                                                    >  2:241024/1‑148 (MQ=255)
                                                                         gTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCTCTTCGCGGCCCGATCGCGac                                                                    <  1:241024/148‑1 (MQ=255)
                                                                          tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc                                                                                                                                  <  2:388649/85‑1 (MQ=255)
                                                                          tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc                                                                                                                                  <  2:388727/85‑1 (MQ=255)
                                                                          tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc                                                                                                                                  >  1:388727/1‑85 (MQ=255)
                                                                          tGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCACGATCGGGCCGCACAGc                                                                                                                                  >  1:388649/1‑85 (MQ=255)
                                                                                   ttGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCt                                                         >  2:494376/1‑149 (MQ=255)
                                                                                   ttGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCt                                                         <  1:494376/149‑1 (MQ=255)
                                                                                                     tGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                              <  1:204169/94‑1 (MQ=255)
                                                                                                     tGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAg                                                                                              >  2:204169/1‑94 (MQ=255)
                                                                                                                                    ccTTGTGTCGCGATCGGGCCGCACAGCGGCCCCACCGGTCTTTGCCCCGATACCGACATCGAGGGGCCAGGTCGCGGCCCGATCGCCACCCAACGCCGCTCCCACAAAAAGACCACCTCCACCTCGAATGCCGGATCACGGCTCCCGccc       >  1:748043/1‑150 (MQ=255)
                                                                                                                                      ttGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCAGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGCCTCCAGCCCGa     >  2:368035/1‑150 (MQ=255)
                                                                                                                                         tgtCGCGATCGGGCCGCACATCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACCCAAGGCCGCTCCCACAAAAAGACCTCGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATcc  >  2:315476/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
TGCATCGAAATGGAGGCACAGGATACCGGCAAGCCGTACCTGAACTGGACGCAGGATCCGCGTAACCCCTACCGTGAGGCCAGTTGAAGCGGGTCACCGCCTGAACAGCTGACGCGGCCCCTGTGGGAGCGGCCTTGTGTCGCGATCGGGCCGCACAGCGGCCCCAGCAGTCTTTGCCGCGATACAGACATCGAGGGGCCGCTTCGCGGCCCGATCGCGACACAAGGCCGCTCCCACAAAAAGACCACGTCCACCTGGAATGCCGGATCACGGCTCCAGCCCGATCC  >  NC_002947/4551864‑4552150

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: