Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,328,949	G→T	13.0%	G66V (GGC→GTC)	PP_3799 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,328,949	0	G	T	13.0%	91.3 / 7.6	39	G66V (GGC→GTC)	PP_3799	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/9);  new base T (2/3);  total (26/13)
Fisher's exact test for biased strand distribution p-value = 3.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-01

CCAGCCCGGCGCTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTC  >  NC_002947/4328802‑4329091                                                                                                                                                    |                                                                                                                                               ccAGCCCGGCGCTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTcgg                                                                                                                                                <  6:81720/148‑1 (MQ=255)  cAGCCCGGCGCTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTcggcg                                                                                                                                              <  1:175038/149‑1 (MQ=255)  cAGCCCGGCGCTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTcggcg                                                                                                                                              >  7:145231/1‑149 (MQ=255)        ggCGCTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGc                                                                                                                                        <  8:84636/149‑1 (MQ=255)          cgcTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCtt                                                                                                                                      <  5:286288/149‑1 (MQ=255)          cgcTGGCGCACAACCCGATGTGCGAAGGCGAGGAACTGGCCATCGGCCAGGTCATGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCtt                                                                                                                                      <  2:375172/149‑1 (MQ=255)            cTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTgg                                                                                                                                    <  1:240236/149‑1 (MQ=255)             tGGCGCACAACCCGATGTGCGAATGCGAGGAACTGCCCAACGGCCAGGTCAGGTCAACGGGCGGTTACTCGGACGCCACCGGTCCCCCCGGGTGGACCGTGGACGTGATGGGTTATGACGACCAGTTGATGGACGGCGGCAGCCTTggg                                                                                                                                   >  6:34376/1‑148 (MQ=255)                 gcACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGcgac                                                                                                                                >  3:127577/1‑148 (MQ=255)                  cacaACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGCTTGCcgccga                                                                                                                              <  6:298226/149‑1 (MQ=255)                   acaACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCgacgac                                                                                                                             >  1:20563/1‑149 (MQ=255)                                    tGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTc                                                                                                              >  5:60984/1‑147 (MQ=255)                                       gAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAgcg                                                                                                         >  7:314320/1‑149 (MQ=255)                                           aaCTGCCCAACGCCCATGTCAAGGGCCCGGGGGGTTTGTCGGCGGGCTGGGGCGCCCCGGGGGTGACCTTGGTCGTGATCGTTTATTCCGAACTGGTGGTGGGGGGCGGAAAGCTTGGCTACGACTCGACGCTGACGTTCAAGCGCCCg                                                                                                     <  4:274898/149‑1 (MQ=255)                                              tGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGc                                                                                                                                        >  2:43025/1‑111 (MQ=255)                                              tGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGc                                                                                                                                        >  2:37173/1‑111 (MQ=255)                                              tGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGc                                                                                                                                        <  1:43025/111‑1 (MQ=255)                                              tGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGTCGCCAAGc                                                                                                                                        <  1:37173/111‑1 (MQ=255)                                              tGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGAc                                                                                                  >  7:223473/1‑149 (MQ=255)                                                gccgttggCCCGGTCATGTGCACGGGCGGTGCCCCGGCCGGCCGCGGGGCCCCGGGGGTGCCCGTGGACGTTATCGTTTATCACCACCCGGGGCGGGTCGGCGTCAAGCTGGGCCACGACTCGACGCTGACGTTCAAGCGCCCGGACgg                                                                                                <  8:285991/143‑1 (MQ=255)                                                 ccACCGGCCAGGTCAAGGGCCCGGGCGGTTTCGCGGTCGCCCGGGGTGCCCCCGGGGTTACCCTGCACGTGATCGTTTATGACTACCAGGTGCTGGTCGGCGGCATGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGc                                                                                               <  8:87731/149‑1 (MQ=255)                                                         cAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTa                                                                                       >  5:326079/1‑149 (MQ=255)                                                                       cggcggTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATg                                                                         >  7:125485/1‑149 (MQ=255)                                                                                         cAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTgg                                                       >  6:311928/1‑149 (MQ=255)                                                                                         cAGCGGTGCCCCCGGGGTGACCCTGGAAGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGCCAAGCTTGGCGACGACTCGACGCTGACCTTCAAGCGCCCGGACGGCGAATACTATGTGCTGTTCGATGCCGGCCCCGGGCACGTgg                                                       >  1:364145/1‑149 (MQ=255)                                                                                            cGGTGCCCCGGGGGAGACCCTGGCCGGGATCGGTTATCCCGAACCGGGGCGGGCCGGCGGCAAGCTGGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTg                                                    <  4:75374/149‑1 (MQ=255)                                                                                                 ccccgGGGGGGACCGTGGAGGTGATGGTTTATTACGAACCGGTGGGGGGCGGCGGCAAGCGTGCCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTg                                               <  5:34376/149‑1 (MQ=255)                                                                                                  ccccGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTg                                               >  7:28528/1‑148 (MQ=255)                                                                                                          tGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGcc                                      >  5:384882/1‑149 (MQ=255)                                                                                                             cccTGGACGTGATCGGTTATGTCGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGCCGTTCCAGCGCCCGGACCGCGAATTGTATGGGCTGTTAGATGCCGGCCCCGGGCACGTGGTTGCAGTGGACTACGCCGa                                    >  2:84643/1‑148 (MQ=255)                                                                                                               cTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGatat                                 >  3:107028/1‑149 (MQ=255)                                                                                                                tGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCAc                                                           >  7:340561/1‑122 (MQ=255)                                                                                                                  gACGTGTTCGGTAATGACGAACAGGTGCTGGTCGCCGGAAAGCTTGGCGACCACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGc                              <  2:117791/149‑1 (MQ=255)                                                                                                                   aCGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGcc                             >  4:226568/1‑149 (MQ=255)                                                                                                                         tCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGTTTGAAGTGGACTACGCCGATATCGCCCAACCa                       >  7:269161/1‑149 (MQ=255)                                                                                                                             ttATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAg                   >  6:393705/1‑149 (MQ=255)                                                                                                                             ttATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAg                   >  7:292252/1‑149 (MQ=255)                                                                                                                             ttATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAg                   >  6:12122/1‑149 (MQ=255)                                                                                                                                    gAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGcc            >  6:323414/1‑149 (MQ=255)                                                                                                                                    gAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGcc            >  1:248841/1‑149 (MQ=255)                                                                                                                                           tGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCa     >  2:125716/1‑149 (MQ=255)                                                                                                                                              tGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTc  >  3:334338/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               CCAGCCCGGCGCTGGCGCACAACCCGATGTGCGAATGCGAGGAACTGGCCAACGGCCAGGTCAAGTGCACCGGCGGTTTCTCGGACGGCAGCGGTGCCCCCGGGGTGACCCTGGACGTGATCGGTTATGACGAACAGGTGCTGGTCGGCGGCAAGCTTGGCGACGACTCGACGCTGACGTTCAAGCGCCCGGACGGCGAATTCTATGTGCTGTTCGATGCCGGCCCCGGGCACGTGGTTGAAGTGGACTACGCCGATATCGCCCAACCATGAGCCGCGCCCAAGTCATTC  >  NC_002947/4328802‑4329091

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