Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,703,248	A→G	20.3%	intergenic (+97/‑70)	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,248	0	A	G	20.3%	63.1 / 12.0	35	intergenic (+97/‑70)	lysS/PP_1497	lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (5/2);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGC  >  NC_002947/1703100‑1703395                                                                                                                                                     |                                                                                                                                                    gCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACaaa                                                                                                                                                     <  8:194612/149‑1 (MQ=255)    gTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGc                                                                                                                                                  >  7:198050/1‑149 (MQ=255)     tCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTAC‑AAGGCCt                                                                                                                                                <  5:165386/149‑1 (MQ=255)       gATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCtt                                                                                                                                               <  2:208910/149‑1 (MQ=255)        aTCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTg                                                                                                                                              <  8:197226/149‑1 (MQ=255)         tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGt                                                                                                                                             <  7:34701/149‑1 (MQ=255)         tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGt                                                                                                                                             >  3:120997/1‑149 (MQ=255)                  tCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCggg                                                                                                                                    <  4:13074/149‑1 (MQ=255)                    aTCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGt                                                                                                                                   <  2:66086/148‑1 (MQ=255)                     tCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGttt                                                                                                                                 <  6:114238/149‑1 (MQ=255)                            cccGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTa                                                                                                                                            <  7:282944/131‑1 (MQ=255)                            cccGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTa                                                                                                                                            >  8:282944/1‑131 (MQ=255)                             ccGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGa                                                                                                                         <  8:198050/149‑1 (MQ=255)                                                            gAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCGTGTAGGAGCGGGTTTTACCGCGAAGGGGCCAGGACTGACGAAACAACCCCag                                                                                            >  1:294271/1‑146 (MQ=37)                                                               cGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGAAGGAGCGGTTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGt                                                                                       <  5:109683/149‑1 (MQ=255)                                                                gAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAACGGGCCAGCACTGACGAACCACCCCCATGAGGTa                                                                                      >  5:265288/1‑149 (MQ=255)                                                                gAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTACTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTa                                                                                      >  4:60551/1‑149 (MQ=255)                                                                  gccgccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATc                                                                                    >  2:237371/1‑149 (MQ=255)                                                                   ccggcTTGCGCGGGGGCTTTTTCTTGCTTGAAGCTTATTGTGGGCAGTACCGGCCTCTTCGCGGGTAAACCCCCTCCTCCAACGGCGTGGTCGGTGCGGTTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATcc                                                                                   <  6:265288/149‑1 (MQ=25)                                                                    cgccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCCCACCAAAAAACGGCTTGGTGGAGCGGGTTTAACCCCGAAggg                                                                                                                     >  5:2196/1‑114 (MQ=37)                                                                        tttCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCCCTCCTACAAAGGCCGTGTAGGAGCGGGTTTAACCGCGAAGGGGCCAGCACAGACGAAACACCCCCCAGAGGTATccccc                                                                                >  3:135116/1‑147 (MQ=37)                                                                                 ggcTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTTAAACCGCTCCTAAAAAGGCCGTGTTGGAGCGGGTTTTACCGCGAAGGGGGCAGCACTGACGAAACACCCCCATGAGGt                                                                                       >  1:145039/1‑131 (MQ=38)                                                                                          ttGCCTGACGCTTTATGTTGTCAGTACCGCCCTCTTCCCGGGTAAACCCGCTCCTCCACGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGccggtac                                                                                                             <  2:246127/100‑6 (MQ=14)                                                                                          ttGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAATCCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGccggtac                                                                                                             >  1:246127/1‑95 (MQ=14)                                                                                                          gttgtCAGTACCGGCCTCTTCGCGGGTAAAAACCCGCCGACAAAGGGCGTGTTGGGGCGGGGTTAACCGCGAAGGGGGGCGGCCGGACGCAACCCCCCCCAGAGGGATGCCCCGGGGCGCCCGCCCGGGCGGACCCCCGCgccgcccgg                                            >  4:178798/1‑146 (MQ=255)                                                                                                                        ccTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCCTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGt                              >  6:160100/1‑149 (MQ=255)                                                                                                                                cgGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGaaa                      <  1:237371/149‑1 (MQ=255)                                                                                                                                 gggTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCAAGCCTGCCCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGCAAg                     >  8:249593/1‑149 (MQ=255)                                                                                                                                  ggTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGaaa                      >  5:33389/1‑147 (MQ=255)                                                                                                                                   gTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCg                   >  3:169144/1‑149 (MQ=255)                                                                                                                                        cccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCt                                                              >  7:73374/1‑101 (MQ=255)                                                                                                                                        cccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCt                                                              <  8:73374/101‑1 (MQ=255)                                                                                                                                          cGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACACGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa            >  4:42809/1‑149 (MQ=255)                                                                                                                                             tCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc          >  8:71271/1‑148 (MQ=255)                                                                                                                                             tCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc          >  6:172321/1‑148 (MQ=255)                                                                                                                                               cTACAAGGCCTTTGTAGGAGCGGGTTTACCCGGGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg       >  5:124949/1‑149 (MQ=255)                                                                                                                                                     aGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg       <  7:312099/143‑1 (MQ=255)                                                                                                                                                     aGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg       >  8:312099/1‑143 (MQ=255)                                                                                                                                                     aGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGc  >  7:207755/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCGTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGC  >  NC_002947/1703100‑1703395

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