Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,703,294:1	+T	17.5%	intergenic (+101/‑66)	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,294	1	.	T	17.5%	115.9 / 26.7	40	intergenic (+101/‑66)	lysS/PP_1497	lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (18/15);  new base T (2/5);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.17e-01

GATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGA  >  NC_002947/1703148‑1703441                                                                                                                                                    |                                                                                                                                                    gATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCC‑tttg                                                                                                                                                 >  2:418736/1‑148 (MQ=38)    ccGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCT‑TTGTAg                                                                                                                                               >  1:439999/1‑150 (MQ=255)        gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGc                                                                                                                                           <  1:418736/150‑1 (MQ=37)        gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGc                                                                                                                                           <  2:353654/150‑1 (MQ=255)          cGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCgg                                                                                                                                         <  1:506136/150‑1 (MQ=255)          cGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                                                         <  2:439999/150‑1 (MQ=255)                 cTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCT‑TTGTAGGAGCgg                                                                                                                                         >  1:348008/1‑143 (MQ=255)                 cTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                                                         <  2:348008/143‑1 (MQ=255)                       ccGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                            >  1:246029/1‑150 (MQ=255)                           acaTGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAgggg                                                                                                                        <  2:115473/150‑1 (MQ=255)                                          ccTGAGTGAGCTGAACGAGCCGCCTTTCTAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAGAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGaaa                                                                                                         <  2:229940/150‑1 (MQ=255)                                                   gCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCAt                                                                                                >  1:68386/1‑150 (MQ=255)                                                     tGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGa                                                                                              >  1:315806/1‑150 (MQ=255)                                                                               ttttCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGc                                                                    <  2:50063/150‑1 (MQ=255)                                                                               ttttCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGc                                                                    <  2:246029/150‑1 (MQ=255)                                                                                    ttGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGAc                                                               <  2:304050/150‑1 (MQ=255)                                                                                    ttGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGAc                                                               <  2:68386/150‑1 (MQ=255)                                                                                        cTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                                                         <  1:16361/72‑1 (MQ=18)                                                                                        cTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCTTTGTAGGAGCgg                                                                                                                                         >  2:16361/1‑72 (MQ=18)                                                                                               tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGcgccgc                                                    >  2:381818/1‑150 (MQ=255)                                                                                                ttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑CTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCgccgcc                                                   >  2:469707/1‑150 (MQ=255)                                                                                                      tgtCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATc                                             <  2:315806/150‑1 (MQ=255)                                                                                                           gTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGGTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTACCTCCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCAc                                        >  1:338132/1‑150 (MQ=255)                                                                                                               cGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATc                                             >  2:196205/1‑141 (MQ=255)                                                                                                               cGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATc                                             <  1:196205/141‑1 (MQ=255)                                                                                                                      ttCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGc                                                                         >  2:292393/1‑106 (MQ=255)                                                                                                                      ttCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGc                                                                         <  1:292393/106‑1 (MQ=255)                                                                                                                           gggTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCa                                                                       >  2:446935/1‑103 (MQ=255)                                                                                                                           gggTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCa                                                                       <  1:446935/103‑1 (MQ=255)                                                                                                                                   ccGCTCCTACAAGGCCT‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa                >  2:16603/1‑150 (MQ=255)                                                                                                                                   ccGCTCCTACAAAGGCC‑TTGTAGGAGCGNNTTTACCCGCGAAGGNGNCNGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa                >  1:366150/1‑150 (MQ=255)                                                                                                                                   ccGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa                >  2:501881/1‑150 (MQ=255)                                                                                                                                   ccGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAACCCGCCCCCGCCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa                >  2:501264/1‑150 (MQ=255)                                                                                                                                   ccGCTCCTAC‑AAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGa                                                                <  2:355642/102‑1 (MQ=25)                                                                                                                                   ccGCTCCTAC‑AAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGa                                                                >  1:355642/1‑102 (MQ=37)                                                                                                                                     gCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATc                                             >  1:94291/1‑119 (MQ=255)                                                                                                                                     gCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATc                                             <  2:94291/119‑1 (MQ=255)                                                                                                                                     gCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc              >  1:74989/1‑150 (MQ=255)                                                                                                                                        ccTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg           <  1:50942/150‑1 (MQ=255)                                                                                                                                          tACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGc         >  1:479471/1‑150 (MQ=255)                                                                                                                                                 gCC‑TTGTAGGAGCGGGTTTACCCCCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGctcga  <  2:366150/150‑5 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCC‑TTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGA  >  NC_002947/1703148‑1703441

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