Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,985,965	A→T	11.1%	D81V (GAC→GTC)	oplB →	5‑oxoprolinase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,985,965	0	A	T	11.1%	120.1 / 5.4	45	D81V (GAC→GTC)	oplB	5‑oxoprolinase B
Reads supporting (aligned to +/- strand):  ref base A (20/20);  new base T (2/3);  total (22/23)
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.78e-01

TCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGG  >  NC_002947/3985820‑3986112                                                                                                                                                  |                                                                                                                                                    tCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACAt                                                                                                                                                  >  4:273540/1‑149 (MQ=255)   cAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATc                                                                                                                                                 >  1:156001/1‑148 (MQ=255)   cAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATc                                                                                                                                                 <  2:156001/148‑1 (MQ=255)        tgtTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                           <  4:248717/149‑1 (MQ=255)             aaCCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGa                                                                                                                                      <  3:273540/149‑1 (MQ=255)                    accACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATc                                                                                                                               <  6:180726/148‑1 (MQ=255)                     ccACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATc                                                                                                                               >  1:189310/1‑148 (MQ=255)                       aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                           >  1:294387/1‑134 (MQ=255)                       aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                           <  2:294387/134‑1 (MQ=255)                       aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTAc                                                                                                                            >  7:246229/1‑149 (MQ=255)                              ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGaa                                                                                                                                           >  5:55246/1‑127 (MQ=255)                              ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGaa                                                                                                                                           <  6:55246/127‑1 (MQ=255)                              ttCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGCACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAAc                                                                                                                      <  1:50947/148‑1 (MQ=255)                                 cAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAc                                                                                                                   <  8:246229/148‑1 (MQ=255)                                  aaCTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGAccc                                                                                                                 >  4:112599/1‑149 (MQ=255)                                      cGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGcc                                                                                                             >  7:255133/1‑149 (MQ=255)                                          gTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTag                                                                                                                            >  3:8323/1‑129 (MQ=255)                                          gTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTAc                                                                                                                            <  4:8323/130‑1 (MQ=255)                                             cGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAg                                                                                                       <  2:54490/148‑1 (MQ=255)                                                  aCTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGc                                                                                                 <  2:189310/149‑1 (MQ=255)                                                              tGATCCGGCAGGGCCAGGACATGCCCATCACCCTGGGCGCGCTTATCCCGGCCATTAAGTTCGTGGCCGAGTTATTCGGCGACGACATCGCCGAAGGCGACGTTATATACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGAc                                                                                     <  6:270457/149‑1 (MQ=255)                                                                         ggCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGCCCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTa                                                                          >  4:235103/1‑149 (MQ=255)                                                                                aCATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCg                                                                   <  2:300395/149‑1 (MQ=255)                                                                                       cATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTAtt                                                            <  3:112599/149‑1 (MQ=255)                                                                                            aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACca                                                                                                                          <  3:237332/82‑1 (MQ=255)                                                                                            aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACca                                                                                                                          >  4:237332/1‑82 (MQ=255)                                                                                            aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAg                                                       >  5:214124/1‑149 (MQ=255)                                                                                             ccTGGGCTCGCTGATCCCGGCCATGAAGTGGGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGTCGACGTGATCTACCACAACGACCCAGCCTACAAGGGCAGACATATTCTCGACTGATGCATGTACAAACCGGTGTATTACCAg                                                       >  1:26565/1‑148 (MQ=255)                                                                                             ccTGGCCTGCCTGATCCCGGCCATGACCGGCGGGGCCGAGTTCTTCGGCCACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAgg                                                      <  8:103747/149‑1 (MQ=255)                                                                                                 ggCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGAc                                                    >  1:44078/1‑147 (MQ=255)                                                                                                    tCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgtgt                                                               >  4:202805/1‑133 (MQ=255)                                                                                                    tCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGgtgt                                                               <  3:202805/133‑1 (MQ=255)                                                                                                      gCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTgg                                             <  8:255133/149‑1 (MQ=255)                                                                                                                  cATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACg                                                                                                                                    >  5:331153/1‑50 (MQ=255)                                                                                                                  cATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACg                                                                                                                                    <  6:331153/50‑1 (MQ=255)                                                                                                                       aGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGgt                                            >  7:178004/1‑133 (MQ=255)                                                                                                                       aGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGgt                                            <  8:178004/133‑1 (MQ=255)                                                                                                                           cGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAg                         <  2:44078/148‑1 (MQ=255)                                                                                                                                 cGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCAcc                  >  5:177283/1‑149 (MQ=255)                                                                                                                                     ttcttcGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGAc              >  7:44884/1‑149 (MQ=255)                                                                                                                                        ttcGGCGACGACATCGCCGAAGGCGACGTGATCTACAACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTATGGACCGTTTGAAAGGCCCACCTGAcct            >  8:264183/1‑147 (MQ=255)                                                                                                                                          cGGCGATGTCGTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAg                                                 <  1:119165/109‑1 (MQ=255)                                                                                                                                          cGGCGATGTCGTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAg                                                 >  2:119165/1‑109 (MQ=255)                                                                                                                                           ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTg                                              >  6:101546/1‑111 (MQ=255)                                                                                                                                           ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTg                                              <  5:101546/111‑1 (MQ=255)                                                                                                                                                 gacATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATcggcgg  >  7:300704/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                    TCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCGACATCGGCGG  >  NC_002947/3985820‑3986112

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑