Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,985,937	A→C	7.8%	K72Q (AAG→CAG) ‡	oplB →	5‑oxoprolinase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,985,937	0	A	C	7.8%	159.9 / 3.3	51	K72Q (AAG→CAG) ‡	oplB	5‑oxoprolinase B
Reads supporting (aligned to +/- strand):  ref base A (23/24);  new base C (2/2);  total (25/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.97e-01

ATGACCCTGACCCTGGAAAAGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTT  >  NC_002947/3985790‑3986083                                                                                                                                                    |                                                                                                                                                    aTGACCCTGACCCTGGAAAAGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGa‑a                                                                                                                                                   <  5:69831/149‑1 (MQ=255)              tgGAAAAGACGGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCTCGGCCACGC‑ACTTCATGGCCg                                                                                                                                        <  4:246594/147‑1 (MQ=255)              tgGAAAAGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCTCGGCCACGC‑ACTTCATGGCCg                                                                                                                                        >  3:246594/1‑147 (MQ=255)                  aaaGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGttct                                                                                                                                  >  7:142813/1‑149 (MQ=255)                     gACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCg                                                                                                                               >  5:238494/1‑149 (MQ=255)                     gACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCg                                                                                                                               >  3:112053/1‑149 (MQ=255)                       ccGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGc                                                                                                                             >  5:259739/1‑149 (MQ=255)                       ccGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGc                                                                                                                             <  2:97666/149‑1 (MQ=255)                                   gtgtgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCACGC‑ACTTCATGGCCggg                                                                                                                                      <  8:123983/128‑3 (MQ=255)                                   gtgtgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCAACCTGGGCTCGCTGATCCCGGCCACGC‑ACTTCATGGCCggg                                                                                                                                      >  7:123983/1‑126 (MQ=255)                                      tgtTCAACCTCGCCCACGACTATTCCAACTCACTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTCGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                              <  5:222004/149‑1 (MQ=255)                                       gtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAg                                                                                                             >  3:56243/1‑149 (MQ=255)                                        ttCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGc                                                                                                                             <  6:121753/132‑1 (MQ=255)                                        ttCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGc                                                                                                                             >  5:121753/1‑132 (MQ=255)                                        ttCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAg                                                                                                             <  1:35909/148‑1 (MQ=255)                                                 gCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGAt                                                                                                   <  6:259739/149‑1 (MQ=255)                                                       gACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGAAATCGCCGAAGGCGACGTGATCTACca                                                                                             >  2:177642/1‑149 (MQ=255)                                                       gACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGGGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACca                                                                                             <  1:242026/149‑1 (MQ=255)                                                            ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                         >  6:45935/1‑148 (MQ=255)                                                               cAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAcc                                                                                     >  6:234215/1‑149 (MQ=255)                                                               cAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAc                                                                                      <  5:102744/148‑1 (MQ=255)                                                                            gACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAggg                                                                        <  5:45935/149‑1 (MQ=255)                                                                              ccACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGaa                                                                                                              <  1:271276/109‑1 (MQ=255)                                                                              ccACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGaa                                                                                                              >  2:271276/1‑109 (MQ=255)                                                                                aCTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGc                                                                    <  6:238494/149‑1 (MQ=255)                                                                                 cTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGAAATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGcc                                                                   <  1:177642/149‑1 (MQ=255)                                                                                        gAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtct                                                            >  5:108647/1‑149 (MQ=255)                                                                                         aaaTGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATtctc                                                           <  8:142813/149‑1 (MQ=255)                                                                                           aTGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGa                                                         >  4:111877/1‑149 (MQ=255)                                                                                               tCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                         <  1:71670/113‑1 (MQ=255)                                                                                               tCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                         >  2:71670/1‑113 (MQ=255)                                                                                                 cTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGActgct                                                    <  5:234215/148‑1 (MQ=255)                                                                                                  tGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGa                                                         <  2:257901/142‑1 (MQ=255)                                                                                                  tGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGa                                                         >  1:257901/1‑142 (MQ=255)                                                                                                    cAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACtgctgc                                                  <  4:112053/147‑1 (MQ=255)                                                                                                    cAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCAt                                                <  6:108647/149‑1 (MQ=255)                                                                                                        gCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAGAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGActgct                                                    >  1:191619/1‑142 (MQ=255)                                                                                                        gCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGActgct                                                    <  2:191619/141‑1 (MQ=255)                                                                                                           aGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACaaa                                         >  4:172194/1‑149 (MQ=255)                                                                                                              aCATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCg                                      <  2:22039/149‑1 (MQ=255)                                                                                                               cATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCTGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCgg                                     <  3:240386/149‑1 (MQ=255)                                                                                                                aTTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgt                                    <  3:172194/149‑1 (MQ=255)                                                                                                                     cATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTAt                                >  4:31119/1‑148 (MQ=255)                                                                                                                      aTCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTAtt                               >  6:26271/1‑148 (MQ=255)                                                                                                                      aTCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTa                              >  7:153801/1‑149 (MQ=255)                                                                                                                              gggCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGAc                       >  1:203847/1‑148 (MQ=255)                                                                                                                               ggCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGa                     <  8:153801/149‑1 (MQ=255)                                                                                                                                     cTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGgt               >  5:114583/1‑149 (MQ=255)                                                                                                                                            cGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTg         >  6:96791/1‑148 (MQ=255)                                                                                                                                             ggCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGa       <  4:56243/149‑1 (MQ=255)                                                                                                                                                 aTGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGt   >  2:3949/1‑149 (MQ=255)                                                                                                                                                  tGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGtt  <  5:96791/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATGACCCTGACCCTGGAAAAGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTT  >  NC_002947/3985790‑3986083

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