Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,552,437	Δ1 bp	14.7%	coding (51/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,437	0	C	.	14.7%	95.9 / 12.9	34	coding (51/978 nt)	PP_3137	group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (13/16);  new base . (3/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

TGCTGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCG  >  NC_002947/3552290‑3552579                                                                                                                                                    |                                                                                                                                                 tGCTGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCt                                                                                                                                                >  1:293380/1‑150 (MQ=255)    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)                                                                                                                                                    |                                                                                                                                                 TGCTGGTAGATGTGCCTGACAATTCACTGGCAGTGGGGGTGCCGGCCATCATCAAGGCCCGGCATCCGGCCGCCGTGGCATCGGTCGAGTGAGGCCTATGGCAACGGTGATCGGGGTGGTGGTAATCGGCCGCAATGAAGGCCAGCGCCTCGA‑GCGCT‑GCCTTGTTTCGCTGGCGCACGGGGCCGACAAAGTCATGTACGTCGACTCGGGTTCCACCGATGGTTCGCTGCAGTTGGCCGCCAACCTTGGCGTGGCAGTACTGGCGCTGGACATGAGTATCCCGTTCACCG  >  NC_002947/3552290‑3552579

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