Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,703,263	T→C	11.1%	intergenic (+70/‑97)	lysS → / → PP_1497	lysine‑‑tRNA ligase/TetR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,703,263	0	T	C	11.1%	112.1 / 5.7	45	intergenic (+70/‑97)	lysS/PP_1497	lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (17/23);  new base C (3/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

TGGTGATGCTGCTGACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCA  >  NC_002947/1703120‑1703402                                                                                                                                                |                                                                                                                                             tGGTGATGCTGCTGACCAACTCGCCGTCGCTCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTcgc                                                                                                                                        <  1:543351/150‑1 (MQ=255)   gTGATGCTGCTGACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCgg                                                                                                                                      >  2:229792/1‑150 (MQ=255)          tgctgACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGTGGGTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAAcc                                                                                                                               >  1:918739/1‑150 (MQ=255)          tgctgACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAAcc                                                                                                                               <  1:229792/150‑1 (MQ=255)           gctgACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAAccc                                                                                                                              >  1:235009/1‑150 (MQ=255)            ctgACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCCCTTCGCGGGTAAACCCg                                                                                                                             >  2:495309/1‑150 (MQ=255)              gACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                           <  1:963766/150‑1 (MQ=255)                  gaCTCGCCGTCTATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGCGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTa                                                                                                                       >  2:17168/2‑150 (MQ=255)                   aCTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                      >  2:415805/1‑150 (MQ=255)                     tCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                    >  1:735052/1‑150 (MQ=255)                     tCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCCCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                    <  1:495309/150‑1 (MQ=255)                       gCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACaagg                                                                                                                  >  2:457435/1‑148 (MQ=255)                          gTCGATCCGCGACGTCATACTGTTCCCCCACATGCCTCCGCAGGCCTTAGTTAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGcc                                                                                                               <  2:735052/150‑1 (MQ=39)                                    gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGc                                                                                                    <  2:235009/150‑1 (MQ=255)                                    gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCCCTTCGCGGGTAAACCCGCTCCTACACGGCCT‑TTGTAGGAGc                                                                                                    <  2:74481/150‑1 (MQ=38)                                     aCGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCg                                                                                                   <  1:574178/150‑1 (MQ=255)                                      cGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCgg                                                                                                  <  2:788042/150‑1 (MQ=255)                                      cGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                  <  2:414684/150‑1 (MQ=255)                                       gTCATCCTGTGCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCT‑TTGTAGGAGCggg                                                                                                 <  1:17168/150‑1 (MQ=255)                                                   ccGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACACGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGaa                                                                                     >  1:525522/1‑150 (MQ=255)                                                     gCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTcgcg                                                                                                                                       >  1:154957/1‑99 (MQ=255)                                                     gCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTcgcg                                                                                                                                       <  2:154957/99‑1 (MQ=255)                                                             gTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCa                                                                           <  2:285380/150‑1 (MQ=255)                                                                cGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCC‑tttgt                                                                                                         >  1:1072551/1‑113 (MQ=255)                                                                cGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCC‑tttgt                                                                                                         <  2:1072551/116‑4 (MQ=255)                                                                  cAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGAc                                                                      <  1:215207/150‑1 (MQ=255)                                                                  cAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACACGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGAc                                                                      <  2:525522/150‑1 (MQ=255)                                                                      ccTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGaaa                                                                  <  1:415805/150‑1 (MQ=255)                                                                       cTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAac                                                                 <  2:323573/150‑1 (MQ=255)                                                                        tgagtgagCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCgg                                                                                                                                      >  1:1101777/1‑81 (MQ=255)                                                                        tgagtgagCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCgg                                                                                                                                      <  2:1101777/81‑1 (MQ=255)                                                                                 tGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTACCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGa                                                       >  1:349790/1‑150 (MQ=255)                                                                                                  gAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCcg                                                                                         <  1:721065/99‑1 (MQ=255)                                                                                                  gAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCcg                                                                                         >  2:721065/1‑99 (MQ=255)                                                                                                           ttttGTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGc                             <  2:542858/150‑1 (MQ=255)                                                                                                                          cTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGAc                                                                      >  2:411521/1‑94 (MQ=255)                                                                                                                          cTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGAc                                                                      <  1:411521/94‑1 (MQ=255)                                                                                                                           tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGaaa                                                                  <  1:97821/97‑1 (MQ=255)                                                                                                                           tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGaaa                                                                  >  2:97821/1‑97 (MQ=255)                                                                                                                           tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                  >  2:1120017/1‑65 (MQ=14)                                                                                                                           tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                  <  1:1120017/65‑1 (MQ=14)                                                                                                                              aTGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCgg          <  2:349790/150‑1 (MQ=255)                                                                                                                               tgttgtCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGc         >  1:603812/1‑150 (MQ=255)                                                                                                                                gttgtCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCa        >  2:567767/1‑150 (MQ=255)                                                                                                                                      agtgctggCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCT‑TTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCa  >  1:814350/7‑150 (MQ=18)                                                                                                                                                |                                                                                                                                             TGGTGATGCTGCTGACCAACTCGCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCA  >  NC_002947/1703120‑1703402

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