Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 564,380:1 +A 14.6% intergenic (+29/‑144) rplQ → / → katA 50S ribosomal protein L17/catalase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947564,3801.A14.6% 158.0 / 24.2 48intergenic (+29/‑144)rplQ/katA50S ribosomal protein L17/catalase
Reads supporting (aligned to +/- strand):  ref base . (17/24);  new base A (3/4);  total (20/28)
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.91e-01

CCGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGA  >  NC_002947/564231‑564527
                                                                                                                                                      |                                                                                                                                                   
ccGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGc                                                                                                                                                      <  1:967106/150‑1 (MQ=255)
  gTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCt                                                                                                                                                   >  2:885373/1‑150 (MQ=255)
           gcTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCgg                                                                                                                                          >  1:408616/1‑150 (MQ=255)
                 tGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGtttttt                                                                                                                                    <  1:496969/150‑1 (MQ=255)
                 tGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAACCCGGGCACTAGGCCCGGttttct                                                                                                                                    <  1:1457360/150‑3 (MQ=255)
                  gcgcATCCTGAAGTGCGGTTTCCGCTCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTg                                                                                                                                   >  2:1380854/1‑150 (MQ=255)
                          tGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACa                                                                                                                           >  1:65640/1‑150 (MQ=255)
                           gAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACAt                                                                                                                          <  2:65640/150‑1 (MQ=255)
                                      tCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTAt                                                                                                               <  1:1110116/150‑1 (MQ=255)
                                         gcgcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGccc                                                                                                                                            >  1:452113/1‑118 (MQ=255)
                                         gcgcTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGccc                                                                                                                                            <  2:452113/118‑1 (MQ=255)
                                              ggCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCTGTTTTTTGTTTATACATCTATTGTTTATTTCTATTg                                                                                                       <  2:491399/150‑1 (MQ=255)
                                               gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTa                                                                                                                <  2:333224/140‑1 (MQ=255)
                                               gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTa                                                                                                                >  1:333224/1‑140 (MQ=255)
                                               gCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGa                                                                                                      <  1:620204/150‑1 (MQ=255)
                                                 gACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATg                                                                                                    <  1:885373/150‑1 (MQ=255)
                                                    aaCGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATTTaa                                                                                                 <  2:405972/150‑1 (MQ=255)
                                                    aaCGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATCGTTTATTTCGATTGATGTaa                                                                                                 <  2:920102/150‑1 (MQ=255)
                                                                        cGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                             <  1:1193210/150‑1 (MQ=255)
                                                                        cGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                             <  1:1047062/150‑1 (MQ=255)
                                                                        cGAGCTGGTTGATCGTCCGGTCGGCGGTGCAGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTcc                                                                             <  1:1209449/150‑1 (MQ=255)
                                                                          aGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCt                                                                           <  2:685856/150‑1 (MQ=255)
                                                                             tGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTa                                                                        >  2:1207657/1‑150 (MQ=255)
                                                                                  gATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCtt                                                                   <  1:1207657/150‑1 (MQ=255)
                                                                                    tCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGc                                                                 >  1:941797/1‑150 (MQ=255)
                                                                                                   tGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                  <  2:1315148/150‑1 (MQ=255)
                                                                                                     cTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctctct                                                >  1:896797/1‑150 (MQ=255)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACa                                             >  2:539123/1‑150 (MQ=255)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGTTTTCTGta                                                                                                                                 >  2:838744/1‑65 (MQ=17)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGTTTTCTGta                                                                                                                                 >  2:838473/1‑65 (MQ=17)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGTTTTCTGta                                                                                                                                 <  1:838744/66‑2 (MQ=17)
                                                                                                        tAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAAAAAACCGGGCACTAGGCCCGGTTTTCTGta                                                                                                                                 <  1:838473/66‑2 (MQ=17)
                                                                                                             gctgctGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGa                                        <  2:896797/150‑1 (MQ=255)
                                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                       >  1:1015665/1‑110 (MQ=255)
                                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTaa                                                                       <  2:1015665/110‑1 (MQ=255)
                                                                                                                      tAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGc                               <  1:539123/150‑1 (MQ=255)
                                                                                                                             ttCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                        >  1:163155/1‑150 (MQ=255)
                                                                                                                                  tctgCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCg                   >  1:971508/1‑150 (MQ=255)
                                                                                                                                      cTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                  >  1:1182092/1‑115 (MQ=255)
                                                                                                                                      cTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGCTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATActctct                                                  <  2:1182092/115‑1 (MQ=255)
                                                                                                                                       aacaaaaaaCCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                        <  1:1094786/135‑1 (MQ=255)
                                                                                                                                       aaCAAAAAACCGGGCACTAGGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGaa                        >  2:1094786/2‑140 (MQ=255)
                                                                                                                                           gAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCCCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGgaga          <  2:941797/150‑1 (MQ=255)
                                                                                                                                            aaaaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAAt               >  1:574004/1‑144 (MQ=255)
                                                                                                                                            aaaaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAAt               <  2:574004/144‑1 (MQ=255)
                                                                                                                                              aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGagag         >  1:945171/1‑148 (MQ=255)
                                                                                                                                              aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGagag         <  2:945171/148‑1 (MQ=255)
                                                                                                                                              aaCCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGcc       >  1:1633294/1‑150 (MQ=255)
                                                                                                                                                   ggCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGa  <  1:914371/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CCGTCAGGGCGGCTACCTGCGCATCCTGAAGTGCGGTTTCCGCGCTGGCGACAACGCGCCTATGGCGTATGTCGAGCTGGTTGATCGTCCGGTCGGCGGTGCTGTAGAAGCTGCTGAGTAAGACGTTCTGTCTGCTACAGAAAACCGGGCCTAGTGCCCGGTTTTTTGTTTATACATCTATTGTTTATTTCTATTGATGTAAGTCATGTAATGAATTTGTCCCTGTAATCTTGCTCGTCAATACTCTCTCTACAGCCGATTAGCCGGCCCTTGAAGACCGAAATGGAGAGCCCATGA  >  NC_002947/564231‑564527

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