Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,840,679 Δ1 bp 11.8% coding (659/666 nt) 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,6790T.11.8% 90.2 / 7.4 34coding (659/666 nt)PP_5119TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (14/15);  new base . (2/2);  total (16/18)
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.67e-01

CCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAATGTT  >  NC_002947/5840537‑5840826
                                                                                                                                               |                                                                                                                                                   
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                                                                             >  2:684183/1‑122 (MQ=255)
                                                                                             ggggCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGa                                                                             <  1:684183/122‑1 (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                                                                                          <  1:843384/90‑1 (MQ=255)
                                                                                                                aTCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACgg                                                                                          >  2:843384/1‑90 (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                                                                                <  2:513019‑M2/87‑70 (MQ=255)
                                                                                                                             gagCCCTGGGCTAGAGCTccagggctcttgtgccatcgatcaggagattgtcatgccccttgcgcaattgatcaccccgcagcagtt                                                                                >  1:513019‑M2/1‑18 (MQ=255)
                                                                                                                                          agagCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAAt     <  2:576268/150‑1 (MQ=255)
                                                                                                                                             gCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAATGtt  <  1:303348/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CCACATCCTCGAGGGCAAAGCGGCAGTCGAGGATCACCAACTTGGGCGAGCCCAGACGCTCGGCCAACTGCTGCGGGGTGATCAATTGCGCAAGGGGCATGACAATCTCCTGATCGATGGCACAAGAGCCCTGGGCTAGAGCTCCAGGGCCTGGTTGAACGGTACATGGAATTCCTGACACAGCGCTTCTACGGCAGTACGGGCATTGTCGGTGACGAAACCGAGCTCCAGCACCAGCACCTGATAGACGCCGCGCTTGAAGGCTTCTTCACCCAAGTGCGCAGAATGTT  >  NC_002947/5840537‑5840826

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: