Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	361,452	C→T	6.9%	T368M (ACG→ATG)	PP_0300 →	metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	361,452	0	C	T	6.9%	173.3 / 2.2	58	T368M (ACG→ATG)	PP_0300	metallopeptidase
Reads supporting (aligned to +/- strand):  ref base C (30/24);  new base T (2/2);  total (32/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 = 9.12e-01

CTTTTTCGGCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAATTTTTATCCGGA  >  NC_002947/361311‑361599                                                                                                                                              |                                                                                                                                                    cTTTTTCGGCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATg                                                                                                                                              <  8:178628/149‑1 (MQ=255)      tCGGCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCGCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAAc                                                                                                                                         >  4:231344/1‑149 (MQ=255)       cGGCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACg                                                                                                                                        >  3:356191/1‑149 (MQ=255)         gCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTc                                                                                                                                      >  8:67845/1‑149 (MQ=255)         gCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTc                                                                                                                                      <  4:356191/149‑1 (MQ=255)           aaTTAGTTGCTTTGGTTCGTTAGGCCACTAATTAGGGCATAACATGGGTGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGa                                                                                                                                    <  3:37047/149‑1 (MQ=255)             ttAGTTTCTTTGGTTCGTTAGGCCACGAATTAGGGAATAACATGGTGTCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTAACAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAAt                                                                                                                                  <  2:265258/149‑1 (MQ=255)                gTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATgg                                                                                                                                >  5:241167/1‑148 (MQ=255)                gTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGt                                                                                                                               >  7:101357/1‑149 (MQ=255)                          ttCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCt                                                                                                                     <  7:67845/149‑1 (MQ=255)                            cGTTAGGCCACTAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTa                                                                                                                    <  8:101357/148‑1 (MQ=255)                                  gCCAGGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGAGTCCCTACCCCCTTCCGCTTATGGCTACCAGCATACAGCGCGAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCGTATCATTc                                                                                                              >  8:113214/1‑148 (MQ=255)                                   ccACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGTCTACCAGCATACAGCTCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCaa                                                                                                            <  5:271939/149‑1 (MQ=255)                                    cACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCaaa                                                                                                           >  1:271712/1‑149 (MQ=255)                                                     cATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGc                                                                                                                             >  4:368708/1‑114 (MQ=255)                                                     cATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGc                                                                                                                             <  3:368708/114‑1 (MQ=255)                                                     cATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTAt                                                                                           >  1:97404/1‑148 (MQ=255)                                                           ggCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAgggg                                                                                    <  2:271712/149‑1 (MQ=255)                                                            gCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCCTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGt                                                                                   >  8:202172/1‑149 (MQ=255)                                                               ggTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTACCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTcc                                                                                <  6:241167/149‑1 (MQ=255)                                                                           gAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGc                                                                                                >  5:55461/1‑121 (MQ=255)                                                                           gAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGc                                                                                                <  6:55461/121‑1 (MQ=255)                                                                           gAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACgg                                                                     >  2:192303/1‑148 (MQ=255)                                                                              gCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACgg                                                                     >  7:225993/1‑145 (MQ=255)                                                                               ttCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGAt                                                                <  3:251575/148‑1 (MQ=255)                                                                               cTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAg                                                                                       <  3:45775/126‑1 (MQ=255)                                                                               cTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAg                                                                                       >  4:45775/1‑126 (MQ=255)                                                                                 ccAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTa                                                                        >  4:136647/1‑139 (MQ=255)                                                                                 ccAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTa                                                                        <  3:136647/139‑1 (MQ=255)                                                                                 ccAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGAt                                                                >  4:251575/1‑147 (MQ=255)                                                                                 ccAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATa                                                               >  3:101535/1‑148 (MQ=255)                                                                                  cAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATaaa                                                             <  4:101535/149‑1 (MQ=255)                                                                                       gCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTAt                                                                             >  4:107774/1‑128 (MQ=255)                                                                                            ccTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTgg                                                   >  3:28854/1‑149 (MQ=255)                                                                                             cTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATaaaa                                                            <  1:144579/139‑1 (MQ=255)                                                                                             cTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATaaaa                                                            >  2:144579/1‑139 (MQ=255)                                                                                                     ttACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCTGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACtt                                          >  2:53720/1‑149 (MQ=255)                                                                                                         gCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCaa                                      <  3:231344/149‑1 (MQ=255)                                                                                                         gCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCaa                                      <  2:97404/149‑1 (MQ=255)                                                                                                           ttATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTgg                                                   <  3:57418/134‑1 (MQ=255)                                                                                                           ttATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTgg                                                   >  4:57418/1‑134 (MQ=255)                                                                                                             aTGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGa                                   >  3:16898/1‑148 (MQ=255)                                                                                                             aTGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGa                                   >  6:184540/1‑148 (MQ=255)                                                                                                              tGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAatgat                                                       >  6:144432/1‑127 (MQ=255)                                                                                                              tGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAatgat                                                       <  5:144432/127‑1 (MQ=255)                                                                                                               ggCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGaaaa                                >  7:32727/1‑149 (MQ=255)                                                                                                                    ccAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCg                             <  4:16898/147‑1 (MQ=255)                                                                                                                             aGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCTGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTg                  <  1:53720/149‑1 (MQ=255)                                                                                                                               cgcCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTAACAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGc                 >  4:320504/1‑148 (MQ=255)                                                                                                                                   aaaTTTCCATATGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCaa                                      <  7:123797/123‑1 (MQ=255)                                                                                                                                   aaaTTTCCATATGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCaa                                      >  8:123797/1‑123 (MQ=255)                                                                                                                                   aaaTTTCCATATGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCaa                                      >  6:379886/1‑123 (MQ=255)                                                                                                                                   aaaTTTCCATATGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCaa                                      <  5:379886/123‑1 (MQ=255)                                                                                                                                    aaTTTCCATACGCACATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAAttt           <  2:71810/149‑1 (MQ=255)                                                                                                                                        tCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAATTTTTa        >  6:285097/1‑148 (MQ=255)                                                                                                                                            tACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAATTTTTATCCgg   <  5:184540/149‑1 (MQ=255)                                                                                                                                            tACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAATTTTTATCCg    >  4:70638/1‑148 (MQ=255)                                                                                                                                             aCGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAATTTTTATCCGGa  >  3:103416/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                    CTTTTTCGGCAATTAGTTGCTTTGGTTCGTTAGGCCACGAATTAGGGCATAACATGGGGGCCGGTCATAAGGACGAGGCTCCAGTTGCGTCCCTACCCCCTTACGCTTATGGCTACCAGCATACAGCGCCAAATTTCCATACGCAGATGCGAACGTCGAATGGTGCCATTCCCTATCATTCAAACCCGCGGTTGCAGTATCAAGGGGTTCCTATGGGTACGGTAGATAAAAATGATGTGGCGCGGACTTTCAATGAAAATCGGGACACTGTTGCCAATTTTTATCCGGA  >  NC_002947/361311‑361599

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