Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	295,583	G→T	8.2%	L49L (CTG→CTT)	PP_0241 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	295,583	0	G	T	8.2%	155.3 / 3.9	49	L49L (CTG→CTT)	PP_0241	transporter
Reads supporting (aligned to +/- strand):  ref base G (19/26);  new base T (2/2);  total (21/28)
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.96e-01

ATGACTATCAAAGCCATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCCACAGG  >  NC_002947/295437‑295713                                                                                                                                                   |                                                                                                                                    aTGACTATCAAAGCCATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTAGGTGATTACCACCCGTTCCGTGCGGGAGCTg                                                                                                                                     >  5:146487/1‑147 (MQ=255)   gACTATCAAAGCCATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTg                                                                                                                                     <  6:146975/145‑1 (MQ=255)   gACTATCAAAGCCATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTg                                                                                                                                     >  5:146975/1‑145 (MQ=255)             gCCATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTggg                                                                                                                       <  2:276100/149‑1 (MQ=255)               cATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCa                                                                                                                     <  3:151767/149‑1 (MQ=255)                aTCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCag                                                                                                                    >  7:247583/1‑149 (MQ=255)                    aCGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCagtgaa                                                                                                                <  6:319238/149‑1 (MQ=255)                     cGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTgg                                                                                                                        >  6:19613/1‑140 (MQ=255)                     cGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTgg                                                                                                                        <  5:19613/140‑1 (MQ=255)                     cGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCagtgaag                                                                                                               >  3:133417/1‑149 (MQ=255)                      gTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCagtgaagt                                                                                                              <  6:146487/149‑1 (MQ=255)                           cAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTg                                                                                                         <  5:165139/149‑1 (MQ=255)                            aaCCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGTAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGc                                                                                                        <  2:258840/149‑1 (MQ=255)                                aGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCttt                                                                                                    <  5:53398/149‑1 (MQ=255)                                aGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCttt                                                                                                    <  8:247583/149‑1 (MQ=255)                                   tCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTGCCGTGCGGGAGCTGGAGTTGCAGGGGCGCAGTGAAGTGATTGGATTtgt                                                                                                  >  8:259220/1‑148 (MQ=255)                                     aaaGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAgg                                                                                                                                 >  8:49246/1‑114 (MQ=255)                                               tgAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCg                                                                                       >  2:186700/1‑147 (MQ=255)                                               tgAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGaa                                                                                     <  4:133417/149‑1 (MQ=255)                                               tgAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGaa                                                                                     <  2:285503/149‑1 (MQ=255)                                                 aaagaaaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAgt                                                                                   >  8:122054/1‑149 (MQ=255)                                                 aaagaaaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAgt                                                                                   >  1:191707/1‑149 (MQ=255)                                                 aaagaaaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAg                                                                                    <  5:194999/148‑1 (MQ=255)                                                   agaaaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAgtgt                                                                                 <  5:306530/149‑1 (MQ=255)                                                    gaaaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTc                                                                                <  7:122054/149‑1 (MQ=255)                                                      aaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCagtgaag                                                                                                               >  1:164153/1‑116 (MQ=255)                                                      aaTCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCagtgaag                                                                                                               <  2:164153/116‑1 (MQ=255)                                                                  cccGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAg                                                                                                                            >  2:148082/1‑91 (MQ=255)                                                                  cccGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAg                                                                                                                            <  1:148082/91‑1 (MQ=255)                                                                              gtaaaTCGACGTGCAGACGGCTTCGTGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTTGAGTTGc                                                                                                                              <  8:284707/75‑1 (MQ=255)                                                                              ggAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTTGAGTTGc                                                                                                                              >  7:289571/1‑77 (MQ=255)                                                                              ggAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTTGAGTTGc                                                                                                                              <  8:289571/77‑1 (MQ=255)                                                                              ggAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACAACCAGTTCCGTGCGGGAGCTTGAGTTGc                                                                                                                              >  7:284707/1‑77 (MQ=255)                                                                               gAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTc                                                    <  2:191707/149‑1 (MQ=255)                                                                               gAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTc                                                    >  2:5975/1‑149 (MQ=255)                                                                                     gACGTGCAGACGGCTTCGGGGAGCGTCACTTCGGTGATTACCACCCGTTCCGCGCGGCAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCTAGCTGTGATGCT‑GCGGGTCGatat                                               <  4:122991/148‑1 (MQ=255)                                                                                       cGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTgg                                            <  1:5975/149‑1 (MQ=255)                                                                                              aCGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTgggg                                          <  1:271103/144‑1 (MQ=255)                                                                                              aCGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTgggg                                          >  2:271103/1‑144 (MQ=255)                                                                                              aCGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGgctgct                                     >  2:345654/1‑149 (MQ=255)                                                                                                   ttCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTg                                  >  6:150381/1‑147 (MQ=255)                                                                                                   ttCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTg                                  >  3:338642/1‑147 (MQ=255)                                                                                                     cGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCGCGCGGGTCGATATTGGGGCTGCTTTGCa                                <  1:186700/148‑1 (MQ=255)                                                                                                      gggTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGcc                             >  7:105110/1‑149 (MQ=255)                                                                                                                  ttCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCtt                                    >  4:192275/1‑130 (MQ=255)                                                                                                                  ttCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCtt                                    <  3:192275/130‑1 (MQ=255)                                                                                                                    cGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGGGAAGTGATTGCCTTTGTGAAGTCGACCGAAGGGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACg                <  5:303063/148‑1 (MQ=255)                                                                                                                             ccaccCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCca      <  5:150381/149‑1 (MQ=255)                                                                                                                             ccaccCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCca      <  8:105110/149‑1 (MQ=255)                                                                                                                             ccaccCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCca      <  1:345654/149‑1 (MQ=255)                                                                                                                              caccCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCca      >  4:174095/1‑148 (MQ=255)                                                                                                                                 ccGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCCACAgg  >  4:42882/1‑149 (MQ=255)                                                                                                                                                   ggAGTTGCAGTTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGgctgct                                     <  8:302052/96‑1 (MQ=255)                                                                                                                                                   ggAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGgctgct                                     >  7:302052/1‑96 (MQ=255)                                                                                                                                                   |                                                                                                                                    ATGACTATCAAAGCCATCAACGTTCGCAACCAGTTCAAAGGTACTGTGAAAGAAATCCTCGAAGGCCCGGTACTGTCGGAAATCGACGTGCAGACGGCTTCGGGTATCGTCACTTCGGTGATTACCACCCGTTCCGTGCGGGAGCTGGAGTTGCAGGTGGGCAGTGAAGTGATTGCCTTTGTGAAGTCGACCGAAGTGTCTATCGCCAAGCTGTGATGCG‑GCGGGTCGATATTGGGGCTGCTTTGCAGCCCTTCGCCGGCACGCCCGCTCCCACAGG  >  NC_002947/295437‑295713

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