Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,201,878	C→T	11.7%	A46V (GCC→GTC)	PP_4583 →	putative Peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,201,878	0	C	T	11.7%	98.4 / 6.0	43	A46V (GCC→GTC)	PP_4583	putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (16/22);  new base T (2/3);  total (18/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.42e-01

AAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGC  >  NC_002947/5201730‑5202026                                                                                                                                                     |                                                                                                                                                     aaGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGcc                                                                                                                                                     >  2:439625/1‑150 (MQ=255)      aCCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGc                                                                                                                                                      >  1:115658/1‑144 (MQ=255)      aCCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGc                                                                                                                                                      <  2:115658/144‑1 (MQ=255)       ccTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAg                                                                                                                                               <  1:439625/150‑1 (MQ=255)        cTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAgg                                                                                                                                              >  2:797429/1‑150 (MQ=255)         ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGc                                                                                                                                             >  2:73921/1‑150 (MQ=255)         ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGc                                                                                                                                             >  1:580745/1‑150 (MQ=255)           cAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCag                                                                                                                                           <  2:582133/150‑2 (MQ=255)                         ccccAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAgcgc                                                                                                                             <  1:421823/150‑1 (MQ=255)                               cccccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTc                                                                                                                       <  2:550746/150‑1 (MQ=255)                                  ccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCgaa                                                                                                                    <  2:580745/150‑1 (MQ=255)                                         cccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCa                                                                                                             >  1:179533/1‑150 (MQ=255)                                          ccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACTCCCGAGGTACTCGCTTACCCGCGAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCaa                                                                                                            <  1:274513/150‑1 (MQ=255)                                            aaGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCa                                                                                                                                            <  2:181594/116‑1 (MQ=255)                                            aaGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGTCGGCCAGGCa                                                                                                                                            >  1:181594/1‑116 (MQ=255)                                                      ggCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCt                                                                                                <  1:208496/150‑1 (MQ=255)                                                       gCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTg                                                                                               <  1:73921/150‑1 (MQ=255)                                                         cACCGACCCGTATGCCTGTCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGga                                                                                             <  2:381350/150‑1 (MQ=255)                                                         cACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGga                                                                                             <  1:608387/150‑1 (MQ=255)                                                         cACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCAGCCTGga                                                                                             <  2:488125/150‑1 (MQ=255)                                                           ccGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGgaga                                                                                           >  1:321703/1‑150 (MQ=255)                                                                      gCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTcc                                                                                <  2:640695/150‑2 (MQ=255)                                                                      gCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCa                                                                                <  2:324446/150‑1 (MQ=255)                                                                        cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAtt                                                                              >  2:787242/1‑150 (MQ=255)                                                                          ggCTGCAACAGCGCGACACCCCCGGGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGc                                                                            <  2:588758/150‑1 (MQ=255)                                                                            cTGCAACAGCGCGACCCCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCg                                                                          <  1:415127/150‑1 (MQ=255)                                                                                aaCAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGgcgc                                                                      >  1:491735/1‑150 (MQ=255)                                                                                        gACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCTTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGc                                                              >  2:379184/1‑150 (MQ=255)                                                                                             ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTggg                                                                >  1:237252/1‑143 (MQ=255)                                                                                             ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTggg                                                                <  2:237252/143‑1 (MQ=255)                                                                                              cccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTAc                                                        >  1:169183/1‑150 (MQ=255)                                                                                                         cGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTacac                                             <  1:739755/150‑1 (MQ=255)                                                                                                              aCCTGCAAGCTGAAAACGCCTACCAGGCAGCGTAGAAGTCCCTCCCGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGCCCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTa                                        <  2:285658/150‑1 (MQ=255)                                                                                                                cTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTacc                                      >  2:687774/1‑150 (MQ=255)                                                                                                                     aGCTGACAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAAGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCAccgc                                 <  2:360071/150‑1 (MQ=255)                                                                                                                          aaaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcg                            >  1:632234/1‑150 (MQ=255)                                                                                                                           aaaCGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcga                           >  1:697786/1‑150 (MQ=255)                                                                                                                                     cAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcgacga                        <  1:71137/143‑1 (MQ=255)                                                                                                                                     cAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcgacga                        >  2:71137/1‑143 (MQ=255)                                                                                                                                        gAAGCCTGCCTGGCCGACCAGGCAGCGTTTCGTGAGCGCCTGGTCGAAGAAAGCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGt              >  1:447903/1‑150 (MQ=255)                                                                                                                                        gAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGt              <  2:179533/150‑1 (MQ=255)                                                                                                                                        gAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGt              >  1:687595/1‑150 (MQ=255)                                                                                                                                             gtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGta                                                   >  2:595613/2‑108 (MQ=255)                                                                                                                                             gtgcctgGTCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGta                                                   <  1:595613/107‑1 (MQ=255)                                                                                                                                              tgcctgGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTAc        <  1:3148/150‑1 (MQ=255)                                                                                                                                                    gCCGACCAGGCACCGTTGCGTGAGCGACTGTTCGAACAAATCATGGGCCGCATCCTGGAGTCCGAAAATGGGGCGGGGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGc  <  1:787242/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGC  >  NC_002947/5201730‑5202026

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑