Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,621,557	Δ1 bp	8.2%	coding (729/957 nt)	PP_4940 ←	glycosyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,621,557	0	G	.	8.2%	151.2 / 7.7	49	coding (729/957 nt)	PP_4940	glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (22/23);  new base . (2/2);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

CTACGTCCAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGT  >  NC_002947/5621414‑5621702                                                                                                                                                |                                                                                                                                                   cTACGTCCAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCAt                                                                                                                                               >  3:265535/1‑149 (MQ=255)    cGTCCAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTc                                                                                                                                            <  7:26987/149‑1 (MQ=255)    cGTCCAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTc                                                                                                                                            <  4:252534/149‑1 (MQ=255)        cAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTg                                                                                                                                         <  2:261051/148‑1 (MQ=255)        cAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGc                                                                                                                                        <  6:293522/149‑1 (MQ=255)           cTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGa                                                                                                                                     >  3:167139/1‑149 (MQ=255)                     aCTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCAc                                                                                                                           >  6:262975/1‑149 (MQ=255)                     aCTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTT‑TGCATGTCTTGCAGACCCAgc                                                                                                                               <  4:287223/144‑1 (MQ=255)                     aCTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTT‑TGCATGTCTTGCAGACCCAgc                                                                                                                               >  3:282699/1‑144 (MQ=255)                     aCTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTT‑TGCATGTCTTGCAGACCCAgc                                                                                                                               >  3:287223/1‑144 (MQ=255)                     aCTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTT‑TGCATGTCTTGCAGACCCAgc                                                                                                                               <  4:282699/144‑1 (MQ=255)                        aaaaCGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGAGAGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAg                                                                                                                                >  5:213714/1‑141 (MQ=255)                         aaaCGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACgg                                                                                                                         <  5:162868/147‑1 (MQ=255)                                      ccAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTTGCCt                                                                                                         <  2:294026/149‑1 (MQ=255)                                       cAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCtt                                                                                                        >  2:70718/1‑149 (MQ=255)                                          aTGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGt                                                                                                     <  1:70718/149‑1 (MQ=255)                                               ggcCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTa                                                                                                 <  4:259369/148‑1 (MQ=255)                                               ggcCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTa                                                                                                 >  3:259369/1‑148 (MQ=255)                                                   aTTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCa                                                                                            <  2:97788/149‑1 (MQ=255)                                                      ttCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGc                                                                                          <  6:287072/148‑1 (MQ=255)                                                        cTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCaa                                                                                        >  2:84018/1‑148 (MQ=255)                                                          gCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAaca                                                                                      >  2:118458/1‑148 (MQ=255)                                                             aTCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAAcacaca                                                                                  <  5:262975/149‑1 (MQ=255)                                                             aTCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAAcacaca                                                                                  <  5:103683/149‑1 (MQ=255)                                                                      gCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGc                                                                         >  5:202384/1‑149 (MQ=255)                                                                          gCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATg                                                                     >  3:64499/1‑149 (MQ=255)                                                                              ggggTTCGCGCCCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCtt                                                                 >  4:274025/1‑149 (MQ=255)                                                                                 gTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CGTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCg                                                               >  8:287152/1‑148 (MQ=255)                                                                                  ttCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGaa                                                             <  6:202384/149‑1 (MQ=255)                                                                                     gcgcgcAAACGACGGAGCTTTTACTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCATACCCAGCGCACGGTTCTCCTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAg                                                            <  1:84018/148‑1 (MQ=255)                                                                                       gcgcAAACGACGGAGCTTTTCCTGAACGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCt                                                        <  4:167139/149‑1 (MQ=255)                                                                                          cAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGa                                                     <  4:265535/149‑1 (MQ=255)                                                                                                             tGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGcc                                  >  6:27464/1‑149 (MQ=255)                                                                                                             tGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGcc                                  >  5:295048/1‑149 (MQ=255)                                                                                                             tGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGcc                                  >  6:222747/1‑149 (MQ=255)                                                                                                              gAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGcc                                                                        <  4:44201/110‑1 (MQ=255)                                                                                                              gAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGcc                                                                        >  3:44201/1‑110 (MQ=255)                                                                                                                               tCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTGGCCCATGCCGACATCGGCACTGACg                  >  6:226301/1‑147 (MQ=255)                                                                                                                                    gAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAAc           <  4:64499/149‑1 (MQ=255)                                                                                                                                    gAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAAc           <  6:295048/149‑1 (MQ=255)                                                                                                                                    gAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAAc           <  1:118458/149‑1 (MQ=255)                                                                                                                                     aaCACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAAc           >  4:108562/1‑148 (MQ=255)                                                                                                                                       caccacGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGaa        <  5:27464/149‑1 (MQ=255)                                                                                                                                       caccacGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGa         >  8:2564/1‑148 (MQ=255)                                                                                                                                          cacGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATa      >  5:28822/1‑148 (MQ=255)                                                                                                                                          cacGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATa      <  6:28822/148‑1 (MQ=255)                                                                                                                                          cacGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATAc     >  7:240560/1‑149 (MQ=255)                                                                                                                                           acGTTGTGCATGTCTCGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACg    <  5:226301/149‑1 (MQ=255)                                                                                                                                              ttGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGt  >  4:262774/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                   CTACGTCCAACTGCGCCGAAACTAAAACGGGAAACCGCCAGATGGCGGCCATTTTCTGCAATCTGCCGGGCTAGCTCGGGGTTCGCGCGCAAACGACGGAGCTTTTCCTGAAGGGTCGGGATGTTGTCATAGAACACCACGTTGTGCATGTCTTGCAGACCCAGCGCACGGTTCT‑CTTCCTCGCCTTGGTCGTATGCCAGCAACACACAGCCACAGGCCATGGCTTCGAAGTTCTTGATCATGTATTCGCCCATGCCGACATCGGCACTGACGAAGAAACGAATACGGT  >  NC_002947/5621414‑5621702

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