Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,840,670:1 +A 10.8% intergenic (‑34/+2) sseA ← / ← PP_5119 3‑mercaptopyruvate sulfurtransferase/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,840,6701.A10.8% 125.5 / 11.5 37intergenic (‑34/+2)sseA/PP_51193‑mercaptopyruvate sulfurtransferase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (15/18);  new base A (2/2);  total (17/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

TGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAA  >  NC_002947/5840522‑5840812
                                                                                                                                                     |                                                                                                                                              
tGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTggg                                                                                                                                               <  2:845826/150‑1 (MQ=255)
     tGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGACCTGGGCTaga                                                                                                                                          >  1:177849/1‑150 (MQ=255)
               ccACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTccagggc                                                                                                                                >  2:456616‑M2/1‑143 (MQ=255)
                cacaTCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGCGCCCTGGGCTCGAGCTccggcgcc                                                                                                                               <  2:607191‑M2/150‑9 (MQ=255)
                           gggCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTgg                                                                                                                            >  1:443574/1‑142 (MQ=255)
                           gggCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTgg                                                                                                                            <  2:443574/142‑1 (MQ=255)
                              cAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCAAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTAc                                                                                                                 <  1:612705/150‑1 (MQ=255)
                                aaGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACAt                                                                                                               <  2:332304/150‑1 (MQ=255)
                                 aGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATg                                                                                                              <  2:452423/150‑1 (MQ=255)
                                           gAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGa                                                                                                    <  1:745050/150‑1 (MQ=255)
                                            aGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGGCCCTGGAGCTCTAGCcagggctctt                                                                                                                              >  1:637492‑M2/1‑114 (MQ=255)
                                            aGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGGCCCTGGAGCTCTAGCcagggctctt                                                                                                                              <  2:637492‑M2/124‑11 (MQ=255)
                                                   ccAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGAGCTCTAGCcagggctcttgtgccatcgatcaggagattg                                                                                                         >  2:77326‑M2/1‑108 (MQ=255)
                                                   ccAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGAGCTCTAGCcagggctcttgtgccatcgatcaggagattg                                                                                                         <  1:77326‑M2/139‑32 (MQ=255)
                                                             cGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGTTACATGGAATTCCTGACACAGCTCTTCTACGGCa                                                                                  <  2:177849/150‑1 (MQ=255)
                                                                gCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTa                                                                               >  2:676667/1‑150 (MQ=255)
                                                                     gACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGc                                                                          >  1:576268/1‑150 (MQ=255)
                                                                                aaCTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGagctctagcccagggctc                                                                                                                               <  2:513019‑M1/87‑19 (MQ=255)
                                                                                aaCTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGagctctagcccagggctc                                                                                                                               >  1:513019‑M1/1‑69 (MQ=255)
                                                                                                   aTTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcac                                            <  1:208420/150‑1 (MQ=255)
                                                                                                   aTTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGATTTTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcac                                            >  1:520386/1‑150 (MQ=255)
                                                                                                    ttGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcacc                                           >  2:578924/1‑150 (MQ=255)
                                                                                                        gcAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTccagcaccagca                                       >  2:647785/1‑150 (MQ=255)
                                                                                                            ggggCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGa                                                               >  2:684183/1‑122 (MQ=255)
                                                                                                            ggggCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGa                                                               <  1:684183/122‑1 (MQ=255)
                                                                                                                 aTGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTccagg                                                                                                                                  <  2:713140‑M2/50‑6 (MQ=255)
                                                                                                                 aTGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTccagg                                                                                                                                  >  1:713140‑M2/1‑45 (MQ=255)
                                                                                                                   gACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGAGTGGCTTCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGGACCAGCACCTGATAGAcg                            >  1:53464/1‑150 (MQ=255)
                                                                                                                     cAATCTCCTGATCGATGGCACAAGAGCCCTGGgctagagctccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttgg                                      <  1:77326‑M1/139‑108 (MQ=255)
                                                                                                                     cAATCTCCTGATCGATGGCACAAGAGCCCTGGgctagagctccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttgg                                      >  2:77326‑M1/1‑32 (MQ=255)
                                                                                                                               aTCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg                                                                            <  1:843384/90‑1 (MQ=255)
                                                                                                                               aTCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg                                                                            >  2:843384/1‑90 (MQ=255)
                                                                                                                                   aTGGCACAAGAGCCAAGGGCTAGAGGGCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACACCCATCCGACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGctt            <  1:15419/150‑1 (MQ=255)
                                                                                                                                        acaAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCa       >  2:303348/1‑150 (MQ=255)
                                                                                                                                        acaAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCa       <  1:676667/150‑1 (MQ=255)
                                                                                                                                          aaGAGCCCTGGgctagagctccagggccttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcct                                <  2:637492‑M1/124‑114 (MQ=255)
                                                                                                                                          aaGAGCCCTGGgctagagctccagggccttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcct                                >  1:637492‑M1/1‑11 (MQ=255)
                                                                                                                                          aaGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCAcc     >  1:32136/1‑150 (MQ=255)
                                                                                                                                            gagCCCTGGGCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt                                                                  <  2:513019‑M2/87‑70 (MQ=255)
                                                                                                                                            gagCCCTGGGCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt                                                                  >  1:513019‑M2/1‑18 (MQ=255)
                                                                                                                                             ggcgccggagctcgagcccagggcgcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcctcgactgccgctttgccctcgaggatgtg   <  2:607191‑M1/144‑143 (MQ=255)
                                                                                                                                              gCCCTGGagctctagcccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagttggccgagcgtctgggctcgcccaagttggtgatcctcgactgccgctttgccctcgaggatgtgg  >  2:456616‑M1/1‑7 (MQ=255)
                                                                                                                                                ccTGGagctctagcccagggctcttgtgccatcgatcaggagattgtcat                                                                                                    <  2:713140‑M1/50‑46 (MQ=255)
                                                                                                                                                ccTGGagctctagcccagggctcttgtgccatcgatcaggagattgtcat                                                                                                    >  1:713140‑M1/1‑5 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
TGCGCTGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAA  >  NC_002947/5840522‑5840812

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