Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,392,933 Δ1 bp 12.9% intergenic (+103/‑83) 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,9330T.12.9% 83.5 / 8.2 31intergenic (+103/‑83)PP_1214/ruvCtranscriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base T (16/11);  new base . (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.85e-01

TCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGC  >  NC_002947/1392786‑1393075
                                                                                                                                                    |                                                                                                                                              
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                                                                                                                            >  6:106360/1‑139 (MQ=17)
                             ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGc                                                                                                                            <  5:106360/139‑1 (MQ=17)
                                    cTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                            >  1:62796/1‑132 (MQ=255)
                                    cTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                            <  2:62796/132‑1 (MQ=255)
                                                  tttGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                          >  1:23463/1‑120 (MQ=255)
                                                  tttGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCgg                                                                                                                          <  2:23463/120‑1 (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                                                 <  6:328329/110‑1 (MQ=255)
                                                                                                                                     gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTc                                                 >  5:328329/1‑110 (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)
                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGTGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGc  <  7:225626/148‑1 (MQ=255)
                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACGCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGc  >  6:169290/1‑148 (MQ=255)
                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGc  <  6:160711/148‑1 (MQ=255)
                                                                                                                                               aaaTGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGc  >  6:96503/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGC  >  NC_002947/1392786‑1393075

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: