Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,426,589	T→C	10.3%	C375C (TGT→TGC)	PP_4765 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,426,589	0	T	C	10.3%	93.5 / 3.7	41	C375C (TGT→TGC)	PP_4765	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (18/17);  new base C (2/2);  total (20/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 = 4.30e-01

TGCTTTGGCCCGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGC  >  NC_002947/5426445‑5426731                                                                                                                                                 |                                                                                                                                                 tGCTTTGGCCCGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCg                                                                                                                                              >  1:39120/1‑149 (MQ=255) tGCTTTGGCCCGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCg                                                                                                                                              <  6:178615/149‑1 (MQ=255)    tttGGCCCGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGc                                                                                                                                            >  4:1393/1‑148 (MQ=255)         cccGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCtgt                                                                                                                                                  <  3:79969/137‑1 (MQ=255)         cccGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCtgt                                                                                                                                                  >  4:79969/1‑137 (MQ=255)           cGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCt                                                                                                                                    >  4:275596/1‑149 (MQ=255)           cGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCt                                                                                                                                    <  4:239195/149‑1 (MQ=255)            gggCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTg                                                                                                                                   >  4:193230/1‑149 (MQ=255)                    aaGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTgg                                                                                                                             >  8:273885/1‑147 (MQ=255)                        cATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGgcg                                                                                                                       <  7:26705/149‑1 (MQ=255)                        cATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGgcg                                                                                                                       >  3:97583/1‑149 (MQ=255)                            gCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGa                                                                                                                   >  8:87556/1‑149 (MQ=255)                             cccATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGTATGAGGGGTTGATGGTCAAGGGCATGCCTGGGGTGTTCTGTGCCGGCGAGATGCTGGACTCGGAGGCGCCGAc                                                                                                                  >  6:229172/1‑149 (MQ=255)                               cATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCg                                                                                                                <  4:97583/149‑1 (MQ=255)                               cATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTCTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCg                                                                                                                >  2:131872/1‑149 (MQ=255)                                 tCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCTGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGCAGGCGCCGACCggg                                                                                                              <  5:58943/149‑1 (MQ=255)                                 tCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCggg                                                                                                              <  2:39120/149‑1 (MQ=255)                                     ccTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGt                                                                                                          >  2:35950/1‑149 (MQ=255)                                         gTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCt                                                                                                      <  3:275596/149‑1 (MQ=255)                                                  cgGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTGCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGGGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGc                                                                                             <  3:241654/149‑1 (MQ=255)                                                             tGAAGCGATGAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCa                                                                                  <  1:185996/149‑1 (MQ=255)                                                               aaGCGATCCGCTGTGCCGGCGGGGTGCGCTTCGAGGGGCTGTATCTGGGGTTGATGGTCAAGTGCATGCCGGGGGTGTTCTGGGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑g                                                                                <  5:199710/148‑1 (MQ=255)                                                                 gCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑gg                                                                            >  3:161199/1‑149 (MQ=255)                                                                          agTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGCGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCt                                                                                        >  1:13533/1‑130 (MQ=255)                                                                          agTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGCGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCt                                                                                        <  2:18953/130‑1 (MQ=255)                                                                          agTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGCGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCt                                                                                        <  2:13533/130‑1 (MQ=255)                                                                          agTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGCGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCt                                                                                        >  1:18953/1‑130 (MQ=255)                                                                                  cggGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACggg                                                           <  3:1393/149‑1 (MQ=255)                                                                                        ccGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAgcgg                                                     >  6:334135/1‑149 (MQ=255)                                                                                        ccGCTTCCTGGGCCTGGATGCGGGGTTGATGGGCACGGGCATGCCGGGGGGGGGCGGGGCCGGCGAGAGGCGGGACTGGGCGGCGCCGACCGGGGGGTAGCTGCGCACGGCGTGCTTTGCCA‑GC‑GGCTGGCGCGCCGGACGGGCAgcgg                                                     <  7:273885/149‑1 (MQ=255)                                                                                         cGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAgcggc                                                    <  1:35950/149‑1 (MQ=255)                                                                                           cTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑g                                                                                >  2:327296/1‑120 (MQ=255)                                                                                           cTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑g                                                                                <  1:327296/120‑1 (MQ=255)                                                                                           cTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAg                                                        >  1:176036/1‑143 (MQ=255)                                                                                                      tGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAg                                       >  1:14131/1‑149 (MQ=255)                                                                                                        gATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAg                                       <  7:87556/147‑1 (MQ=255)                                                                                                             ggggTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCt                                <  3:193230/149‑1 (MQ=255)                                                                                                                gTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCGGCGGCGGATTGGGTGCAGGCTTCCTGa                              >  2:199348/1‑148 (MQ=255)                                                                                                                        tCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTc                                  >  6:315540/1‑136 (MQ=255)                                                                                                                        tCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTc                                  <  5:315540/136‑1 (MQ=255)                                                                                                                          aaGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGTTTGGGTGCAGGCTTCCTGCCAGATTGCCgg                   >  8:289407/1‑149 (MQ=255)                                                                                                                                      ggggTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCGGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGa       <  1:199348/149‑1 (MQ=255)                                                                                                                                              cTGTGGCGGGGAGATGGTGGTCTGGGAGGCGCCGTCCGGGGGGTATCTGCTCACGGCGAGCTTTGCCATGCTTTTTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGc  <  1:279583/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                 TGCTTTGGCCCGGGCCATCAAGGCATTGCCCATCACCCTGGTACGCACGCGGCCGCTGGATGAAGCGATCAGCAGTGCCGGCGGGGTCCGCTTCGAGGGCCTGGATGAGGGGTTGATGGTCAAGGGCATGCCGGGGGTGTTCTGTGCCGGCGAGATGCTGGACTGGGAGGCGCCGACCGGGGGGTATCTGCTCACGGCGTGCTTTGCCA‑GC‑GGCTTGCGCGCCGGACGGGCAGCGGCGGATTGGGTGCAGGCTTCCTGACAGATTGCCGGGCTCATGCGGGAGCGGC  >  NC_002947/5426445‑5426731

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑