Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,531,110	A→C	10.1%	I117L (ATA→CTA)	PP_2219 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,531,110	0	A	C	10.1%	133.5 / 2.4	46	I117L (ATA→CTA)	PP_2219	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base C (2/3);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGATGAAA  >  NC_002947/2530967‑2531258                                                                                                                                                |                                                                                                                                                       agctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAg                                                                                                                                                     <  8:323675/147‑1 (MQ=255)   ctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCCTa                                                                                                                                                      <  5:267200/144‑1 (MQ=255)   ctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCCTa                                                                                                                                                      >  6:267200/1‑144 (MQ=255)   ctgGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCAt                                                                                                                                                 <  2:310022/149‑1 (MQ=255)         tGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCCTATGCATCgccag                                                                                                                                           <  1:52458/149‑5 (MQ=255)            cGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGca                                                                                                                                         >  5:36674/1‑148 (MQ=255)                          tgAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGTCGATGCCTATGCATCgcca                                                                                                                                            <  5:95865/131‑4 (MQ=255)                          tgAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGTCGATGCCTATGCATCgcca                                                                                                                                            >  6:95865/1‑128 (MQ=255)                                      gCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑ttttt                                                                                                             >  1:99640/1‑148 (MQ=255)                                       cAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCTCATGAATCGCCCA‑AATGCACA‑TTTTTTc                                                                                                           <  6:90252/149‑1 (MQ=255)                                        aaaTTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCa                                                                                                          <  2:99640/149‑1 (MQ=255)                                                  ccATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATc                                                                                                >  1:338832/1‑149 (MQ=255)                                                          aGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGAt                                                                                             >  4:8397/1‑144 (MQ=255)                                                          aGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGAt                                                                                             <  3:8397/144‑1 (MQ=255)                                                            ttGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAg                                                                                      >  2:268738/1‑149 (MQ=255)                                                              gCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAgg                                                                                     >  5:190348/1‑148 (MQ=255)                                                                 ccACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACa                                                                                 <  2:338832/149‑1 (MQ=255)                                                                  cACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAg                                                                                >  1:226829/1‑149 (MQ=255)                                                                   aCGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGAACAAGGAACAGt                                                                               >  8:144389/1‑149 (MQ=255)                                                                     gcggcgTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGa                                                                                              <  4:190115/132‑1 (MQ=255)                                                                     gcggcgTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGa                                                                                              >  3:190115/1‑132 (MQ=255)                                                                         cgTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGAt                                                                                             >  6:248090/1‑129 (MQ=255)                                                                         cgTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGAt                                                                                             <  5:248090/129‑1 (MQ=255)                                                                             gcgcCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGTCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCACAATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGcc                                                                       <  1:268738/148‑1 (MQ=255)                                                                             gcgcCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCAt                                                                     >  1:344614/1‑149 (MQ=255)                                                                             gcgcCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCAt                                                                     <  6:190348/149‑1 (MQ=255)                                                                             gcgcCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAAATCAGCTGGCGATGAATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATATATAAATCGATGCCCAAGGAACAGTTGCTGGCCAt                                                                     >  1:90024/1‑149 (MQ=255)                                                                               gcCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGt                                                                   >  8:8909/1‑149 (MQ=255)                                                                                  ttAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACATTTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTc                                                                  >  5:73404/1‑148 (MQ=255)                                                                                    aGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGc                                                               <  7:144389/148‑1 (MQ=255)                                                                                    aGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGc                                                               <  6:36674/148‑1 (MQ=255)                                                                                     gACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACATTTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCgac                                                               <  6:73404/148‑3 (MQ=255)                                                                                        tGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGAt                                                           >  4:75527/1‑148 (MQ=255)                                                                                             cccccTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTAc                                                     <  8:212616/149‑1 (MQ=255)                                                                                                     aaaaGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGc                                             >  3:158247/1‑149 (MQ=255)                                                                                                       aaGAGTTCCTGTTGGCGCTTGAAGAATTCAGCTGGCGATGCATAGTCATCGACGGCACATGAATCGCCAA‑AATGCACA‑TTTTTTCATAGAGAAATAGATGAACAATGAACAGTTGATGGCCATGTAGGCCGATCACTACATGAACGcc                                            >  8:89136/1‑148 (MQ=255)                                                                                                                     gcgcTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCctt                              <  4:158247/148‑1 (MQ=255)                                                                                                                           gaagaaTTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCg                                                                 >  7:46082/1‑107 (MQ=255)                                                                                                                           gaagaaTTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCg                                                                 <  8:46082/107‑1 (MQ=255)                                                                                                                           gaagaaTTCCGCTGGCGATGCATACGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTa                         <  2:90024/147‑1 (MQ=255)                                                                                                                                 ttCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAgctg                 >  5:76886/1‑149 (MQ=255)                                                                                                                                 ttCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAgctg                 >  7:177888/1‑149 (MQ=255)                                                                                                                                    cGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCtt                                   >  4:58163/1‑128 (MQ=255)                                                                                                                                    cGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCtt                                   <  3:58163/128‑1 (MQ=255)                                                                                                                                      cTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCa            <  1:302797/149‑1 (MQ=255)                                                                                                                                              gCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGATg     >  1:216159/1‑148 (MQ=255)                                                                                                                                                aTAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGATGaaa  <  2:216159/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                       AGCTGGGGTGCCGACCGGCGGATTGTGAGCAGTGGAAGCAAATTCGCACCCATTTGCAGTTGCTCCACGCGGCGTGGCGCCTTAGACTGCCGCCCCCTGTAAAAGAGTGCCGGTTGGCGCTTGAAGAATTCCGCTGGCGATGCATAGGCATCGACGGCACATGAATCGCCCA‑AATGCACA‑TTTTTTCATAGAGAAATCGATGACCAAGGAACAGTTGCTGGCCATGTCGGCCGATGACTACATGAACGCCGACCAGCTTGCCTTCTTTACTGAGCTGCTGCAGGCGATGAAA  >  NC_002947/2530967‑2531258

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