Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	718,854	G→T	9.5%	Q355K (CAG→AAG)	PP_0613 ←	amidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	718,854	0	G	T	9.5%	117.6 / 2.2	41	Q355K (CAG→AAG)	PP_0613	amidase family protein
Reads supporting (aligned to +/- strand):  ref base G (22/15);  new base T (2/2);  total (24/17)
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.29e-01

TCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTCC  >  NC_002947/718706‑718993                                                                                                                                                     |                                                                                                                                            tccatccaTTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTg                                                                                                                                             >  1:66949/1‑149 (MQ=255)   catccatTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTgg                                                                                                                                            <  3:222565/148‑1 (MQ=255)   catccatTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCACAATCAATTCTGCCCGTGCCTCAAGCGCCTg                                                                                                                                             <  7:315049/147‑1 (MQ=255)   catccatTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGAGGAAAGGCATTCATCTTCTAAATCAATTCTGCCCGTGCCTCAAGCGCCTgg                                                                                                                                            <  7:82229/148‑1 (MQ=255)               cATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGcc                                                                                                                                 >  1:237138/1‑147 (MQ=255)                aTCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAg                                                                                                                              <  6:386577/149‑1 (MQ=255)                   cccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTgg                                                                                                                                            <  6:95661/132‑1 (MQ=255)                   cccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTgg                                                                                                                                            >  5:95661/1‑132 (MQ=255)                   cccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGc                                                                                                                            >  8:292292/1‑148 (MQ=255)                   cccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGc                                                                                                                            <  7:292292/148‑1 (MQ=255)                    ccGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGcgc                                                                                                                          >  7:246721/1‑147 (MQ=255)                       aCCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTTGGCGAATACACCCAGGcgcttg                                                                                                                       <  8:246721/149‑6 (MQ=255)                        cccAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCt                                                                                                                      <  7:313367/149‑1 (MQ=255)                              cggcACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCAc                                                                                                                <  2:237138/149‑1 (MQ=255)                                         gtCCGGCTTCGAACGCCCCCAACGGCAGCATCGGCAATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGc                                                                                                      >  1:255267/1‑148 (MQ=255)                                         gtCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGc                                                                                                      >  7:316699/1‑148 (MQ=255)                                          tCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCATGTTCGGCa                                                                                                     >  3:37390/1‑148 (MQ=255)                                            cGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATc                                                                                                  >  5:261136/1‑149 (MQ=255)                                              gCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCa                                                                                                <  5:21163/149‑1 (MQ=255)                                                     aCGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGc                                                                                         <  7:6759/149‑1 (MQ=255)                                                         caccaACGGCAGCATCGGCGATACCCGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCg                                                                                     >  3:331246/1‑149 (MQ=255)                                                            caACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGt                                                                                  <  8:33903/149‑1 (MQ=255)                                                               cGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCaa                                                                               <  2:66949/149‑1 (MQ=255)                                                               cGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCaa                                                                               <  8:316699/149‑1 (MQ=255)                                                               cGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAAATCTCCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCaa                                                                               <  7:65416/149‑1 (MQ=255)                                                                       tCGGCGCTACCAGCACGTCATATCGCTGATGAAAGGCATTCATGGTCGCAATCAATTCTGCGCGTGCCTGAAGAGCCTGGGGGTAGTCGCCAAGGCTGATCTGTGCAGGTGGTTCGGACATCCAGCGCAGCGCCGGGTCCAggacctct                                                                       >  6:109841/1‑142 (MQ=255)                                                                       tCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCt                                                                       <  6:261136/149‑1 (MQ=255)                                                                            gATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTg                                                                  <  2:171911/149‑1 (MQ=255)                                                                                ccAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCa                                                              >  3:113452/1‑149 (MQ=255)                                                                                   gCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCgg                                                                                    >  6:263661/1‑124 (MQ=255)                                                                                           tatCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAgg                                                    >  1:145684/1‑148 (MQ=255)                                                                                                                 tCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGcgccgcgccg                             >  8:79769/1‑149 (MQ=255)                                                                                                                     cGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCa                          >  5:3989/1‑148 (MQ=255)                                                                                                                        aatcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCa                                                                            >  4:122369/1‑95 (MQ=255)                                                                                                                        aatcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCa                                                                            <  3:122369/95‑1 (MQ=255)                                                                                                                        aatcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCCCCGCGCCGGCAAAc                       >  3:28170/1‑148 (MQ=255)                                                                                                                         atcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACca                     >  1:8957/1‑149 (MQ=255)                                                                                                                          tcaatTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACcac                    >  6:35690/1‑149 (MQ=255)                                                                                                                             atTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAgg                 <  6:3989/149‑1 (MQ=255)                                                                                                                              ttCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAggg                >  5:140610/1‑149 (MQ=255)                                                                                                                                cTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGgtg              >  7:218600/1‑149 (MQ=255)                                                                                                                                 tGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGgtgt             >  3:170153/1‑149 (MQ=255)                                                                                                                                  gCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGtt            >  3:141570/1‑149 (MQ=255)                                                                                                                                   cccGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGgtg              <  2:138621/146‑1 (MQ=255)                                                                                                                                     cGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTg           >  3:251241/1‑147 (MQ=255)                                                                                                                                     cGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTg           >  1:95059/1‑147 (MQ=255)                                                                                                                                             aaGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTcc  <  4:28170/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                            TCCATCCATTGCGCCATCCCCGACCCAGGCGGCACGTTGTGTCCGGCTTCGAACGCCACCAACGGCAGCATCGGCGATACCAGCACGTCATATCGCTGATGAAAGGCATTCATCTTCGCAATCAATTCTGCCCGTGCCTCAAGCGCCTGGGTGTATTCGCCAAGGCTGATCTGTGCACCTTGTTCGGCAATCCAGCGCAGCCCCGGGTCAAGCAACGCTTTCTGCTCATCGCTCAAGGCACTGGCCAACCGCGCCGCGCCGGCAAACCACAGGGTGTTGAAGGTTTCC  >  NC_002947/718706‑718993

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