Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,985,963	C→T	11.1%	D80D (GAC→GAT)	oplB →	5‑oxoprolinase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,985,963	0	C	T	11.1%	116.6 / 5.2	45	D80D (GAC→GAT)	oplB	5‑oxoprolinase B
Reads supporting (aligned to +/- strand):  ref base C (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 = 2.01e-01

GGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCG  >  NC_002947/3985815‑3986102                                                                                                                                                     |                                                                                                                                            ggCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGcgac                                                                                                                                             <  3:96561/149‑1 (MQ=255)      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                                                                                                                                 <  2:294387/134‑1 (MQ=255)                            aCGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGATGTCGTCGCCGaa                                                                                                                                 >  1:294387/1‑134 (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                                                                                                                >  4:237332/1‑82 (MQ=255)                                                                                                 aCCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACca                                                                                                                <  3:237332/82‑1 (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                                       >  2:119165/1‑109 (MQ=255)                                                                                                                                               cGGCGATGTCGTCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAg                                       <  1:119165/109‑1 (MQ=255)                                                                                                                                                ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTg                                    >  6:101546/1‑111 (MQ=255)                                                                                                                                                ggCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTg                                    <  5:101546/111‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                            GGCGCTCCAGTGTGTTCAACCTCGCCCACGACTATTCCAACTCGCTGTTCGACCACTTGCCGGAAATGATCCTGCAGGGCCAGGACATTCCCATCCACCTGGGCTCGCTGATCCCGGCCATGAAGTGCGTGGCCGAGTTCTTCGGCGACGACATCGCCGAAGGCGACGTGATCTACCACAACGACCCCGCCTACAAGGGCAGCCATATTCTCGACTGCTGCATGTACAAACCGGTGTATTACCAGGACGAGCTGGTGTTCTGGACCGTTTGCAAGGGCCACCTGACCG  >  NC_002947/3985815‑3986102

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