Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,702,766	A→G	8.5%	E373G (GAG→GGG)	lysS →	lysine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,702,766	0	A	G	8.5%	131.1 / 3.3	47	E373G (GAG→GGG)	lysS	lysine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base G (2/2);  total (24/23)
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.06e-01

TCGACTCGATCCTCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCGCCAACGCC  >  NC_002947/1702622‑1702911                                                                                                                                                 |                                                                                                                                                  tCGACTCGATCCTCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTg                                                                                                                                               >  6:286004/1‑149 (MQ=255)   gACTCGATCCTCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTgg                                                                                                                                              <  2:153679/148‑1 (MQ=255)   gACTCGATCCTCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGt                                                                                                                                             >  8:13089/1‑149 (MQ=255)             tCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAg                                                                                                                                   >  5:226143/1‑149 (MQ=255)             tCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAg                                                                                                                                   >  7:38423/1‑149 (MQ=255)                    cAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGaa                                                                                                                             <  6:226143/148‑1 (MQ=255)                       cccGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGc                                                                                                                         <  8:38423/149‑1 (MQ=255)                       cccGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGc                                                                                                                         <  1:8477/149‑1 (MQ=255)                               tGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACtt                                                                                                                  >  2:145400/1‑148 (MQ=255)                               tGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACtt                                                                                                                  >  6:291424/1‑148 (MQ=255)                                 aCCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGc                                                                                                                         >  5:191435/1‑139 (MQ=255)                                 aCCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGc                                                                                                                         <  6:191435/139‑1 (MQ=255)                                     ccgccgACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTAcc                                                                                                           >  5:197527/1‑149 (MQ=255)                                       gccgACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGa                                                                                                         <  4:294196/149‑1 (MQ=255)                                       gccgACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGa                                                                                                         <  1:145400/149‑1 (MQ=255)                                        ccgACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGTTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAg                                                                                                        <  4:99067/149‑1 (MQ=255)                                           aCCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTAc                                                                                                     <  5:320977/149‑1 (MQ=255)                                               gCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGGGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGt                                                                                                 >  6:196237/1‑149 (MQ=255)                                               gCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTAGAAGGGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTAACGAGTAccggg                                                                                                 >  8:192620/1‑146 (MQ=255)                                                   taggTCGACCGCGCCCGCGAGATTGCCAGGATGGGCGGCGCCAAGGTGGTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTgg                                                                                                                                              <  6:263671/97‑1 (MQ=255)                                                   gACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTgg                                                                                                                                              >  5:263671/1‑100 (MQ=255)                                                    aCGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGaa                                                                                            <  3:72544/149‑1 (MQ=255)                                                                 ccGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTgg                                                                               >  4:145515/1‑149 (MQ=255)                                                                   gcgAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTc                                                                                               >  7:210255/1‑131 (MQ=255)                                                                   gcgAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTc                                                                                               <  8:210255/131‑1 (MQ=255)                                                                            ccAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTGCCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCaa                                                                     >  5:277461/1‑148 (MQ=255)                                                                            ccAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCaa                                                                     <  6:277461/148‑1 (MQ=255)                                                                            ccAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAac                                                                    >  3:246775/1‑149 (MQ=255)                                                                              aGGAAGGCCGGCTCCAAGGGGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGGGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAacga                                                                  <  7:192620/149‑1 (MQ=255)                                                                              aGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGGGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAacga                                                                  <  5:196237/149‑1 (MQ=255)                                                                               ggAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAacgac                                                                 <  7:13089/149‑1 (MQ=255)                                                                                     ccGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAAc                                                           <  5:291424/149‑1 (MQ=255)                                                                                     ccGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAAc                                                           <  4:246775/149‑1 (MQ=255)                                                                                                    tGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTaccgac                                            >  1:109332/1‑149 (MQ=255)                                                                                                        ccATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGtt                                        <  2:109332/149‑1 (MQ=255)                                                                                                             aaGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAg                                   <  3:145515/149‑1 (MQ=255)                                                                                                                       cAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATc                          >  8:205305/1‑148 (MQ=255)                                                                                                                       cAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATc                          <  7:205305/148‑1 (MQ=255)                                                                                                                       cAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCg                         >  2:313135/1‑149 (MQ=255)                                                                                                                             gCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGcc                   >  8:55035/1‑149 (MQ=255)                                                                                                                                     tgatTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATc           <  5:286004/149‑1 (MQ=255)                                                                                                                                      gatTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCg          >  3:267327/1‑149 (MQ=255)                                                                                                                                      gatTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCg          >  5:98424/1‑149 (MQ=255)                                                                                                                                             gAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCGCCAACGc   >  1:221487/1‑149 (MQ=255)                                                                                                                                              aaGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCGc         >  2:66629/1‑142 (MQ=255)                                                                                                                                              aaGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCGc         <  1:66629/142‑1 (MQ=255)                                                                                                                                              aaGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCGCCAACGcc  >  6:134679/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TCGACTCGATCCTCAAGTACAACCCGGAGCTGACCGCCGCCGACCTGCAGGACGTCGACCGCGCCCGCGAGATTGCCAGGAAGGCCGGCGCCAAGGTGCTGGGCCATGAAGGCCTGGGCAAGCTGCAGGTGATGATTTTCGAAGAGCTGGTCGAGCACAAGCTGGAACAGCCGCACTTCATTACCGAGTACCCGTTCGAAGTGTCGCCGTTGGCCCGTCGCAACGACGACAACCCGGCCGTTACCGACCGTTTCGAGCTGTTCATCGGTGGCCGCGAAATCGCCAACGCC  >  NC_002947/1702622‑1702911

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