Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,593,269	T→A	13.6%	I102N (ATC→AAC)	PP_3171 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,593,269	0	T	A	13.6%	84.3 / 4.1	35	I102N (ATC→AAC)	PP_3171	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/14);  new base A (2/3);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACG  >  NC_002947/3593123‑3593401                                                                                                                                                   |                                                                                                                                     ggCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTAt                                                                                                                                      <  2:91478/147‑1 (MQ=255)        aCGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctggg                                                                                                                             <  8:69554/149‑1 (MQ=255)        aCGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGATACGAcctggg                                                                                                                             <  5:82558/149‑1 (MQ=255)          gTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcc                                                                                                                                 <  8:152665/143‑1 (MQ=255)          gTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGAcc                                                                                                                                 >  7:152665/1‑143 (MQ=255)          gTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACACTCAATACATTGCACCCAGGTCGTATCGCcctgggcc                                                                                                                           >  1:187029/1‑149 (MQ=255)                    aTCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTg                                                                                                                  <  2:35962/148‑1 (MQ=255)                    aTCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTg                                                                                                                  <  4:75824/148‑1 (MQ=255)                       aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgcgc                                                                                                               >  8:187770/1‑148 (MQ=255)                       aaCGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgcgc                                                                                                               >  5:189468/1‑148 (MQ=255)                                   aaCGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCgg                                                                                                  <  7:187770/149‑1 (MQ=255)                                       ccACCCGACACATTCGTGTGGGTGCTGGAGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACAATGCACCCAGGTCGTAAAGACCTGGGCCTTGGCCGTGCGCCAGGCACGACCGGGcc                                                                                               >  4:245009/1‑148 (MQ=255)                                              acacATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgc                                                                                                                 >  4:86511/1‑123 (MQ=255)                                              acacATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTgc                                                                                                                 <  3:86511/123‑1 (MQ=255)                                                                 ggCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgccgc                                                                    >  1:160135/1‑149 (MQ=255)                                                                  gCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGAGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGgcgcagcg                                                                   >  5:127221/1‑149 (MQ=255)                                                                                 ccAACCATGCGCCGTTGGTGGTGGCCGAAGAGTTTGGCACCCTCAATACATTTCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                    <  6:239762/149‑1 (MQ=255)                                                                                 ccAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTc                                                    >  4:317322/1‑149 (MQ=255)                                                                                  cAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCg                                                   >  6:211479/1‑149 (MQ=255)                                                                                   aaCCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCgg                                                  >  2:190896/1‑149 (MQ=255)                                                                                   aaCCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCgg                                                  >  3:14921/1‑149 (MQ=255)                                                                                       aTGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAg                                              <  3:317322/149‑1 (MQ=255)                                                                                       aTGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAg                                              <  6:189468/149‑1 (MQ=255)                                                                                         gcgcCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGa                                            >  6:141112/1‑149 (MQ=255)                                                                                                       ggCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGc                              >  4:14572/1‑149 (MQ=255)                                                                                                         ccGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGt                            <  2:160135/149‑1 (MQ=255)                                                                                                               aGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTAc                      <  3:14572/149‑1 (MQ=255)                                                                                                               aGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTAc                      <  4:14921/149‑1 (MQ=255)                                                                                                                     gCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGcc                <  3:245009/149‑1 (MQ=255)                                                                                                                      cACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGGTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGct                >  7:135809/1‑147 (MQ=255)                                                                                                                           tCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAAc          >  5:184066/1‑149 (MQ=255)                                                                                                                             aaTACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACg         <  5:211479/148‑1 (MQ=255)                                                                                                                                  aTTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCTGACAACGGCAAAc   <  8:135809/149‑1 (MQ=255)                                                                                                                                   ccaggcccAGGTCGATACGACCTGGGTCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACa            <  1:31121/134‑1 (MQ=255)                                                                                                                                   ccaggcccAGGTCGATACGACCTGGGTCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACa            >  2:31121/6‑139 (MQ=255)                                                                                                                                   ccagaCCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACg  >  1:213612/5‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                     GGCTGCCACGTCGCTGGTGATCAACGAAATCGCCAACGCCACCCGACACATTCGTGTGGGTGCTGGCGGTATCATGCTGCCCAACCATGCGCCGTTGGTGGTGGCCGAACAGTTTGGCACCCTCAATACATTGCACCCAGGTCGTATCGACCTGGGCCTGGGCCGTGCGCCAGGCACGTCCGGGCCGACCGTTCGGGCATTGCGCGGCGCCGCGCCAGAGCGGGATTTCGGCCAGGACATCGCCGAGCTGCGTGATTACCTGGCCGACAACGGCAAACG  >  NC_002947/3593123‑3593401

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