Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,113,808	A→G	9.3%	Y151C (TAC→TGC)	PP_3619 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,113,808	0	A	G	9.3%	112.4 / 3.1	43	Y151C (TAC→TGC)	PP_3619	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

ATAAAAAAACTCTCTGACATATACCTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCA  >  NC_002947/4113672‑4113952                                                                                                                                         |                                                                                                                                                 ataaAAAAACTCTCTGACATATACCTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAg                                                                                                                                      <  4:196385/149‑1 (MQ=255)       aaaCTCTCTGACATATACCTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGc                                                                                                                                <  2:268268/149‑1 (MQ=255)       aaaCTCTCTGACATATACCTCAATAGAGCCACCAAAAAAATATATAACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGc                                                                                                                                >  1:268268/1‑149 (MQ=255)             tctGACATATACCTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCg                                                                                                                           >  7:221464/1‑148 (MQ=255)                     ataCCTCAATAGAGCCACCTAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCATCGCCTTCGccc                                                                                                                  >  6:127101/1‑149 (MQ=255)                        ccTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAAATCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGACCtg                                                                                                                            >  8:22037/1‑135 (MQ=255)                         cTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTAt                                                                                                                                  >  6:343545/1‑129 (MQ=255)                         cTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTAt                                                                                                                                  <  5:343545/129‑1 (MQ=255)                         cTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCccg                                                                                                              >  2:4131/1‑149 (MQ=255)                         cTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTAGGCGGCTCCTACTTCACCACCAGCTATGCCGACGCATTAGCCCCccg                                                                                                              >  4:289825/1‑149 (MQ=255)                                       caAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTGCTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACgg                                                                                                  >  5:152577/1‑147 (MQ=255)                                       caAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTGCTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACgg                                                                                                  >  5:146829/1‑147 (MQ=255)                                             aaTATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGc                                                                                          <  8:147037/149‑1 (MQ=255)                                               tatatTACAAACTCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCAc                                                                                        <  3:289825/149‑1 (MQ=255)                                                 tatTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCAcc                                                                                       <  1:4131/148‑1 (MQ=255)                                                   ttACAAACGCATAAATACACTCATAGTTATTCATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCaa                                                                                     >  1:39040/1‑148 (MQ=255)                                                       aaaCGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGa                                                                                >  1:198131/1‑149 (MQ=255)                                                            cATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGcc                                                                           <  1:27434/149‑1 (MQ=255)                                                              tAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTGCTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCtg                                                                         <  6:152577/149‑1 (MQ=255)                                                              tAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTGCTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCtg                                                                         <  6:146829/149‑1 (MQ=255)                                                               aaaTACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCtgt                                                                        >  6:30312/1‑149 (MQ=255)                                                                 aTACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACGACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTgtgg                                                                      <  5:30312/149‑1 (MQ=255)                                                                 aTACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCAc                                                                                        <  2:80514/131‑1 (MQ=255)                                                                     acTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTgtgggtg                                                                   >  4:176915/1‑148 (MQ=255)                                                                         aTAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAg                                                              >  4:2227/1‑149 (MQ=255)                                                                                ttgattgGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTACTACTTCACAACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCATGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGaa                                                       >  3:216944/1‑149 (MQ=255)                                                                                   attgGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCaccacc                                                                                                                                        >  8:86599/1‑65 (MQ=255)                                                                                   attgGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCaccacc                                                                                                                                        <  7:86599/65‑1 (MQ=255)                                                                                          tAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAAc                                             >  7:121463/1‑149 (MQ=255)                                                                                                      cGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGc                                 <  7:128480/149‑1 (MQ=255)                                                                                                        gTATTTTCAAACGCAGAGATTTCGTCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGGCCCCCGCCGGTAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCAt                               <  8:221464/149‑1 (MQ=17)                                                                                                        gTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCATCAAAAAGACTGCCTCTGGGTCGACAGTAACGAACCCACCGAACCCGGCGTCAAGCAt                               >  6:321997/1‑149 (MQ=255)                                                                                                        gTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCAt                               <  3:2227/149‑1 (MQ=255)                                                                                                          aTTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGt                             >  3:102431/1‑149 (MQ=255)                                                                                                                aaaCGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCAcc                                                 >  6:228973/1‑123 (MQ=255)                                                                                                                aaaCGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCAcc                                                 <  5:228973/123‑1 (MQ=255)                                                                                                                aaaCGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCg                        <  8:121463/148‑1 (MQ=255)                                                                                                                 aaCGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCTCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCg                        >  8:25535/1‑147 (MQ=255)                                                                                                                       caGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCTCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTgg                <  4:216944/149‑1 (MQ=255)                                                                                                                       caGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAATAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTgg                <  5:127101/149‑1 (MQ=255)                                                                                                                                    cTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTgcgc    >  4:355672/1‑148 (MQ=255)                                                                                                                                     tCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCa  >  5:178916/1‑149 (MQ=255)                                                                                                                                     tCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCa  >  4:140645/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                 ATAAAAAAACTCTCTGACATATACCTCAATAGAGCCACCAAAAAAATATATTACAAACGCATAAATACACTCATAGTTATTGATTGGAATAACACTGGAGGCGGTATTTTCAAACGCACAGAATTCGGCGGCTCCTACTTCACCACCAGCTATGCCCTCGCCTTCGCCCCCCGCCGGGAAGACGGCAGCCTGCACCAAAAAGACTGCCTGTGGGTGGACAGCAACGAACCCACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCA  >  NC_002947/4113672‑4113952

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