Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,063,205	T→C	10.2%	intergenic (+62/‑145)	PP_1838 → / → PP_1839	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,063,205	0	T	C	10.2%	127.7 / 5.0	50	intergenic (+62/‑145)	PP_1838/PP_1839	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/22);  new base C (2/3);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

GCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAT  >  NC_002947/2063065‑2063351                                                                                                                                               |                                                                                                                                                   gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCtt                                                                                                                                            <  2:161994/149‑1 (MQ=255) gCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAAcc                                                                                                                                              >  3:259464/1‑147 (MQ=38)   tGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgt                                                                                                                                          >  7:190096/1‑149 (MQ=255)    ggAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTtgtg                                                                                                                                         >  2:93259/1‑149 (MQ=255)             tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCaag                                                                                                                                     >  3:93579/1‑142 (MQ=255)             tatCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGGGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCaag                                                                                                                                     <  4:93579/144‑3 (MQ=255)              atCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgt                                                                                                                               <  8:36224/149‑1 (MQ=255)               tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                              <  8:190096/149‑1 (MQ=255)               tCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCtgtg                                                                                                                              <  4:242233/149‑1 (MQ=255)               tCCAGACCGCGACCAAGAGCCTGCTCGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTATGACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCCCGCCCACAGGCATTGCCCAGGGTTAGGACCTTGTGCAGGTCCagt                                                                                                                               >  1:298058/1‑149 (MQ=255)                ccAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGGCCTGTGc                                                                                                                             >  5:127111/1‑149 (MQ=255)                 cAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCtgtgg                                                                                                                             <  4:259464/148‑2 (MQ=21)                    aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCCCGGGCTTTCCCGCTCCCACAGGCATTGCACCGGGTTCGGACCATGTGCAGGTCCTGTgcgcgc                                                                                                                         >  4:236611/1‑149 (MQ=255)                    aCCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCAT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGc                                                                                                                         <  1:69321/149‑1 (MQ=18)                     ccGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCg                                                                                                                        >  3:45540/1‑149 (MQ=255)                      cgcgACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCgg                                                                                                                       <  6:156739/149‑1 (MQ=18)                            cAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGc                                                                                                                 >  3:248965/1‑149 (MQ=255)                              agagCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCGCCCCCAGGCCTTGCACAGGGTTCGGGCCTTGTGCAGGTGCTGTGCGCGCGGGCCAGccc                                                                                                               >  3:146470/1‑149 (MQ=255)                                    tgctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTTCCCGCGCCCCCAGGCATTGCCCAGGGGTAGGACCATGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                                         >  3:67382/1‑149 (MQ=255)                                      ctgGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACTGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAAt                                                                                                         <  5:88268/147‑1 (MQ=255)                                            aaCACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCGCCCACAGGCATTGCACAGGGTTAGGACCTTGTGCAGGTACTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCg                                                                                                 >  8:176725/1‑149 (MQ=255)                                                       cGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTcc                                                                                      <  1:93259/149‑1 (MQ=255)                                                        gtggtgCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCCCAGGCATTGCACAGGGTTAGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                     >  8:100691/1‑149 (MQ=255)                                                        gtggtgCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTccc                                                                                     <  4:45540/149‑1 (MQ=255)                                                           gtgCACGCCACACCGAAGTA‑AACCGGGGCCGCGT‑GGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGCACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGa                                                                                  <  5:303888/149‑1 (MQ=255)                                                              cACGGCACACCGAAGTA‑AACCGGGGCCGCGTATGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACATTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACt                                                                                >  1:33803/1‑149 (MQ=255)                                                                      acCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                           >  5:81218/1‑145 (MQ=255)                                                                      acCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACa                                                                           <  6:81218/145‑1 (MQ=255)                                                                        cGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATcc                                                                     >  3:269080/1‑149 (MQ=255)                                                                         gAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATccc                                                                    <  4:248965/149‑1 (MQ=255)                                                                                     gggCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGTGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGt                                                         <  5:140122/149‑1 (MQ=255)                                                                                       gCCGCGT‑TGCGGCCCATGCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTGCGGACCTGGTCCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCa                                                       <  1:43949/149‑1 (MQ=255)                                                                                                 cGGCCCATTCGCGGGCTTGCCGGCTCCCACCGGCATGGCACAGGGTTCGGACCTTTTGCTGGGCCTGTGCTCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCTAGGTCAGTTATTCa                                               <  5:220218/148‑1 (MQ=38)                                                                                                 cGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGCGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCa                                               <  6:127111/148‑1 (MQ=255)                                                                                                     ccATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCTGGGCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCAt                                          <  7:135292/149‑1 (MQ=255)                                                                                                                cTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTgg                               <  2:33803/149‑1 (MQ=255)                                                                                                                   gCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGgtgt                            >  3:74106/1‑149 (MQ=255)                                                                                                                     ccGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTc                          >  7:157520/1‑149 (MQ=255)                                                                                                                       gCTCGCACAGGACCTGCACAAGGTCCGAACCCTGTGCAATGCCTGTGGGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGt                        >  1:328441/1‑149 (MQ=18)                                                                                                                            cacaGGCATTGCACAGGGATCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAACCCGGCGAATGGGCCCAGCACAATATTCCCCGACTCAACATTATCCCCGCCCCAACGTCAGTTATACAGGCATCACACACCTGGTGTTCGTAGACt                   >  4:103200/1‑149 (MQ=255)                                                                                                                              caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACt                   >  7:326843/1‑147 (MQ=255)                                                                                                                              caGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTga                 <  1:139141/149‑1 (MQ=255)                                                                                                                               aGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCGGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgag                <  4:269080/149‑1 (MQ=255)                                                                                                                                 gctttGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTgaga               <  3:236611/145‑1 (MQ=255)                                                                                                                                         caGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACTCCTGGTGTTCGTTGTCTGAGAGGGTCCGt       >  7:334959/1‑148 (MQ=255)                                                                                                                                          aGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGtt      >  2:25067/1‑148 (MQ=255)                                                                                                                                            ggTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCaa   >  6:156525/1‑149 (MQ=255)                                                                                                                                             gTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCaa   <  3:103200/148‑1 (MQ=255)                                                                                                                                             gTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt  <  7:176725/149‑1 (MQ=255)                                                                                                                                             gTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt  <  4:74106/149‑1 (MQ=255)                                                                                                                                              ttCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCACGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAt  >  3:344473/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                   GCTGGAAGCGTATATCCAGACCGCGACCAAGAGCCTGCTGGATAACACCGCCGGCGTGGTGCACGGCACACCGAAGTA‑AACCGGGGCCGCGT‑TGCGGCCCATTCGCGGGCTTGCCCGCTCCCACAGGCATTGCACAGGGTTCGGACCTTGTGCAGGTCCTGTGCGAGCGGGCAAGCCCGCGAATGGGCCCCGCACAATTTTCCCGGACTGAACATTATCCCCGCCCCAAGGTCAGTTATTCAGGCATCACACACCTGGTGTTCGTTGACTGAGAGGGTCCGTTCAAT  >  NC_002947/2063065‑2063351

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