Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	857,147	Δ1 bp	7.3%	coding (238/699 nt)	PP_0738 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	857,147	0	C	.	7.3%	175.8 / 6.8	55	coding (238/699 nt)	PP_0738	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base . (2/2);  total (28/27)
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.70e-01

TGCCCCTTTGTCAGCAACGGCAGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCACTCGCGGCAA  >  NC_002947/857002‑857293                                                                                                                                                  |                                                                                                                                                    tGCCCCTTTGTCAGCAACGGCAGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAg                                                                                                                                                  >  2:110395/1‑149 (MQ=255)    cccTTTGTCAGCAACGGCAGCGATTTTCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGATTCAGCCGGTCCAGCg                                                                                                                                                >  4:269311/1‑148 (MQ=255)           tCAGCAACGGCAGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAg                                                                                                                                       <  2:358289/149‑1 (MQ=255)              gcgaCGGCAGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCa                                                                                                                                    >  8:211627/4‑149 (MQ=255)                    gCAGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCa                                                                                                                               >  2:342482/1‑148 (MQ=255)                      aGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAAt                                                                                                                            >  8:266571/1‑149 (MQ=255)                      aGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAAt                                                                                                                            <  4:164148/149‑1 (MQ=255)                          cTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGt                                                                                                                        <  2:44134/149‑1 (MQ=255)                           tttGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAAt                                                                                                                            <  3:22492/144‑1 (MQ=255)                           tttGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAAt                                                                                                                            >  4:22492/1‑144 (MQ=255)                             tGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGtt                                                                                                                       <  6:228265/147‑1 (MQ=255)                             tGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTg                                                                                                                      <  7:260936/148‑1 (MQ=255)                                tCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCAc                                                                                                                  >  5:157052/1‑149 (MQ=255)                                  aGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGc                                                                                                                <  2:328289/149‑1 (MQ=255)                                      aaTGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGACTCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCg                                                                                                            >  4:221106/1‑149 (MQ=255)                                      aaTGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCg                                                                                                            <  7:211627/149‑1 (MQ=255)                                      aaTGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGT‑CAGCGAAC‑ACCAGGCATg                                                                                                                                  >  5:175711/1‑126 (MQ=255)                                      aaTGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGT‑CAGCGAAC‑ACCAGGCATg                                                                                                                                  <  6:175711/126‑1 (MQ=255)                                      aaTGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGT‑CAGCGAAC‑ACCAGGCATg                                                                                                                                  >  7:135360/1‑126 (MQ=255)                                      aaTGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGT‑CAGCGAAC‑ACCAGGCATg                                                                                                                                  <  8:135360/126‑1 (MQ=255)                                             tGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGAc                                                                                                       <  3:8088/147‑1 (MQ=255)                                                    aGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGACAAATCTGTTg                                                                                                                      >  7:332311/1‑125 (MQ=255)                                                    aGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGACAAATCTGTTg                                                                                                                      <  8:332311/125‑1 (MQ=255)                                                      cTGGTTGGTGGTGACGATCGCTTCGAGCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGttctt                                                                                            <  8:67017/149‑1 (MQ=255)                                                         gttggtGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGttcttc                                                                                           <  1:342482/147‑1 (MQ=255)                                                                gtgACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACt                                                                                   <  2:368498/148‑1 (MQ=255)                                                                    cGATCTCTTCGTTCATCTCACTGTCACTGACATTGTCTGGCAGGTATTCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCAc                                                                              <  3:269311/149‑1 (MQ=255)                                                                    cGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCa                                                                               <  8:42082/148‑1 (MQ=255)                                                                    cGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCAc                                                                              <  2:33122/149‑1 (MQ=255)                                                                      aTCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAACTACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACg                                                                             <  1:110395/149‑1 (MQ=255)                                                                              gATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGGTGTTCTTCCTGCAACTGCCACGGCAAACGc                                                                     >  2:149923/1‑148 (MQ=255)                                                                               aTCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGcaa                                                                   >  3:110075/1‑149 (MQ=255)                                                                                    cTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAAtt                                                              >  2:210540/1‑149 (MQ=255)                                                                                           tcactgACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTgg                                                            <  2:61869/144‑1 (MQ=255)                                                                                           tcactgACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTgg                                                            >  1:61869/1‑144 (MQ=255)                                                                                           tcactgACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGcc                                                       <  7:266571/149‑1 (MQ=255)                                                                                               tgACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAg                                                   >  5:174086/1‑149 (MQ=255)                                                                                               tgACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAg                                                   >  2:66091/1‑149 (MQ=255)                                                                                               tgACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCAACGGCAACCGCAACAATTGGTTGCCGTAg                                                   >  7:324024/1‑149 (MQ=255)                                                                                                 aCATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGc                                                  <  4:110075/148‑1 (MQ=255)                                                                                                           gCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTgg                                                            >  6:143877/1‑128 (MQ=255)                                                                                                           gCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTgg                                                            <  5:143877/128‑1 (MQ=255)                                                                                                           gCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCaa                                        <  6:157052/148‑1 (MQ=255)                                                                                                            cAGGTAATCGGTGGTGCCGGTGCTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAat                                      >  7:181453/1‑149 (MQ=255)                                                                                                                  aTCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCa                                <  6:174086/149‑1 (MQ=255)                                                                                                                             cGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTa                     >  6:147633/1‑149 (MQ=255)                                                                                                                                  ttGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCtt                <  1:66091/149‑1 (MQ=255)                                                                                                                                    gagaATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGc              >  2:22616/1‑149 (MQ=255)                                                                                                                                    gagaATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGc              >  7:45367/1‑149 (MQ=255)                                                                                                                                    gagaATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGc              >  3:139309/1‑149 (MQ=255)                                                                                                                                      gaATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCa            >  2:37458/1‑149 (MQ=255)                                                                                                                                      gaATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCa            >  2:43084/1‑149 (MQ=255)                                                                                                                                            gCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCACTcgcg      >  3:197893/1‑149 (MQ=255)                                                                                                                                                gTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCACTCGCGGCaa  >  6:374545/1‑149 (MQ=255)                                                                                                                                                gTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCACTCGCGGCaa  <  1:210540/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TGCCCCTTTGTCAGCAACGGCAGCGCTTTGCTCAGGCAATGCTCTGCGGCCAGCTGGTTGGTGGTGACGATCGCTTCGATCAGCTCACTGTCACTGACATTGTCTGGCAGGTAATCGGTGGTGCCGGTGTTGAGAATCAGCCCGTCCAGCGAAC‑ACCAGGCATGGCAAATCTGTTGGCACGCATCGGTGACGTGTTGTTCTTCCTGCAACTGCCACGGCAAACGCAACAATTGGTTGCCGTAGCGTGCCCCCAAATCATCCAGTGCATCGCTATCCTTGCCACTCGCGGCAA  >  NC_002947/857002‑857293

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