Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,390,996	G→C	8.5%	G268A (GGC→GCC)	aspS →	aspartate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,390,996	0	G	C	8.5%	139.4 / 2.5	46	G268A (GGC→GCC)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (19/23);  new base C (2/2);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAGGTGTTCGCCGGCCCGGC  >  NC_002947/1390855‑1391134                                                                                                                                                |                                                                                                                                           aCGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGTTCGCGGTGTAGACCAGGAGAGCAAGATCCGCAAGCTGTTCAACGAATTGGTCGATCTGGAAATCGGCCAAt                                                                                                                                       >  1:115835/1‑149 (MQ=255)   gAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTccc                                                                                                                                   <  5:106110/149‑1 (MQ=255)    aaGACCTGCGTGCCGACCTTCAGCCGGAATTCACCCAGATCGAAATCGAGACCTGTTTCCTCTATGAAAGCTAGATCATGGGCCTGACC‑‑GAAAGCATTATCCGCAAGCTTTTCATGGAAGTGCTCGATCTGGAATTCGCCGAATTccag                                                                                                                                  <  6:350809/149‑3 (MQ=255)            cGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGa                                                                                                                           <  1:289491/148‑1 (MQ=255)               gCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCtt                                                                                                                       >  2:333216/1‑149 (MQ=255)                  gACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCg                                                                                                                     >  4:304735/1‑148 (MQ=255)                      gTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAagag                                                                                                                <  7:203528/149‑1 (MQ=255)                       gcAGCCGGAATTCACCCAGATCTACATCGTGACCAGTTTGCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GACAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTGGGCGAATTCCCGCACATGACCTTCGAAGAgg                                                                                                               <  8:307336/148‑1 (MQ=255)                         aggcGGAATTCAGCCAGTTCGACATCGATACCAGTTTCCTCTATGAATGGGATATCATGTGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCTAATTCCCGGACATGACCTTCGAAGAGGcc                                                                                                             <  4:217504/146‑1 (MQ=255)                         aGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTccc                                                                                                                                   >  4:314239/1‑127 (MQ=255)                         aGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTccc                                                                                                                                   <  3:314239/127‑1 (MQ=255)                         aGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGcc                                                                                                             <  8:82622/149‑1 (MQ=255)                         aGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAACCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGc                                                                                                              >  3:9412/1‑148 (MQ=255)                             ggAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATgc                                                                                                         <  1:333216/149‑1 (MQ=255)                             ggAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATgc                                                                                                         <  4:272403/149‑1 (MQ=255)                                     cccAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTGCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTAc                                                                                                 <  3:290237/149‑1 (MQ=255)                                              aCATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCgg                                                                                                                                            >  5:316715/1‑97 (MQ=255)                                              aCATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCgg                                                                                                                                            <  6:316715/97‑1 (MQ=255)                                                 tCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATTCGCCGCTACGGCTCCGACAAg                                                                                     >  5:15368/1‑149 (MQ=255)                                                 tCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAg                                                                                     <  3:304735/149‑1 (MQ=255)                                                       ccAGTTTCCTAGATGAAAGCGAGATCATGGGCCTGAAC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGAc                                                                               <  4:9412/149‑1 (MQ=255)                                                        cAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGc                                                                                    >  7:240160/1‑143 (MQ=255)                                                        cAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGc                                                                                    <  8:240160/143‑1 (MQ=255)                                                                     gAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCt                                                                 >  8:75335/1‑149 (MQ=255)                                                                      aaaGCGAGATCATGGGCCTGACC‑‑GAAACCATGATCCGCTAGCTGTGCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACTGCTCCGACAAGCCAGACCTGCGTATTCCGCt                                                                 <  2:115835/148‑1 (MQ=255)                                                                         gCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGaa                                                             >  1:192715/1‑149 (MQ=255)                                                                           gagaTCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGCCGAATTcc                                                                                                                                    <  8:187254/76‑1 (MQ=255)                                                                           gagaTCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGCCGAATTcc                                                                                                                                    >  7:187254/1‑76 (MQ=255)                                                                             gaTCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACAt                                                                                                                             >  4:313877/1‑81 (MQ=255)                                                                             gaTCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACAt                                                                                                                             <  3:313877/81‑1 (MQ=255)                                                                             gaTCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTgg                                                         <  2:192715/149‑1 (MQ=255)                                                                              agcaggggCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGGGCGCGCTCTGGAATTCGGCGATTTCCCGCAGATGAGCTTCGAGGGGGCCAGGCGCCGCTACGGCTCCGACACGCCAGACCTGCGGAGTGAGCTGGCACTGGt                                                        >  4:69000/6‑149 (MQ=255)                                                                                       cTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTg                                                                <  4:45607/132‑1 (MQ=255)                                                                                       cTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTg                                                                >  3:45607/1‑132 (MQ=255)                                                                                            c‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGca                                                                                                                                >  5:148800/1‑63 (MQ=255)                                                                                            c‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGca                                                                                                                                <  6:148800/63‑1 (MQ=255)                                                                                                         gccGCAAGCGGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGAGAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGAt                               <  4:96268/148‑1 (MQ=255)                                                                                                            gCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGc                                                                                    >  2:357071/1‑93 (MQ=255)                                                                                                            gCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGc                                                                                    <  1:357071/93‑1 (MQ=255)                                                                                                            gCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTc                            >  6:242010/1‑149 (MQ=255)                                                                                                                      cAAGGAAGTGCTCGATCTGGAATTCGGCTAATTCCCTCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGATCCGACAAGCCAGACCTGAGTATTCCGCTTTAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAg                   <  5:199773/148‑1 (MQ=255)                                                                                                                       aaGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAGgt                 >  5:6329/1‑149 (MQ=255)                                                                                                                          gAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCg                           <  2:334297/136‑1 (MQ=255)                                                                                                                          gAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCg                           >  1:334297/1‑136 (MQ=255)                                                                                                                            aGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAGGTGTTCg            >  6:142163/1‑149 (MQ=255)                                                                                                                                cTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAGGTGTTCGccgg        >  2:286426/1‑149 (MQ=255)                                                                                                                                      cgggAATTCGCCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAGGTGTTCGccggcccggc  >  8:187037/3‑149 (MQ=255)                                                                                                                                        caaaTTCGGCGAATTCCCGCACATGACCTTCGACGAGGCCACGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAg                   <  4:276023/128‑1 (MQ=255)                                                                                                                                                |                                                                                                                                           ACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTGACC‑‑GAAAGCATGATCCGCAAGCTGTTCAAGGAAGTGCTCGATCTGGAATTCGGCGAATTCCCGCACATGACCTTCGAAGAGGCCATGCGCCGCTACGGCTCCGACAAGCCAGACCTGCGTATTCCGCTGGAACTGGTCGACGTTGCCGACCAGCTGAAGGATGTCGACTTCAAGGTGTTCGCCGGCCCGGC  >  NC_002947/1390855‑1391134

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