Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,392,927 C→G 12.5% intergenic (+97/‑89) PP_1214 → / → ruvC transcriptional regulator/crossover junction endodeoxyribonuclease RuvC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,392,9270CG12.5% 80.3 / 2.4 31intergenic (+97/‑89)PP_1214/ruvCtranscriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base C (15/11);  new base G (2/2);  total (18/13)
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.77e-01

TCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAC  >  NC_002947/1392781‑1393073
                                                                                                                                                   |                                                                                                                                                  
tCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCa                                                                                                                                                   >  4:318286/1‑148 (MQ=255)
     tCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAAtgtg                                                                                                                                             <  8:201410/149‑1 (MQ=255)
             ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTa                                                                                                                                     >  2:272991/1‑149 (MQ=255)
              cAGATTTCCGACGAGATCATGGAATACCTCGGCTAAGGTTATTTCGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAc                                                                                                                                    <  3:10572/149‑1 (MQ=255)
                     ccGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                              <  2:250068/148‑1 (MQ=255)
                                 tGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAcc                                                                                                                 <  1:272991/149‑1 (MQ=255)
                                 tGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAcc                                                                                                                 <  3:318286/149‑1 (MQ=255)
                                  ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGc                                                                                                                          <  5:106360/139‑1 (MQ=17)
                                  ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGc                                                                                                                          >  6:106360/1‑139 (MQ=17)
                                         cTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                          >  1:62796/1‑132 (MQ=255)
                                         cTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                          <  2:62796/132‑1 (MQ=255)
                                                       tttGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                        <  2:23463/120‑1 (MQ=255)
                                                       tttGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                        >  1:23463/1‑120 (MQ=255)
                                                                                       tAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaa                                                                                                                                                  <  5:143680/62‑1 (MQ=255)
                                                                                       tAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCaa                                                                                                                                                  >  6:143680/1‑62 (MQ=255)
                                                                                                               gggTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAcc                                   >  6:79038/1‑149 (MQ=255)
                                                                                                                ggTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                                  >  3:36528/1‑149 (MQ=255)
                                                                                                                      cccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTc                            >  2:230271/1‑149 (MQ=255)
                                                                                                                       ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                           >  1:58206/1‑149 (MQ=255)
                                                                                                                           tCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                          <  2:203103/146‑1 (MQ=18)
                                                                                                                           tCCTACAGGTACTGCACCACATTTGAAAAGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                          >  1:203103/1‑146 (MQ=18)
                                                                                                                                 aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCg                 >  2:307351/1‑149 (MQ=255)
                                                                                                                                 aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCg                 >  5:160711/1‑149 (MQ=255)
                                                                                                                                 aGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCg                 <  5:79038/149‑1 (MQ=255)
                                                                                                                                     aTTACGCCCATTTCAATTGTGGGTCTGTACGTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCt              <  4:36528/148‑1 (MQ=39)
                                                                                                                                          gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTc                                               >  5:328329/1‑110 (MQ=255)
                                                                                                                                          gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTc                                               <  6:328329/110‑1 (MQ=255)
                                                                                                                                          gCCCATATCAAATCTGCTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGATGTTTGAACAGATGCATCATTTTTTCGACCGCCTTAAACACACTGATTCTAGGAATCGACCCCGTATCGCGAAACACCGGCTATGGCg        >  2:136813/1‑149 (MQ=25)
                                                                                                                                           cccATTTAAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGAGCCGGCGCAGTTTGTAGAGATGCATAATTTTTTCGAAGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTATCTCATCACCGGCTAtttcgt       >  8:225626/1‑144 (MQ=38)
                                                                                                                                            ccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTg      >  6:48147/1‑149 (MQ=255)
                                                                                                                                                 ttCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAc  >  8:52289/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTAC  >  NC_002947/1392781‑1393073

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