Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,546,666	T→C	8.3%	C589C (TGT→TGC)	PP_2238 →	metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,546,666	0	T	C	8.3%	129.1 / 2.8	48	C589C (TGT→TGC)	PP_2238	metallopeptidase
Reads supporting (aligned to +/- strand):  ref base T (24/20);  new base C (2/2);  total (26/22)
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.85e-01

TCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCG  >  NC_002947/2546518‑2546808                                                                                                                                                      |                                                                                                                                               tCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGATCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGctgt                                                                                                                                                >  1:12841/1‑149 (MQ=255)        gatcgaGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGCCTGTTGc                                                                                                                                         >  6:42784/1‑149 (MQ=255)               ccGGGTACGTACGGCCCGGGGCTGTGGGGCGTTGGCATTGTTAACCTTGTGGTCAACGGTTCAGGGGGGTAGTTGTTCTATTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGt                                                                                                                                   <  1:282109/149‑1 (MQ=255)                    tacgtacCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGCCTGTTGCCCGAGTTGCTGa                                                                                                                             <  5:42784/149‑1 (MQ=255)                          ccGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAag                                                                                                                        <  7:166905/148‑1 (MQ=255)                             cccGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGCCCAaggag                                                                                                                     >  5:298730/1‑148 (MQ=255)                              ccGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAaggag                                                                                                                     >  1:251265/1‑147 (MQ=255)                              ccGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAaggag                                                                                                                     >  2:15241/1‑147 (MQ=255)                                 gggCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGt                                                                                                                <  1:15241/149‑1 (MQ=255)                                  ggCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTg                                                                                                               <  5:294185/149‑1 (MQ=255)                                   gCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTg                                                                                                               >  4:264068/1‑148 (MQ=255)                                    cTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGcc                                                                                                                                        <  2:329203/122‑1 (MQ=255)                                    cTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGcc                                                                                                                                        >  1:329203/1‑122 (MQ=255)                                         gggCGGTCTCATCGAGAACCTGGTGGGCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGAATTCCTGCAGTTCGAAACCTTGACGCTGAGCCCTAGCGACACCCGCTGTCTGTTTCCCGAGTTGCTGACCAAGGAGGAGGTTGAGTGGc                                                                                                        <  6:146500/149‑1 (MQ=255)                                          ggCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCt                                                                                                       <  3:264068/149‑1 (MQ=255)                                              ttCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAAc                                                                                                   <  2:251265/149‑1 (MQ=255)                                                  cATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCt                                                                                               >  1:297078/1‑149 (MQ=255)                                                  cATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGCCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGc                                                                                                >  1:176234/1‑148 (MQ=255)                                                      gagaACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTg                                                                                                          >  6:338492/1‑134 (MQ=255)                                                      gagaACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTg                                                                                                          <  5:338492/134‑1 (MQ=255)                                                         aaCCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTAc                                                                                             <  7:113349/144‑1 (MQ=255)                                                         aaCCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTAc                                                                                             >  8:113349/1‑144 (MQ=255)                                                           ccTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCt                                                                                                                               >  1:346020/1‑108 (MQ=255)                                                           ccTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCt                                                                                                                               <  2:346020/108‑1 (MQ=255)                                                                       ccGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAgc                                                                          >  2:150106/1‑149 (MQ=255)                                                                                      g‑agTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGc                                                           <  1:150106/149‑1 (MQ=255)                                                                                          tGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGCCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTcc                                                        <  2:176234/149‑1 (MQ=255)                                                                                             cTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAg                                                     >  3:294505/1‑149 (MQ=255)                                                                                                   tGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAaggaggagg                                                                                                                 <  2:3968/83‑1 (MQ=255)                                                                                                          cTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAg                                                     >  6:89970/1‑136 (MQ=255)                                                                                                          cTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAg                                                     <  5:89970/136‑1 (MQ=255)                                                                                                          cTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGc                                        >  5:36475/1‑149 (MQ=255)                                                                                                               gTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggc                                   >  4:132538/1‑149 (MQ=255)                                                                                                                ttCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggct                                  <  2:297078/149‑1 (MQ=255)                                                                                                                   gAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGga                               <  2:12841/149‑1 (MQ=255)                                                                                                                   gAAACCTTGACGCTGTGCCCTATCGACACCAGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGAtc                                               >  4:298235/1‑132 (MQ=255)                                                                                                                    aaaCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGgag                              <  3:132538/149‑1 (MQ=255)                                                                                                                     aaCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggctg                                 >  2:241802/1‑145 (MQ=255)                                                                                                                      aCCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGAt                            >  1:150311/1‑149 (MQ=255)                                                                                                                         ttGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggc                                   <  2:229656/139‑1 (MQ=255)                                                                                                                         ttGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCctggc                                   >  1:229656/1‑139 (MQ=255)                                                                                                                                     ccTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTAGCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCt             <  2:150311/149‑1 (MQ=255)                                                                                                                                      cTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTg            >  5:301181/1‑149 (MQ=255)                                                                                                                                      cTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTg            >  5:277451/1‑149 (MQ=255)                                                                                                                                             cacCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCAAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCa     <  4:294505/149‑1 (MQ=255)                                                                                                                                                ccGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCg  >  5:244557/1‑149 (MQ=255)                                                                                                                                                ccGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCg  >  4:116497/1‑149 (MQ=255)                                                                                                                                                ccGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCg  >  1:50926/1‑149 (MQ=255)                                                                                                                                                 cGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCg  >  3:164459/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                               TCAGCTCGATCGAGCCGGGTACGTACCGCCCGGGGCTGTGGGGCGTTCGCATCGAGAACCTGGTGGTCAACCGAGAAGCGGGGAAG‑AGTGCCTTTGGTGACTTCCTGCAGTTCGAAACCTTGACGCTGTGCCCTATCGACACCCGCTGTCTGTTGCCCGAGTTGCTGACCAAGGAGGAGGTGGAGTGGCTGAACGGCTACCACGCCTGCGTACATGAGCGCCTGGCGCCCTTGCTCCAAGGGGATGCGCTAGCCTGGCTGGAGATGCGTACTGCACCGCTGTAACGCAGCG  >  NC_002947/2546518‑2546808

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