Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,840,675 G→T 12.5% L221L (CTC→CTA PP_5119 ← TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,840,6750GT12.5% 81.4 / 3.3 32L221L (CTC→CTAPP_5119TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (14/14);  new base T (2/2);  total (16/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 = 8.67e-01

TGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAAT  >  NC_002947/5840527‑5840823
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                         <  2:637492‑M2/124‑11 (MQ=255)
                                       aGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGGCCCTGGAGCTCTAGCcagggctctt                                                                                                                                         >  1:637492‑M2/1‑114 (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)
                                                                                              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                                                                          <  1:684183/122‑1 (MQ=255)
                                                                                                       ggggCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGa                                                                          >  2:684183/1‑122 (MQ=255)
                                                                                                            aTGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTccagg                                                                                                                                             >  1:713140‑M2/1‑45 (MQ=255)
                                                                                                            aTGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTccagg                                                                                                                                             <  2:713140‑M2/50‑6 (MQ=255)
                                                                                                              gACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGAGTGGCTTCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGGACCAGCACCTGATAGAcg                                       >  1:53464/1‑150 (MQ=255)
                                                                                                                          aTCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg                                                                                       >  2:843384/1‑90 (MQ=255)
                                                                                                                          aTCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg                                                                                       <  1:843384/90‑1 (MQ=255)
                                                                                                                              aTGGCACAAGAGCCAAGGGCTAGAGGGCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACACCCATCCGACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGctt                       <  1:15419/150‑1 (MQ=255)
                                                                                                                                   acaAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCa                  <  1:676667/150‑1 (MQ=255)
                                                                                                                                   acaAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCa                  >  2:303348/1‑150 (MQ=255)
                                                                                                                                     aaGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCAcc                >  1:32136/1‑150 (MQ=255)
                                                                                                                                       gagCCCTGGGCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt                                                                             >  1:513019‑M2/1‑18 (MQ=255)
                                                                                                                                       gagCCCTGGGCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt                                                                             <  2:513019‑M2/87‑70 (MQ=255)
                                                                                                                                                    agagCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAAt  <  2:576268/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCCATAGTCCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAAT  >  NC_002947/5840527‑5840823

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