Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,613,319	T→C	10.0%	L59L (TTA→CTA)	PP_4084 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,613,319	0	T	C	10.0%	106.0 / 3.4	40	L59L (TTA→CTA)	PP_4084	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/19);  new base C (2/2);  total (19/21)
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.15e-01

AATCCTCAATCTGGGTGATGAAGCCCACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTC  >  NC_002947/4613171‑4613451                                                                                                                                                     |                                                                                                                                     aaTCCTCAATCTGGGTGATGTAGCCCACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTAt                                                                                                                                      >  3:2873/1‑149 (MQ=255)                 gatgaAGCCCACGAGGACTCGACGCTGAACATGCTGCGCTTTCCAGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCa                                                                                                                        >  8:232094/1‑147 (MQ=255)                  atgaAGCCCACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAg                                                                                                                     <  1:43704/149‑1 (MQ=255)                      aGCCCACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGgcgc                                                                                                                 >  3:302397/1‑149 (MQ=255)                         ccACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGa                                                                                                              >  4:95933/1‑149 (MQ=255)                         ccACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGa                                                                                                              >  5:235520/1‑149 (MQ=255)                          cACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGAc                                                                                                             >  6:14614/1‑149 (MQ=255)                          cACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGAc                                                                                                             <  5:14614/149‑1 (MQ=255)                                   tCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGcc                                                                                                    >  4:314682/1‑149 (MQ=255)                                        gCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGaagaag                                                                                                                           <  8:309585/121‑1 (MQ=255)                                            aaCATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGaaa                                                                                             <  1:237020/147‑1 (MQ=255)                                            aaCATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTCGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGt                                                                                           >  6:61903/1‑149 (MQ=255)                                             aCATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGt                                                                                           >  4:337224/1‑148 (MQ=255)                                               aTGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGAc                                                                                        <  3:314682/149‑1 (MQ=255)                                                   tgcGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTACTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGc                                                                                    >  3:218211/1‑149 (MQ=255)                                                   tgcGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTACTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGc                                                                                    >  3:222583/1‑149 (MQ=255)                                                     cgcTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGcc                                                                                   <  4:2873/148‑1 (MQ=255)                                                       cTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTg                                                                                <  3:95933/149‑1 (MQ=255)                                                       cTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTg                                                                                <  5:274372/149‑1 (MQ=255)                                                                 aGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGa                                                                              <  3:180740/141‑1 (MQ=255)                                                                 aGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGa                                                                              >  4:180740/1‑141 (MQ=255)                                                                   gCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAg                                                                     >  6:119508/1‑148 (MQ=255)                                                                   gCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGc                                                                    <  6:235520/149‑1 (MQ=255)                                                                            gTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCg                                                           <  3:337224/149‑1 (MQ=255)                                                                                   aTCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTACTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTc                                                     <  4:218211/148‑1 (MQ=255)                                                                                   aTCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTACTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTc                                                     <  4:222583/148‑1 (MQ=255)                                                                                      gctcgcAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCAtt                                                                                                         >  1:335340/1‑93 (MQ=255)                                                                                      gctcgcAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCAtt                                                                                                         <  2:335340/93‑1 (MQ=255)                                                                                        tcgcAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGaaa                                               >  1:284870/1‑149 (MQ=255)                                                                                          gcAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATg                                             <  4:302397/149‑1 (MQ=255)                                                                                                 tGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATc                                       >  3:162750/1‑148 (MQ=255)                                                                                                      cAACCACCATGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTTTTACGGAAGAAGTCAATGGCGCTGACCATTACGCCCGCGAAAGTTACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGCGAAATg                                             <  5:48559/137‑1 (MQ=255)                                                                                                      cAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGCAGCTGACCATTACGCCCGCGAAAGTCACACGCCTTGGAAGCAATGACCATCCCCTCGACAGTCGGCAAATg                                             >  6:48559/1‑137 (MQ=255)                                                                                                       aaCCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTc                                                            <  2:104654/121‑1 (MQ=255)                                                                                                       aaCCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTc                                                            >  1:104654/1‑121 (MQ=255)                                                                                                         ccaccaAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATg                              <  6:73059/149‑1 (MQ=255)                                                                                                                       gATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGcc                <  7:232094/149‑1 (MQ=255)                                                                                                                               cTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTAc         >  8:89820/1‑148 (MQ=255)                                                                                                                                     ggagCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAg      >  5:303226/1‑145 (MQ=255)                                                                                                                                     ggagCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAg      <  6:303226/145‑1 (MQ=255)                                                                                                                                     ggagCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTc  <  4:162750/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                     AATCCTCAATCTGGGTGATGAAGCCCACGAGGACTCGACGCTGAACATGCTGCGCTTTCCCGACAGGCAAACCACGTTGGTGATCGCTCGCAGCCCTGTTGCAACCACCAAGACCCTGGATGAGGCCTTGGAGGAGCAACTGAACGTATTACGGAAGAAGTCAAAGGCGCTGACCATTACGCCCGCGAAAGTGACCCGCCTTGGAAGCAATGAGCATCCCGTCGACAGTCGGGAAATGGCCATCCAGTTCATGGTCGGCGACAAGCCCCACTACCAGCTTC  >  NC_002947/4613171‑4613451

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