Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,328,512	G→C	8.0%	A21A (GCG→GCC)	PP_2927 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,328,512	0	G	C	8.0%	158.5 / 2.9	50	A21A (GCG→GCC)	PP_2927	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base C (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CGAAAACCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCG  >  NC_002947/3328368‑3328659                                                                                                                                                 |                                                                                                                                                     cgaaAACCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGc                                                                                                                                                  >  2:218012/1‑149 (MQ=255)    acaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTg‑g                                                                                                                                              <  8:39888/147‑1 (MQ=255)    aaaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTg‑g                                                                                                                                              >  4:330891/1‑149 (MQ=255)     aaCCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑Gc                                                                                                                                             >  4:356167/1‑149 (MQ=255)         cTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAAt                                                                                                                                         <  3:356167/149‑1 (MQ=255)              aCCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATTGGC‑CAGCGTCGCCGCTGCGCCCATGTcgc                                                                                                                                    <  4:344045/149‑3 (MQ=38)                  accaccGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGAAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATtctc                                                                                                                                >  4:93091/1‑149 (MQ=255)                     accGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGc                                                                                                                             <  4:245063/149‑1 (MQ=255)                      ccGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCg                                                                                                                            <  2:289241/149‑1 (MQ=255)                       cGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGa                                                                                                                           <  4:181450/149‑1 (MQ=255)                         ccTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGAt                                                                                                                          >  7:352293/1‑148 (MQ=255)                           ttCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGaa                                                                                                                       >  1:330421/1‑149 (MQ=255)                              ccGGAGCGGACGTCCCTCGTTGGCGATGCCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGAc                                                                                                                                                    >  2:342739/1‑118 (MQ=255)                              ccGGAGCGGACGTCCCTCGTTGGCGATGCCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGAc                                                                                                                                                    <  1:342739/118‑1 (MQ=255)                               cGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTgg                                                                                                                   <  3:261778/149‑1 (MQ=255)                                 gAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTgggg                                                                                                                 <  8:352293/149‑1 (MQ=255)                                            ccTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGc                                                                                                      >  2:178928/1‑149 (MQ=255)                                            ccTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCCACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGc                                                                                                      >  2:130836/1‑149 (MQ=255)                                            ccTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCCACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGc                                                                                                      >  2:135807/1‑149 (MQ=255)                                             cTCGTTGGCGTTGGACTTCTGGCTGAGGAGGGCAGGTCATGGACGACGCAATGCTTCACGATTTCATTGGCAAGCTGGTGGGCTACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGcc                                                                                                     <  6:233812/149‑1 (MQ=255)                                                gTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCAt                                                                                                   <  1:79793/148‑1 (MQ=255)                                                gTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATg                                                                                                  >  3:278976/1‑149 (MQ=255)                                                      gATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGa                                                                                             >  1:351320/1‑148 (MQ=255)                                                            cttcttGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAg                                                                                      >  6:193828/1‑149 (MQ=255)                                                             ttcttGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAgg                                                                                     >  6:177250/1‑149 (MQ=255)                                                             ttcttGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAgg                                                                                     >  5:290004/1‑149 (MQ=255)                                                                  gCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCgg                                                                                  <  7:44023/147‑1 (MQ=255)                                                                  gCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGt                                                                                 >  3:169654/1‑148 (MQ=255)                                                                          gggCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGa                                                                        <  5:177250/149‑1 (MQ=255)                                                                          gggCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGa                                                                        <  1:218012/149‑1 (MQ=255)                                                                                gTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGc                                                                                                                             <  7:336200/90‑1 (MQ=255)                                                                                gTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGc                                                                                                                             >  8:336200/1‑90 (MQ=255)                                                                                 tCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTg                                                                 >  4:25298/1‑149 (MQ=255)                                                                                      gACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCgg                                                             <  1:178928/148‑1 (MQ=255)                                                                                       aCGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGaa                                                           <  4:278976/149‑1 (MQ=255)                                                                                         gAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGaa                                                           >  6:19042/1‑147 (MQ=255)                                                                                            gCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGAc                                                      <  3:93091/149‑1 (MQ=255)                                                                                              aaGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCg                                                    >  1:184619/1‑149 (MQ=255)                                                                                                 gCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCCACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCt                                                 <  1:130836/149‑1 (MQ=255)                                                                                                 gCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCCACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCt                                                 <  1:135807/149‑1 (MQ=255)                                                                                                      aCGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCAc                                            >  8:107834/1‑149 (MQ=255)                                                                                                              aTGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCt                                    >  4:255319/1‑149 (MQ=255)                                                                                                                gggCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTgg                                  >  6:283935/1‑149 (MQ=255)                                                                                                                     aGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCga                             <  5:19042/149‑1 (MQ=255)                                                                                                                     aGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCga                             <  6:290004/149‑1 (MQ=255)                                                                                                                        tggtggGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGATCTGGGGCTGTGCCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGaa                          >  5:282620/1‑149 (MQ=255)                                                                                                                          gtggGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATg                        >  2:5594/1‑149 (MQ=255)                                                                                                                          gtggGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATg                        <  3:330891/149‑1 (MQ=255)                                                                                                                                    tGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGcc              <  3:25298/149‑1 (MQ=255)                                                                                                                                         gcAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAggc         >  8:186121/1‑149 (MQ=255)                                                                                                                                              ggcgACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCAt    >  3:147143/1‑149 (MQ=255)                                                                                                                                                cgACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAATGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTAGACTGCCCAGGCGGCATCg  >  2:273522/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                     CGAAAACCCTTGAACCGACCACCGCCTTCCCGGAGCGGACGTCCCTCGTTGGCGATGGCCTTCTTGCAGAGGAGGGCAGGTCATGGACGAGGCAAGGCTTCACGATTTCATGGGCAAGCTGGTGGGCGACATGGGCGCAGCGGCGACGCTG‑GCCAATGTCATTCTCGGCGATGAACTGGGGCTGTACCGGGCCATGGCCGACAGCCAGGCGGTCACCCCTGAGCAACTGGCGGAAAAGACCGGCTGCCACCCGCGGCTGGTGCGAGAATGGTTGAATGCCCAGGCGGCATCG  >  NC_002947/3328368‑3328659

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