Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,001,815	C→G	8.1%	A180G (GCC→GGC)	PP_2625 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,001,815	0	C	G	8.1%	146.7 / 2.3	48	A180G (GCC→GGC)	PP_2625	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base G (2/2);  total (24/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.66e-01

ATCAGTTGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACA  >  NC_002947/3001667‑3001960                                                                                                                                                      |                                                                                                                                                  aTCAGTTGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGc                                                                                                                                                   <  8:42657/149‑1 (MQ=255)         cATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCcgggcggg                                                                                                                                           <  2:292056/149‑1 (MQ=255)          aTACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGTTAGGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCgggg                                                                                                                                          <  4:52773/149‑1 (MQ=255)           tACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGcgggcaggcc                                                                                                                                         >  2:45331/1‑147 (MQ=255)                        ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCa                                                                                                                               <  5:165413/146‑1 (MQ=255)                         cGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                           <  6:341440/149‑1 (MQ=255)                               cAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGAt                                                                                                                     <  5:342192/149‑1 (MQ=255)                               cAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGAt                                                                                                                     <  5:336408/149‑1 (MQ=255)                                 gCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCacttgg                                                                                                                   <  1:45331/149‑6 (MQ=38)                                     aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAgg                                                                                                                 >  2:104839/1‑147 (MQ=255)                                     aTTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                               <  8:67077/149‑1 (MQ=255)                                         ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCt                                                                                                            >  8:208899/1‑148 (MQ=255)                                         ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTg                                                                                                           <  1:104839/149‑1 (MQ=255)                                                        aGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCg                                                                                             >  4:336249/1‑148 (MQ=255)                                                              gCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                      <  7:208899/149‑1 (MQ=255)                                                              gCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                      >  8:211382/1‑149 (MQ=255)                                                               cTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGTCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGACCAGAACAGGAAGCTGATCGGTGCCCAGCGAATCGAgg                                                                                     >  2:260465/1‑149 (MQ=255)                                                                 gggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATTCCGGGCGGGGACTGCCCGACATCGAGCAGATCAGGAAGCTGATCGGTGCACAGCGTATCGAGGtt                                                                                   >  5:230774/1‑149 (MQ=255)                                                                   gTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGTACGTCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgt                                                                                 <  1:260465/149‑1 (MQ=255)                                                                       aGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgagt                                                                             >  5:246805/1‑149 (MQ=255)                                                                         gTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgagtga                                                                           >  8:303665/1‑149 (MQ=255)                                                                            cGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATc                                                                         <  1:152132/148‑1 (MQ=255)                                                                               cTGAGCGGGGAGGGCTGGTCTACCCGTGAGAGCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACg                                                                      >  1:38712/1‑149 (MQ=255)                                                                                 gAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCTCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGc                                                                    >  7:344190/1‑148 (MQ=255)                                                                                  aGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGa                                                                  <  6:246805/149‑1 (MQ=255)                                                                                    cGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAgg                                                                <  8:30392/149‑1 (MQ=255)                                                                                       ggAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                               >  8:60885/1‑99 (MQ=255)                                                                                       ggAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                               <  7:60885/99‑1 (MQ=255)                                                                                                    aCCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCtt                                                >  1:203478/1‑149 (MQ=255)                                                                                                      ccGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATc                                                                                                                             >  6:316653/1‑70 (MQ=21)                                                                                                      ccGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATc                                                                                                                             <  5:316653/70‑1 (MQ=21)                                                                                                       cGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGa                                             >  4:212766/1‑149 (MQ=255)                                                                                                             a‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATgcgc                                        >  5:292251/1‑148 (MQ=255)                                                                                                                     ggCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                >  3:334317/1‑149 (MQ=255)                                                                                                                        aaGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGcc                             <  8:167586/149‑1 (MQ=255)                                                                                                                         aGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGccc                            >  4:185678/1‑149 (MQ=255)                                                                                                                          gTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGcccc                           <  8:344190/149‑1 (MQ=255)                                                                                                                               ggAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGa                      >  2:98386/1‑149 (MQ=255)                                                                                                                                 agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGAt                     >  3:77266/1‑148 (MQ=255)                                                                                                                                  gaACGGCCAAGTGCTCGATGCCGGGCGGGACCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGAt                     >  1:28818/1‑147 (MQ=255)                                                                                                                                    aCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGc                 <  2:38712/149‑1 (MQ=255)                                                                                                                                        ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCggg             >  4:313530/1‑149 (MQ=255)                                                                                                                                           aGTGCTCGATGCCGGGCGGTGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCACAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGTCCCCCTTGATCAGCCGGgcgc          <  6:230774/149‑1 (MQ=255)                                                                                                                                              gCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGcgcg         >  2:135661/1‑147 (MQ=255)                                                                                                                                              gCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGaa       <  7:192536/149‑1 (MQ=255)                                                                                                                                                 cGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca    >  8:273986/1‑149 (MQ=255)                                                                                                                                                  gATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGca    <  3:336249/148‑1 (MQ=255)                                                                                                                                                  gATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGc     <  6:114735/147‑1 (MQ=255)                                                                                                                                                   aTGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca  <  3:333285/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ATCAGTTGCATACAGGCCATCTCCCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGA‑CGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACA  >  NC_002947/3001667‑3001960

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