Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,282,287	C→T	11.1%	E27K (GAG→AAG)	PP_2883 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,282,287	0	C	T	11.1%	92.7 / 3.5	36	E27K (GAG→AAG)	PP_2883	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/15);  new base T (2/2);  total (19/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.39e-01

CTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAA  >  NC_002947/3282141‑3282431                                                                                                                                                   |                                                                                                                                                  cTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCaa                                                                                                                                                 <  1:14759/149‑1 (MQ=255)   gCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAgg                                                                                                                                               >  7:159721/1‑149 (MQ=255)      aGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAggtgg                                                                                                                                            >  4:121760/1‑149 (MQ=255)              aGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCA‑CCACCTTGAGGTGGCTGCACGcc                                                                                                                                   <  8:226145/149‑1 (MQ=255)                                tCGGGGGCTTGGCAACGAACGCAATCCGCCAGGCCTGCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAg                                                                                                                 <  6:318850/149‑1 (MQ=255)                                             aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCCCGCCAAGGGCCACATTGCGCAGGCACTGGTAGACa                                                                                                    >  2:139854/1‑149 (MQ=255)                                             aacgaacgCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACa                                                                                                    >  6:241192/1‑149 (MQ=255)                                                      aaTGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCg                                                                                           >  8:163189/1‑149 (MQ=255)                                                           gcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGtt                                                                                         >  6:260580/1‑146 (MQ=255)                                                           gcCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGtt                                                                                         <  5:260580/146‑1 (MQ=255)                                                            ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgcgc                                                                                                                   <  7:34840/119‑1 (MQ=255)                                                            ccAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTgcgc                                                                                                                   >  8:34840/1‑119 (MQ=255)                                                                gCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAAc                                                                                 <  3:121760/149‑1 (MQ=255)                                                                  cTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAg                                                                                                                 >  8:5750/1‑115 (MQ=255)                                                                       gTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCtgtg                                                                          <  7:163189/149‑1 (MQ=255)                                                                        tATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTgg                                                                         <  5:275461/149‑1 (MQ=255)                                                                            cAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                      <  5:241192/148‑1 (MQ=255)                                                                              gggCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTc                                                                      >  2:62150/1‑146 (MQ=255)                                                                                 ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGAt                                                                 >  5:275747/1‑148 (MQ=255)                                                                                 ccAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGAt                                                                 >  3:295302/1‑148 (MQ=255)                                                                                               cAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGc                                                  >  5:208869/1‑149 (MQ=255)                                                                                                    cTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGTAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGtt                                             >  6:294322/1‑149 (MQ=255)                                                                                                       gATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTg                                            <  4:295302/147‑1 (MQ=255)                                                                                                         tGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGgcgc                                        <  6:275747/149‑1 (MQ=255)                                                                                                            cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                     >  5:220823/1‑149 (MQ=255)                                                                                                            cGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtga                                     >  1:51040/1‑149 (MQ=255)                                                                                                             gACATCGATCAGGCCGCCATCGCGGCGTGTAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCtgat                                    <  5:294322/149‑1 (MQ=255)                                                                                                                   gATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGAtctc                              >  3:39104/1‑149 (MQ=255)                                                                                                                       aGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGc                           <  1:62150/148‑1 (MQ=255)                                                                                                                            gccATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATc                     <  1:139854/149‑1 (MQ=255)                                                                                                                                   cggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTAGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg              <  8:74127/149‑1 (MQ=255)                                                                                                                                    ggcgTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGa             <  2:51040/149‑1 (MQ=255)                                                                                                                                    ggcgTGCA‑CCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGAGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCg              >  6:38533/1‑148 (MQ=255)                                                                                                                                     gcgTGCA‑CCACCTTGAGGTGGCTGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGa             <  5:38533/148‑1 (MQ=255)                                                                                                                                             ccACCTTGAGGTGGCTGCACGCCGAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCtt      >  4:309399/1‑148 (MQ=255)                                                                                                                                               aCCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGaaa  >  8:195421/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CTGCCAGTGTTCGAGGCCGGCCTGTTCGCCGTCGGGGGCTTGGCAACGAACGCAATGCGCCAGGCCTTCGGTATCCAGGGCCAACGTCAGTTGCCAGCCCTGGATGCCGACATCGATCAGGCCGCCATCGCGGCGTGCAGCCACCTCAAGGTGG‑TGCACGCCAAGGGCCACATTGCGCAGGCACTGGTAGACATCCCGGGCGTTCGGTGAACGTCTGTGGGTCAAGATCAGAAATTCGCCGGCGTGTTGGCGCTGATGATCTCGGCCTCATCCGGGCCGATGTTCTTGAAA  >  NC_002947/3282141‑3282431

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