Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,241,537 C→T 32.2% intergenic (+23/‑69) PP_1080 → / → grxD molybdopterin‑binding oxidoreductase/monothiol glutaredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,241,5370CT32.2% 31.9 / 21.9 31intergenic (+23/‑69)PP_1080/grxDmolybdopterin‑binding oxidoreductase/monothiol glutaredoxin
Reads supporting (aligned to +/- strand):  ref base C (8/13);  new base T (5/5);  total (13/18)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGC  >  NC_002947/1241390‑1241675
                                                                                                                                                   |                                                                                                                                          
aCCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTTCCGGTCGAAGCCGCTTGAGGAACGGCAAGGCCGAGCAAGgcgc                                                                                                                                          >  2:5638/1‑150 (MQ=255)
aCCTGCAGATCGCCCTGGCGCAACCAGGAGTCACTGGTTATGACCTGACTGATTATCGCGTTCGGGACGTGGTTTACACGCACGCCGAAGACGAGGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAGGGCCTGTCACGgcgc                                                                                                                                          >  1:578316/1‑150 (MQ=255)
              cTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATa                                                                                                                            >  2:482562/1‑150 (MQ=255)
                   gcAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATg                                                                                                                       >  2:709729/1‑150 (MQ=255)
                      aCCCGGCGTAAGGGGCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGa                                                                                                                    <  2:335652/150‑1 (MQ=255)
                                        tgacctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGccttgc                                                                                                                                <  2:3832/120‑6 (MQ=255)
                                        tgacctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCGCCGTGCTCGGccttgc                                                                                                                                >  1:3832/1‑115 (MQ=255)
                                        tgacctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACtgt                                                                                                  <  1:279579/150‑1 (MQ=255)
                                           cctgaccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTc                                                                                               <  1:482562/150‑1 (MQ=255)
                                                ccGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGAGTACCCGTGCAGGCTGAGTGCTCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATa                                                                                          <  1:412817/150‑1 (MQ=255)
                                                      aGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAg                                                                                    <  1:784401/150‑1 (MQ=255)
                                                           tttCGGGACCAGGTTTCCGGAAACGCTGCGGTGAATGGGGTACCCGTGCAGGTCGAGGCGTCTGGACGTGGGTCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAg                                                                               <  1:412979/150‑1 (MQ=255)
                                                            ttCGGGACCAGGTTTCCGGAAACGCTGCGGTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAgg                                                                              <  1:400039/150‑1 (MQ=255)
                                                                                     tGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGc                                                                                   >  1:2630/1‑120 (MQ=255)
                                                                                     tGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGc                                                                                   <  2:2630/120‑1 (MQ=255)
                                                                                          tCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCaaa                                                <  2:578316/150‑1 (MQ=255)
                                                                                                        gTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGGGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCcaa                                  <  1:5638/150‑1 (MQ=255)
                                                                                                              gTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                  <  1:289888/112‑1 (MQ=255)
                                                                                                              gTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                  >  2:289888/1‑112 (MQ=255)
                                                                                                                cacggtgtcgcattgtatcggaaaGCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCa                          <  1:286677‑M1/126‑1 (MQ=255)
                                                                                                                    gCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATCCAATGTGATACCAGACCAACAACCTGGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCt                      <  1:412936/150‑1 (MQ=255)
                                                                                                                       gCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGAATGTGTCGGATAGTTCAGACGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGACAACAAAACAATTATgat                   >  2:289353/1‑148 (MQ=255)
                                                                                                                            aGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACa              >  2:85433/1‑150 (MQ=255)
                                                                                                                                 tcggaaaGCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACAt                                                                      <  2:294921‑M1/82‑1 (MQ=255)
                                                                                                                                 tcggaaaGCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACAt                                                                      >  1:294921‑M1/8‑89 (MQ=255)
                                                                                                                                    gaaaGCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGc                                     >  1:603961‑M1/5‑119 (MQ=255)
                                                                                                                                    gaaaGCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGc                                     <  2:603961‑M1/115‑1 (MQ=255)
                                                                                                                                    gaaaGCCGAGCGCCGTGCTCGGCCTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                  <  2:395045‑M1/86‑1 (MQ=255)
                                                                                                                                    gaaaGCCGAGCGCCGTGCTCGGCCTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGat                                                                  >  1:395045‑M1/5‑90 (MQ=255)
                                                                                                                                     aaGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCt     >  1:747864/1‑150 (MQ=255)
                                                                                                                                        GCCGAGCGCCGTGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGc  >  2:286677‑M1/1‑150 (MQ=255)
                                                                                                                                        GCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGa                                                                                            >  2:177783‑M1/1‑60 (MQ=255)
                                                                                                                                        GCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGa                                                                                            <  1:177783‑M1/60‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
ACCTGCAGATCGCCCTGGCGCAACCCGGCGTGAGTGCCAATGACCTGACCGATGAGCGCTTTCGGGACCAGGTTTCCGGAAACGCTGCGCTCAATGGCGTACCCGTGCAGGTCGAGGCCGCTTGAGGAACGGCAAGGCCGAGCACGGCGCTCGGCTTTCCGATACAATGCGACACCGTGCCGACGACTGTGTCGGATAGTTCAGCCGAGGTGCTACATGGATATCATCGAAACGATCAAAGAGCAGATTGCCAACAACACCATTCTGCTTTACATGAAAGGCTCGC  >  NC_002947/1241390‑1241675

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