Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,985,941	G→T	7.7%	C73F (TGC→TTC)	oplB →	5‑oxoprolinase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,985,941	0	G	T	7.7%	162.3 / 3.4	52	C73F (TGC→TTC)	oplB	5‑oxoprolinase B
Reads supporting (aligned to +/- strand):  ref base G (23/25);  new base T (2/2);  total (25/27)
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.23e-01

TGGAAAAGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAA  >  NC_002947/3985803‑3986088                                                                                                                                            |                                                                                                                                                    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                                                                                                                                    >  3:112053/1‑149 (MQ=255)        gACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCg                                                                                                                                    >  5:238494/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                                                                                                                                  >  5:121753/1‑132 (MQ=255)                           ttCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGc                                                                                                                                  <  6:121753/132‑1 (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)                                                                                                                                         aGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTg     <  3:111877/148‑1 (MQ=255)                                                                                                                                            gCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCaa  <  5:26271/148‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                    TGGAAAAGACCGGGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGA‑AGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAA  >  NC_002947/3985803‑3986088

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