Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,002,014	G→C	10.3%	G652A (GGC→GCC)	PP_1786 →	glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,002,014	0	G	C	10.3%	105.7 / 4.3	39	G652A (GGC→GCC)	PP_1786	glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (20/15);  new base C (2/2);  total (22/17)
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.64e-01

ATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTACGGCAA  >  NC_002947/2001868‑2002159                                                                                                                                                   |                                                                                                                                                  aTACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGgct                                                                                                                                                 <  6:323781/149‑1 (MQ=255)         agcagTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGCCTGCTTcaaa                                                                                                                                         <  3:170472/149‑4 (MQ=255)         agcagTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGCCTGCTTcaa                                                                                                                                          >  4:170472/1‑146 (MQ=255)          gcagTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGCCTGCTTcaaac                                                                                                                                        <  4:88667/149‑5 (MQ=255)                        tGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGtttt                                                                                                                           <  6:29419/148‑1 (MQ=255)                         gCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGtttt                                                                                                                           >  6:22238/1‑147 (MQ=255)                                             aaGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGt                                                                                                     >  6:281522/1‑149 (MQ=255)                                                             gctTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTc                                                                                            >  6:259560/1‑142 (MQ=255)                                                             gctTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTc                                                                                            <  5:259560/142‑1 (MQ=255)                                                                   gCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGATTTGAGCAGCACATCCAATTACCGTGACGTTGTCAAAGATGgag                                                                                  >  3:255747/1‑146 (MQ=255)                                                                          gATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGtt                                                                         >  8:328018/1‑148 (MQ=255)                                                                             atCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCt                                                                     >  2:303341/1‑149 (MQ=255)                                                                              tCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGgaga                                                                                 >  3:114946/1‑136 (MQ=255)                                                                              tCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGgaga                                                                                 <  4:114946/136‑1 (MQ=255)                                                                                aaTATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGc                                                                  >  5:191571/1‑149 (MQ=255)                                                                                      gCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCacta                                                             <  5:22238/148‑1 (MQ=255)                                                                                            ttGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGagaa                                                                                >  4:162860/1‑123 (MQ=255)                                                                                            ttGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGagaa                                                                                <  3:162860/123‑1 (MQ=255)                                                                                              gCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGCTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGATTCCTGGCCACTACTgaag                                                       >  5:113872/1‑146 (MQ=255)                                                                                               cACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTgaagaag                                                    >  8:165121/1‑148 (MQ=255)                                                                                                 cGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAAt                                                 <  6:191571/149‑1 (MQ=255)                                                                                                 cGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAAt                                                 <  5:281522/149‑1 (MQ=255)                                                                                                  gATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATg                                                <  7:328018/149‑1 (MQ=255)                                                                                                        aaaGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCg                                          >  3:11270/1‑149 (MQ=255)                                                                                                              aaTGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCtt                                     <  4:11270/148‑1 (MQ=255)                                                                                                                      tAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTaaaaaaa                             <  7:165121/148‑1 (MQ=255)                                                                                                                           gCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCt                       >  6:20675/1‑149 (MQ=255)                                                                                                                             gAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGa                     >  8:244/1‑149 (MQ=255)                                                                                                                               gCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATa                   <  7:244/149‑1 (MQ=255)                                                                                                                                   aaaTGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACTGCCAATTACCGTGACGTTGTCAAAGATGGTGATGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGca                >  7:276274/1‑147 (MQ=255)                                                                                                                                   aaaTGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGc                        >  2:155328/1‑140 (MQ=255)                                                                                                                                   aaaTGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGc                        <  1:155328/140‑1 (MQ=255)                                                                                                                                   aaaTGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGc                 <  8:30632/147‑1 (MQ=255)                                                                                                                                     aTGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAg             >  7:228151/1‑149 (MQ=255)                                                                                                                                          ccgAAGCAGCCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTa        >  8:77010/3‑149 (MQ=255)                                                                                                                                           tGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTAc       >  3:115526/1‑149 (MQ=255)                                                                                                                                             aaGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTACgg     >  4:310974/1‑149 (MQ=255)                                                                                                                                              aGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTACgg     <  3:345776/148‑1 (MQ=255)                                                                                                                                                cAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATg                                                <  4:224749/103‑1 (MQ=255)                                                                                                                                                cAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATg                                                >  3:224749/1‑103 (MQ=255)                                                                                                                                                  ggCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTACGGCaa  >  3:254740/1‑147 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ATACGCTCAGCAGTTTACCCAGTTGCCGGCCATGGACAGTGTGCAAGGTTATTGGTCGCTGCTTGAGCAAGCGGATATCAATATTGCGGTGTTGCACGATGACAAAGTCAATGCCTGTAAAAGCGAGCTGAAATGGTTTGAAGCAGGCTGCTTCGGTATCCCATCGGTTTTGAGCAGCACAGCCAATTACCGTGACGTTGTCAAAGATGGAGAAGATGGTTTCCTGGCCACTACTGAAGAAGAATGGTATCGGGCTTTAAAAAAACTCGCTGATAGCCAAGCCTTACGGCAA  >  NC_002947/2001868‑2002159

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