Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,392,936 T→C 17.4% intergenic (+106/‑80) 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,9360TC17.4% 47.9 / 4.5 23intergenic (+106/‑80)PP_1214/ruvCtranscriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base T (12/7);  new base C (2/2);  total (14/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

CAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGC  >  NC_002947/1392788‑1393083
                                                                                                                                                     |                                                                                                                                                   
cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggt                                                                                                                                                     >  1:238820/1‑149 (MQ=255)
cAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGAGGGTAAAAACGCTCATACAGGTATTACGCCCATTTAAAATgtggt                                                                                                                                                     >  8:61696/1‑149 (MQ=255)
 aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggtg                                                                                                                                                    <  2:238820/149‑1 (MQ=255)
 aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggtg                                                                                                                                                    <  2:361782/149‑1 (MQ=255)
 aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATgtggt                                                                                                                                                     >  1:331183/1‑148 (MQ=255)
         gATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAc                                                                                                                                              >  8:95261/1‑147 (MQ=255)
         gATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTAc                                                                                                                                              <  7:95261/147‑1 (MQ=255)
         gATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCt                                                                                                                                            <  6:179338/149‑1 (MQ=255)
             tCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                        >  2:253971/1‑149 (MQ=255)
              ccGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                        <  5:251833/148‑1 (MQ=255)
              ccGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                                        <  2:331183/148‑1 (MQ=255)
               cgacgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGCGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGa                                                                                                                                      <  1:55738/149‑1 (MQ=18)
                 acgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGc                                                                                                                                    <  7:99775/149‑1 (MQ=18)
                    agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                                                 <  1:253971/149‑1 (MQ=255)
                                                                                tAGCGCAGGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCTTTTCAAATGTGGTCCAGTCCCTGTAGGAGCGGGTTAACCCACGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTa                                                                     <  2:277803/149‑1 (MQ=37)
                                                                                  ccgcgTGTGCCGGCCCTTGCGAGTGTAAACCCTCTCCTACTGGTTTTACGCCCATTTAAATTGTGGTGCAGTACCTGTAGGAGCGGTTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                                    <  7:61696/147‑1 (MQ=37)
                                                                                                    tCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                                  >  7:174910/1‑84 (MQ=255)
                                                                                                            aaaCCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTTCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACACCgg                                         >  5:31804/1‑149 (MQ=255)
                                                                                                              aCCCGCTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCTCTTTGTAGAGATGCATCATTTTTTCGGCGGGCTTGAAGACTCTGATTCTAGGTATCGACCCCGGCt                                       >  7:247965/1‑149 (MQ=14)
                                                                                                                ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                                     >  3:82300/1‑149 (MQ=255)
                                                                                                                ccGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                                     >  4:80437/1‑149 (MQ=255)
                                                                                                                   cTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCAGCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgcg                                   >  4:272820/1‑148 (MQ=255)
                                                                                                                   cTCCTACAGGTACTGCACCACATTTGAAATGGGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgcg                                   >  7:157005/1‑148 (MQ=18)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGAccc                                            >  6:184443/1‑123 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCg                  >  8:48915/1‑149 (MQ=255)
                                                                                                                                   gCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCg                  >  4:84328/1‑149 (MQ=255)
                                                                                                                                     ccATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACACCGGATCGCGCATCACCGGCTATGGCGTg                >  3:69631/1‑149 (MQ=255)
                                                                                                                                       cATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTgg               >  6:301789/1‑149 (MQ=255)
                                                                                                                                                    gtgCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGc  <  5:301789/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCCCATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGGCGTGGTACGCCAGACCGC  >  NC_002947/1392788‑1393083

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