Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,994,761	G→C	14.6%	V496V (GTG→GTC)	PP_5251 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,994,761	0	G	C	14.6%	135.9 / 14.2	54	V496V (GTG→GTC)	PP_5251	transporter
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base C (5/3);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCG  >  NC_002947/5994614‑5994907                                                                                                                                                    |                                                                                                                                                    gcagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGa                                                                                                                                                    >  6:72779/1‑149 (MQ=255)             ggAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAAc                                                                                                                                       <  4:305238/149‑1 (MQ=255)                 cTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTg                                                                                                                                     >  2:59167/1‑147 (MQ=255)                 cTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGgtggtg                                                                                                                                          >  3:163510/1‑142 (MQ=255)                 cTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGgtggtg                                                                                                                                          <  4:163510/142‑1 (MQ=255)                              gcaagcgcaaTGGCGCATTGTCGTGGTGGGTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGCGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCgg                                                                                                                      <  8:274133/149‑1 (MQ=255)                              gcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCgg                                                                                                                      <  8:239702/149‑1 (MQ=255)                                    gcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAgg                                                                                                                <  1:59167/149‑1 (MQ=255)                                      aaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAt                                                                                                              >  5:68391/1‑149 (MQ=255)                                       aTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGaa                                                                                                                                        <  8:167196/122‑1 (MQ=255)                                       aTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATACGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGaa                                                                                                                                        >  7:167196/1‑122 (MQ=255)                                            gcATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCa                                                                                                        >  4:39880/1‑149 (MQ=255)                                              aTTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGAt                                                                                                                                >  4:97473/1‑123 (MQ=255)                                              aTTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGAt                                                                                                                                <  3:97473/123‑1 (MQ=255)                                                    gTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATg                                                                                                >  6:42258/1‑149 (MQ=255)                                                     tGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAg                                                                                                       <  2:83324/141‑1 (MQ=255)                                                     tGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTAGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAg                                                                                                       >  1:83324/1‑141 (MQ=255)                                                       ggggTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATggg                                                                                              <  6:68391/148‑1 (MQ=255)                                                               tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                        >  3:113903/1‑146 (MQ=255)                                                               tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                        <  4:113903/146‑1 (MQ=255)                                                               tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTg                                                                                     >  7:120525/1‑149 (MQ=255)                                                                   ggACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGa                                                                                 >  1:60547/1‑149 (MQ=255)                                                                     aCAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGAtt                                                                               <  5:42258/149‑1 (MQ=255)                                                                             gaaccGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCAGCGACATCACCACCGTGTCCCCC‑CTGGTGAACTGAAGATCGAGCCGAGGGTCAGGATGCTGCAGCCCGAAGGCGTAGAGGTGCGGATTCGTCATc                                                                        >  4:137849/4‑148 (MQ=255)                                                                                 cGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTg                                                                   <  2:60547/149‑1 (MQ=255)                                                                                  ggCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTgg                                                                  >  4:145935/1‑149 (MQ=255)                                                                                          ggCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                                    >  3:31427/1‑123 (MQ=255)                                                                                          ggCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGc                                                                                    <  4:31427/123‑1 (MQ=255)                                                                                           gCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTggg                                                          <  8:120525/148‑1 (MQ=255)                                                                                            cTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTggg                                                          >  6:33981/1‑147 (MQ=255)                                                                                            cTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGa                                                        <  6:281619/149‑1 (MQ=255)                                                                                                       cttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                          <  1:298333/88‑1 (MQ=255)                                                                                                       cttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgctg                                                                                                          >  2:298333/1‑88 (MQ=255)                                                                                                       cttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                        >  1:269325/1‑106 (MQ=255)                                                                                                       cttcGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                        <  2:269325/106‑1 (MQ=255)                                                                                                           gAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCa                                         >  6:147765/1‑149 (MQ=255)                                                                                                                                  gCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCa                  >  1:93533/1‑149 (MQ=255)                                                                                                                                       aTCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgc             >  4:132477/1‑149 (MQ=255)                                                                                                                                        tCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgcg            >  1:117410/1‑149 (MQ=255)                                                                                                                                         cacca‑cGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCg                                <  2:52650/128‑1 (MQ=255)                                                                                                                                         cacca‑cGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCg                                >  1:52650/1‑128 (MQ=255)                                                                                                                                          accaccGGTTTACCC‑TTGTTTAACTGAAGATCGAGGCTCGTGTCATGATGCTGCATACCTATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGTGGAGAATAAGGCACATCATTGATGCCGGCCCTTACGGCCCATTCGCggg          <  7:47034/149‑1 (MQ=255)                                                                                                                                          accaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgcg            <  2:58635/147‑1 (MQ=255)                                                                                                                                          accaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTcgcg            >  1:58635/1‑147 (MQ=255)                                                                                                                                          accaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCgg           <  5:72779/148‑1 (MQ=255)                                                                                                                                          accaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCgg           <  3:137849/148‑1 (MQ=255)                                                                                                                                          accaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCgg           <  3:39880/148‑1 (MQ=255)                                                                                                                                          accaccGAGTGACCC‑GTGATAAACTAAAGATCGAGGCGCGGATCAGAATGCTGCAGACCGATGGGATTGAGGTGCGGATTCGGCATCTATTGGATGGTGGGGAGCATAAGGAACACAATTGATGCCGGCCATTACGGACCATTAGCggt          >  8:47034/1‑148 (MQ=255)                                                                                                                                           ccaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGt         >  2:138156/1‑149 (MQ=255)                                                                                                                                           accacGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTAtt                                                                    >  2:248210/5‑91 (MQ=255)                                                                                                                                           accacGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTAtt                                                                    <  1:248210/87‑1 (MQ=255)                                                                                                                                           accacGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCggg          >  6:32282/5‑149 (MQ=255)                                                                                                                                                 ggTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAAccc   <  1:138156/149‑1 (MQ=255)                                                                                                                                                  gTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCg  >  1:239048/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGTGAACCCG  >  NC_002947/5994614‑5994907

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