Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	524,352	G→C	21.2%	intergenic (+216/‑600)	tyrS → / → PP_16SC	tyrosine‑‑tRNA ligase/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	524,352	0	G	C	21.2%	62.7 / 14.1	33	intergenic (+216/‑600)	tyrS/PP_16SC	tyrosine‑‑tRNA ligase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base C (2/5);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

TTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCA  >  NC_002947/524203‑524497                                                                                                                                                      |                                                                                                                                                    tttGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCcgcg                                                                                                                                                     >  2:254402/1‑150 (MQ=255)     gtggtTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCgcggg                                                                                                                                                 <  2:114953/150‑1 (MQ=255)             tgttgtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                         >  1:161674/1‑150 (MQ=255)                 gtAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCCCGCGCGGCGCTCGGtt                                                                                                                                     <  2:764549/150‑2 (MQ=255)                      gggagggggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATa                                                                                                                                <  2:95589/150‑1 (MQ=255)                            gggCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGTCATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCt                                                                                                                          >  1:529659/1‑150 (MQ=255)                              gCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGa                                                                                                                        >  1:294884/1‑150 (MQ=255)                                tGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAGACCTAGCGCCGCCCGCGCGTCGCTCGGTTTCATAGGCGGTGaat                                                                                                                      <  1:217743/150‑2 (MQ=21)                                  ttttGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACt                                                                                                                    <  2:593797/150‑1 (MQ=255)                                    ttGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCAGCGCCTATGAGACCGAGCGCCGCCCGCGCGGCGCTCGGTTTCATAGGCGGTGaatgtct                                                                                                                  <  1:256792/150‑6 (MQ=25)                                     tGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTc                                                                                                                 <  2:294884/150‑1 (MQ=255)                                          gggTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGc                                                                                                            <  2:131685/150‑1 (MQ=255)                                                               tCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGCAACTGTAGTGGCGAACACTAGGTAGC‑‑CTCGAtt                                                                                     <  2:446943/150‑2 (MQ=255)                                                               tCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATa                                                                                     >  2:298554/1‑150 (MQ=255)                                                                  gCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACgg                                                                                  >  2:99529/1‑150 (MQ=255)                                                                  gCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCCCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACgg                                                                                  <  1:99529/150‑1 (MQ=255)                                                                       aaCCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCttt                                                                             >  1:64158/1‑150 (MQ=255)                                                                          ccTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTcc                                                                          <  2:161674/150‑1 (MQ=255)                                                                               cAGCGCCTGTGAGATCGAGCGCCGCCCGCGTGCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTGTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTta                                                                     >  1:352748/1‑150 (MQ=255)                                                                                              cGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCttt                                                      >  2:564526/1‑149 (MQ=255)                                                                                                      gcccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATCGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTaaa                                              <  2:64158/150‑1 (MQ=255)                                                                                                      gcccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGaaa                                                                                                                      >  1:168133/1‑80 (MQ=255)                                                                                                      gcccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGaaa                                                                                                                      <  2:168133/80‑1 (MQ=255)                                                                                                        ccgcGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAgg                                            >  1:25438/1‑150 (MQ=255)                                                                                                               cGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGTCGCTCGGTCTCATAGGCGCTGAACCTTTCTTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGc                                     <  1:254402/150‑1 (MQ=255)                                                                                                                          aTTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTATCTTCACCTTATATTTCACTTTCTTGAAATTAAAGGTGGACGCCCGCTCGAATc                          <  2:529659/150‑1 (MQ=17)                                                                                                                                      tGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGCAGCATGATAC‑‑ACGTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATTCCCTTATAATgc              >  2:27037/1‑150 (MQ=14)                                                                                                                                          aCCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGCAACTGTCGTGGCGAACACTTGGCAGC‑‑CCCGATACGGTCTTTTCCCCTTATATTTCACt                                                            >  2:714221/1‑100 (MQ=11)                                                                                                                                          aCCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCACTGCAACTGTCGTGGCGAACACTTGGCAGC‑‑CCCGATACGGTCTTTTCCCCTTATATTTCACt                                                            <  1:714221/100‑1 (MQ=11)                                                                                                                                            cGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGc                                                                                                            >  2:270600/1‑52 (MQ=255)                                                                                                                                            cGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGc                                                                                                            <  1:270600/52‑1 (MQ=255)                                                                                                                                            cGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATccc                        <  1:475365/134‑1 (MQ=21)                                                                                                                                            cGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATccc                        >  2:475365/1‑134 (MQ=255)                                                                                                                                            cGAGCGCCGCCCGCGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCAGTATAATGCTCCCCa        >  1:360239/1‑150 (MQ=11)                                                                                                                                                  ccgcGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCa  >  1:544441/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TTTGGTGGTTGCTGTTGTAGGGGGAGGGGGCTGTTTTGAGAGGGTGTTGGCGGTGAGATGTGTCCGCCACAACCCTTGCAGCGCCTGTGAGATCGAGCGCCGCCCGCGTCCGGCTTGAGACATTCACCGCCTATGAAACCGAGCGCCGCGCGGGCGGCGCTCGGTCTCATAGGCGCTGAAACTGTCGTGGCGAACACTTGGCAGC‑‑CTCGATACGGTCTTTTCCCCTTATATTTCACTTTCTTGAAATTAAAGGTTGACGCCCTCTCGAATCCCCTTATAATGCGCCCCACTTCCA  >  NC_002947/524203‑524497

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