Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,972,594	G→A	13.3%	T427M (ACG→ATG)	PP_2600 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,972,594	0	G	A	13.3%	72.2 / 3.8	30	T427M (ACG→ATG)	PP_2600	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (11/15);  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 = 8.49e-01

TGTGGGTCACCCTGCTTGCGGTTGAACGGCGTGGAACCGCGGCCAAACGCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCTCCTGGCCCTGCGGCACTG  >  NC_002947/2972447‑2972735                                                                                                                                                    |                                                                                                                                              tgtgGGTCACCCTGCTTGCGGTTGAACGGCGTGGAACCGCGGCCAAACGCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCg                                                                                                                                               >  1:163814/1‑148 (MQ=255)                        gAACGGCGTGGAACCGCGGCCAAACGCCGGATCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGt                                                                                                                       >  1:160855/1‑149 (MQ=255)                        gAACGGCGTGGAACAGCGGCCAAACGCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCg                                                                                                                        <  6:227557/148‑1 (MQ=255)                                      cgcgGCCAAACGCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCATAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGc                                                                                                         >  6:192862/1‑149 (MQ=255)                                             aaaCGCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCATAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCg                                                                                                  >  3:154831/1‑149 (MQ=255)                                                 gCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAg                                                                                              <  6:147854/149‑1 (MQ=255)                                                              cGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGt                                                                                                                             <  8:59232/105‑1 (MQ=255)                                                              cGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCcgggcgg                                                                                 <  2:197289/149‑1 (MQ=255)                                                                         tGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGccc                                                                                                       <  6:239895/116‑1 (MQ=255)                                                                         tGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGccc                                                                                                       >  5:239895/1‑116 (MQ=255)                                                                           ggTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAgg                                                                     <  2:163814/148‑1 (MQ=255)                                                                           ggTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCATAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGc                                                                    <  5:192862/149‑1 (MQ=255)                                                                                   aGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAgcgg                                                            <  1:76730/149‑1 (MQ=255)                                                                                    gTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAgcggc                                                           >  8:206064/1‑149 (MQ=255)                                                                                    gTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCATAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAgcggc                                                           <  4:154831/149‑1 (MQ=255)                                                                                       aGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGgcgc                                                         >  4:257502/1‑148 (MQ=255)                                                                                       aGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGgcgc                                                         <  2:138150/148‑1 (MQ=255)                                                                                               tGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCga                                                <  7:206064/149‑1 (MQ=255)                                                                                                  cgcgcAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCcg                                                                                      >  3:62195/1‑108 (MQ=255)                                                                                                  cgcgcAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCcg                                                                                      <  4:62195/108‑1 (MQ=255)                                                                                                   gcgcAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCgatgat                                            >  7:135101/1‑149 (MQ=255)                                                                                                                  gTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATc                                           <  5:203416/135‑1 (MQ=255)                                                                                                                  gTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATc                                           >  6:203416/1‑135 (MQ=255)                                                                                                                     aaTGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAg                           <  3:257502/148‑1 (MQ=255)                                                                                                                                 cTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCcgggcgg                                                                                 <  4:82902/82‑1 (MQ=255)                                                                                                                                 cTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCcgggcgg                                                                                 >  3:82902/1‑82 (MQ=255)                                                                                                                                     gCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCTCCTGGCCCTg          <  2:160855/149‑1 (MQ=255)                                                                                                                                      cGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCTCCTGGCCCTGc         <  3:134276/149‑1 (MQ=255)                                                                                                                                           tgctCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCTCCTGGCCCTGCGGCAc    >  8:156724/1‑149 (MQ=255)                                                                                                                                             ctCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCTCCTGGCCCTGCGGCACTg  >  3:105402/1‑149 (MQ=255)                                                                                                                                              tCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCACCAGGCCCTCCGGCACTg  >  3:4626/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                              TGTGGGTCACCCTGCTTGCGGTTGAACGGCGTGGAACCGCGGCCAAACGCCGGGTCCTGGCCGTTACGGGCATGGGTATTGAAGTCAGTGACCGTGGCGCGCAGCCCGGACGGGTCAATGCCGCAGGCCTGGGCGAGCTGCTCGAGCGTAGCGCCGCGTTTCAGGTAGCCGTTGCGCACATGCGGCCCTACCGGCAGCGGTGCCGGGCGGGCATGCCCCAGGCCATAGCGGCGCAAGAAGCGATGATCGCAAACCAGCCACGAGCACACCTCCTGGCCCTGCGGCACTG  >  NC_002947/2972447‑2972735

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