Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,392,917	A→G	8.1%	intergenic (+87/‑99)	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_002947	1,392,917	0	A	G	8.1%	131.6 / 2.9	50	intergenic (+87/‑99)	PP_1214/ruvC	transcriptional regulator/crossover junction endodeoxyribonuclease RuvC
Reads supporting (aligned to +/- strand):  ref base A (22/23);  new base G (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

ACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGG  >  NC_002947/1392769‑1393065                                                                                                                                                     |                                                                                                                                                      acgTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTa                                                                                                                                                       <  1:198420/149‑1 (MQ=255)        gAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CAt                                                                                                                                               >  1:60780/1‑149 (MQ=255)        gAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CAt                                                                                                                                               <  6:321757/149‑1 (MQ=255)        gAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CAt                                                                                                                                               >  2:190834/1‑149 (MQ=255)         aaCGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CAtt                                                                                                                                              <  4:22486/149‑1 (MQ=255)         aaCGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CAtt                                                                                                                                              <  2:358531/149‑1 (MQ=255)          aCGACTATGCCAATGCCCAGATAACCGACGAGATCATCGAAAACCTCGGCTAAGGTGAGTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CAttt                                                                                                                                             <  8:164199/149‑1 (MQ=39)             gcTATTCCAATGCCCAGATTTGCGACGAGATCATGGAAATCCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTACCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCaa                                                                                                                                          <  6:305190/148‑1 (MQ=255)                aTTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTACAACAGGTATTACGCC‑CATTTCAAAtg                                                                                                                                       >  7:196557/1‑149 (MQ=255)                     aaTGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATgtggtg                                                                                                                                  <  1:190834/149‑1 (MQ=255)                          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATgtggtg                                                                                                                                  >  6:184837/1‑144 (MQ=255)                          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATgtggtg                                                                                                                                  <  5:184837/144‑1 (MQ=255)                          ccAGATTTCCGACGAGATCATGGAAAACCTCGGCGAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTAAAGGTATTACGCC‑CATTTCAAATTTGGTGGAGt                                                                                                                              >  1:237498/1‑148 (MQ=255)                            aGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTAcc                                                                                                                           >  3:290046/1‑149 (MQ=255)                               tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGt                                                                                                                        <  6:111073/149‑1 (MQ=255)                               tttCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGt                                                                                                                        <  2:60780/149‑1 (MQ=255)                                 tCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAg                                                                                                                      >  2:181065/1‑149 (MQ=255)                                   cgacgaGAACATGGAAAACCTCTTGTAATGTGATTTGGCATGTTGTGCCAAATATGAGCGTTTGATATCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTACTGCACCACATTTGAAATG‑GGCGTAATACCTGTAGGa                                                                                                                    <  8:268483/149‑1 (MQ=9)                                     acgaGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                  <  1:181065/149‑1 (MQ=255)                                        agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                               <  7:59844/149‑1 (MQ=255)                                        agaTCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCggg                                                                                                               <  7:22202/149‑1 (MQ=255)                                            cATTGAAAAGCTCGGCTCAGGTGATTTGGCATTTTGTGCGAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGGC‑CATTTCCAATGTGGTGCAGTACCTGTAGGAGCGGGttt                                                                                                            <  1:226901/148‑1 (MQ=38)                                              tGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTGAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAcc                                                                                                         >  2:115681/1‑149 (MQ=255)                                               ggAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGGATTACGCC‑CATTTCAAATGGGGTGCAGTACCTGTAGGAGCGGGTTTAacc                                                                                                        >  6:188158/1‑149 (MQ=255)                                                gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTAccc                                                                                                        <  4:290046/148‑1 (MQ=255)                                                gAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCcg                                                                                                       <  2:300286/149‑1 (MQ=255)                                                 aaaaCCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                  <  5:140942/137‑1 (MQ=255)                                                 aaaaCCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGc                                                                                                                  >  6:140942/1‑137 (MQ=255)                                                     ccTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTC‑AATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                 <  1:111762/149‑1 (MQ=255)                                                      cTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCTACTACAGGTTTTAAGCC‑CATTTAAAATGTGGGGCAGTACCTGTAGGAGCGGGGTTACCCGCGcaac                                                                                                 >  1:283912/1‑148 (MQ=255)                                                               ggTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGgcg                                                                                         >  7:269211/1‑148 (MQ=255)                                                                          aTTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTTAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTAACGGTAGGAGCGGGTTTACCCGCGAAATGGCCGGGGCAGTTTGTtgag                                                                             >  5:311830/1‑149 (MQ=37)                                                                             tgcgcccACCTCTGCGCGTTTCACCGCGCCTGGGCCGGCCCTTGCGCGGTCAACCCCGCCCCAACCGTTTAACCCCC‑CATTTCAATTGGGGTCCCTTCCCTTTTGCCGGGGGTTAACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGAt                                                                           <  8:196557/143‑1 (MQ=255)                                                                                         tGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAAtgtg                                                                                                                                     >  3:358069/1‑78 (MQ=255)                                                                                         tGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAAtgtg                                                                                                                                     <  4:358069/78‑1 (MQ=255)                                                                                          gAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTAAGCC‑CCTTTAAAATGTGGTGGAATAACTGTAGGAGCGGGTTTAACCGCGAAAGGGCCGGCGCACTTTGTTTAAATGCATAAgtttttg                                                              >  2:61826/1‑147 (MQ=21)                                                                                                      gcgcCTGTGCCGGCCCTTTCGCGGTTAAACCCGCTCCTACCGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGGGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTAt                                                  <  5:188158/148‑1 (MQ=255)                                                                                                           tgtgCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGt                                                                                                                        <  6:176450/73‑1 (MQ=255)                                                                                                           tgtgCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGt                                                                                                                        >  5:176450/1‑73 (MQ=255)                                                                                                                  gCCGTTTGGCGGGTAAACCCGCTCCTACAGTTATTACGCC‑ATTTTCAAATTTGTTGATGAACCTGTAGGTGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTa                                     <  1:61826/149‑1 (MQ=37)                                                                                                                   ccgtttCGCGGGAAAACCCGCTCATACAGTTATTACGCC‑CATTCAAATTGTGGTGCAGTACCTTTCGTAGGGGTTTTACCCGCTAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAg                                    <  2:237498/146‑1 (MQ=255)                                                                                                                                    ccGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcg                   >  4:317374/1‑149 (MQ=255)                                                                                                                                     cGCTCCTACAGGTACTGCACCACATTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTcgc                  >  3:156649/1‑149 (MQ=18)                                                                                                                                          cTACAGCTATTACGCC‑CATTTCAAATGTGCTGCAGTACCTGTAGGAGCGCGTTTACCATCGAAAGGGACGGCGCAGTTTGTAGAGATGAATCAATTTTTAGACGGGATTTATTACTCTGATTCTACATATCGCCCCCGGATCGCTCAt              >  5:122211/1‑148 (MQ=255)                                                                                                                                           tACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCa            >  7:201251/1‑149 (MQ=255)                                                                                                                                           tACAGGTACTGCACCACATTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGActct                                            <  6:266350/117‑1 (MQ=14)                                                                                                                                           tACAGGTACTGCACCACATTTGAAATG‑GGCGTAATACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGActct                                            >  5:266350/1‑117 (MQ=14)                                                                                                                                             cAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATc                               <  5:212381/128‑1 (MQ=255)                                                                                                                                             cAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATc                               >  6:212381/1‑128 (MQ=255)                                                                                                                                             cAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCAcc          >  8:105444/1‑149 (MQ=255)                                                                                                                                                  aTTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCTACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGcaa     >  7:230206/1‑147 (MQ=255)                                                                                                                                                     aCGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATgg  >  2:154062/1‑149 (MQ=255)                                                                                                                                                     aCGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATg   <  8:201251/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACGTGCAGAACGTCTATTCCAATGCCCAGATTTCCGACGAGATCATGGAAAACCTCGGCTAAGGTGATTTGGCATTTTGTGCCAACTCTGAGCGTTTGATAGCGCCTGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGTATTACGCC‑CATTTCAAATGTGGTGCAGTACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCGCAGTTTGTAGAGATGCATCATTTTTTCGACGGGCTTTATGACTCTGATTCTAGGTATCGACCCCGGCTCGCGCATCACCGGCTATGG  >  NC_002947/1392769‑1393065

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: ‑