Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,498,379	C→A	13.7%	A100D (GCC→GAC)	trpS →	tryptophan‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,498,379	0	C	A	13.7%	78.5 / 2.9	30	A100D (GCC→GAC)	trpS	tryptophan‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base A (2/2);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACCCGCTGCGCATCCAGCGCTCGCGCCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGATATCCTGATGTTCAAC  >  NC_002947/1498238‑1498524                                                                                                                                              |                                                                                                                                                  acCCGCTGCGCATCCAGCGCTCGCGCCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAg                                                                                                                                           <  2:46941/150‑1 (MQ=255)         cgcATCCAGCGCTCGCGCCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCctgctg                                                                                                                                   <  2:635989/150‑1 (MQ=255)              ccAGCGCTCGCGCCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGaa                                                                                                                                 >  2:33447/1‑147 (MQ=255)              ccAGCGCTCGCGCCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGaa                                                                                                                                 <  1:33447/147‑1 (MQ=255)                       gcgcCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGACGCCAAAGGCCTGCTGAACCGCGCGCATGc                                                                                                                     >  2:553172/1‑150 (MQ=255)                       gcgcCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGACGCCAAAGGCCTGCTGAACCGCGCGCATGc                                                                                                                     >  2:553138/1‑150 (MQ=255)                           cTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCgcgcgc                                                                                                                         >  2:402414/1‑142 (MQ=255)                           cTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCgcgcgc                                                                                                                         <  1:402414/142‑1 (MQ=255)                             ggAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACaa                                                                                                               <  2:13513/150‑1 (MQ=255)                              gAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAg                                                                                                              <  2:575912/150‑1 (MQ=255)                                aaTCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCCGTCCGACATCCCGGAAATCCTCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGc                                                                                                            >  1:90581/1‑150 (MQ=255)                                      cgccACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGt                                                                                                      <  1:674586/150‑1 (MQ=255)                                       gccACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTg                                                                                                     <  2:471458/150‑1 (MQ=255)                                         cACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAg                                                                                                                                           <  2:193281/110‑1 (MQ=255)                                         cACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAg                                                                                                                                           >  1:193281/1‑110 (MQ=255)                                                     cGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGACGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCga                                                                                       <  1:553138/150‑1 (MQ=255)                                                     cGGTGGCCTGGACCCGGACAAGGTGACCCTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGACGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCga                                                                                       <  1:553172/150‑1 (MQ=255)                                                                  ccGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGAc                                                                          >  1:534446/1‑150 (MQ=255)                                                                    ggACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGAccc                                                                        >  2:232140/1‑150 (MQ=255)                                                                                 ttCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTg                                                           >  2:226929/1‑150 (MQ=255)                                                                                                   aTCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGc                                         <  1:232140/150‑1 (MQ=255)                                                                                                                ccGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGa                            >  2:167112/1‑150 (MQ=255)                                                                                                                  gAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATg                          >  2:589422/1‑150 (MQ=255)                                                                                                                          cTGGCTGCTGCCCTACATCGGCGCCAAAGGCCTGCTGACCCGCGAGTATGCCTCCAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGCCCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGa                  <  1:115410/150‑1 (MQ=255)                                                                                                                          cTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGa                  <  1:351526/150‑1 (MQ=255)                                                                                                                            ggCTGCTGACCTGCGCGGCCGCCAAAGGCCTGCTGACCCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGata                <  2:287728/150‑1 (MQ=255)                                                                                                                               tgctgACCTTCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGt                                              >  1:734334/1‑117 (MQ=255)                                                                                                                               tgctgACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGt                                              <  2:734334/117‑1 (MQ=255)                                                                                                                                      cTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGATATCCTGATGtt      >  1:472397/1‑150 (MQ=255)                                                                                                                                          gTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGATATCCTGATGTTCAAc  <  1:226929/150‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                  ACCCGCTGCGCATCCAGCGCTCGCGCCTGGAAATCGCCGCCACCTGGCTGGCCGGTGGCCTGGACCCGGACAAGGTGACCTTCTATCGCCAGTCCGACATCCCGGAAATCCCCGAGCTGACCTGGCTGCTGACCTGCGTTGCCGCCAAAGGCCTGCTGAACCGCGCGCATGCCTACAAGGCCTCGGTGGACAAGAACGTCGAGAACGGCGAAGACCCGGACGCCGGCGTGACGATGGGCCTGTTCAGCTACCCGGTGCTGATGGCCGCCGATATCCTGATGTTCAAC  >  NC_002947/1498238‑1498524

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