Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,552,442:1	+G	15.2%	coding (56/978 nt)	PP_3137 →	group 2 family glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,552,442	1	.	G	15.2%	98.2 / 16.4	33	coding (56/978 nt)	PP_3137	group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (13/15);  new base G (3/2);  total (16/17)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01

TGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCG  >  NC_002947/3552293‑3552588                                                                                                                                                       |                                                                                                                                                    tGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGa                                                                                                                                                      <  1:822604/150‑1 (MQ=255)        tgtgCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGATTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑Gc                                                                                                                                             <  2:378320/150‑1 (MQ=255)            ccTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTg                                                                                                                                         <  2:169298/150‑1 (MQ=255)             cTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCG‑CTCGAGGCGCT‑g                                                                                                                                              <  1:173325/144‑1 (MQ=255)             cTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCG‑CTCGAGGCGCT‑g                                                                                                                                              >  2:173325/1‑144 (MQ=255)                        cTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGgcgc                                                                                                                             >  2:16799/1‑150 (MQ=255)                         tGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCa                                                                                                                            >  1:670670/1‑150 (MQ=255)                                gggggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGcc                                                                                                                     <  2:293380/150‑1 (MQ=255)                                  gggTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGa                                                                                                                   >  2:468491/1‑150 (MQ=255)                                                atcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGAT‑GCTTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACg                                                                                                     >  1:615559/1‑150 (MQ=255)                                                 tcaAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGt                                                                                                    >  1:706021/1‑150 (MQ=255)                                                     ggCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGAc                                                                                                <  2:615559/150‑1 (MQ=255)                                                             aTCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGACGACAAAGTCATGTACGTCGACTCGGGTTc                                                                                        <  1:501624/150‑1 (MQ=255)                                                                 ggCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCAcc                                                                                    <  2:706021/150‑1 (MQ=255)                                                                   ccgccgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGAGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGa                                                                                  >  1:492438/1‑150 (MQ=255)                                                                    cgccgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGTGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGAt                                                                                 <  2:739571/150‑1 (MQ=255)                                                                       cgTGGCATCGGTCGAGTGAGGCCTATGGCAACGGAGATCGGGGAGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTCTATCAACGCCGACCGGGCGCGACAAAGTCATTTACGTCGACCCGGGTTCCACCGATGGt                                                                              >  2:663077/1‑150 (MQ=255)                                                                                 gTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTAATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGt                                                                    >  1:611854/1‑150 (MQ=255)                                                                                      gTGAGGCATATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCCCT‑GCCTTGTCTCGCTGGCGCACGCGGCCGACAAAGTCATGTACGTCGACTGGCGGCCCACCGATGGTTCGCTGCAGTTGgcc                                                               <  1:663077/150‑1 (MQ=255)                                                                                                    aaCGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGt                                                                                          <  1:678008/109‑1 (MQ=255)                                                                                                    aaCGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGt                                                                                          >  2:678008/1‑109 (MQ=255)                                                                                                    aaCGGTGATCGGGGGGGTGGTAATCGGGCGCAACGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGt                                                 <  2:306910/150‑1 (MQ=255)                                                                                                               gggTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTgg                                      <  1:16799/150‑1 (MQ=255)                                                                                                                   ggtggtAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCt                                  <  1:468491/150‑1 (MQ=255)                                                                                                                      ggtAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGa                               >  1:745941/1‑150 (MQ=255)                                                                                                                        tAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACa                             <  2:745941/150‑1 (MQ=255)                                                                                                                           tCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGTCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGa                          <  2:670670/150‑1 (MQ=255)                                                                                                                                   aaTGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGt                                                                                                    <  2:533368/68‑1 (MQ=255)                                                                                                                                   aaTGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGt                                                                                                    >  1:533368/1‑68 (MQ=255)                                                                                                                                         aggcAGCG‑CTCGAGGCGCTGGCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAc             >  2:415468/4‑150 (MQ=255)                                                                                                                                          gCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCAccg           >  2:127567/1‑150 (MQ=255)                                                                                                                                           gcAGCG‑CTCGAGGCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTc                                                                                                   >  2:633762/2‑61 (MQ=38)                                                                                                                                           gcAGCG‑CTCGAGGCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTc                                                                                                   <  1:633762/60‑1 (MQ=255)                                                                                                                                                   tCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCAcgcg  >  2:586634/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                    TGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCGCCGCACGCG  >  NC_002947/3552293‑3552588

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