Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,994,755	Δ1 bp	16.7%	coding (1482/1599 nt)	PP_5251 →	transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,994,755	0	C	.	16.7%	99.4 / 16.5	36	coding (1482/1599 nt)	PP_5251	transporter
Reads supporting (aligned to +/- strand):  ref base C (10/20);  new base . (2/4);  total (12/24)
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.57e-01

CTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGT  >  NC_002947/5994608‑5994900                                                                                                                                                    |                                                                                                                                                   cTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCaccacc                                                                                                                                                   >  1:13568/1‑149 (MQ=255)     cagcagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGTTGGAGCATCTGCGACATCACCACCGGGt                                                                                                                                               <  1:203799/149‑1 (MQ=255)       gcagGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGa                                                                                                                                             >  6:216395/1‑149 (MQ=255)             ggcgCGGTAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gtg                                                                                                                                      <  1:122350/149‑1 (MQ=255)             ggcgCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCC‑ggtg                                                                                                                                     >  4:352170/1‑146 (MQ=255)              gcgcGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gtgg                                                                                                                                     <  7:178584/149‑1 (MQ=255)                gcgGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑gtggt                                                                                                                                    >  7:149340/1‑148 (MQ=255)                gcgGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGgtggtg                                                                                                                                   <  2:79453/149‑1 (MQ=255)                gcgGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGgtggtg                                                                                                                                   <  3:352170/149‑1 (MQ=255)                  gggAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGaa                                                                                                                                 <  4:291384/149‑1 (MQ=255)                   ggAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAAc                                                                                                                                <  1:284466/149‑1 (MQ=255)                    gAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCA‑CGGGTCACCCGGTGGTGatgt                                                                                                                               <  3:45992/149‑4 (MQ=255)                          gAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAg                                                                                                                           <  5:197158/147‑1 (MQ=255)                            aGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCg                                                                                                                       <  1:275373/149‑1 (MQ=255)                              cGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAg                                                                                                                     <  8:149340/149‑1 (MQ=255)                              cGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAg                                                                                                                     <  5:51666/149‑1 (MQ=255)                                   cgcaagcgcaaTGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGgcgc                                                                                                                 <  2:13568/148‑1 (MQ=255)                                              tGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgc                                                                                                     <  1:44753/149‑1 (MQ=255)                                              tGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGAtgc                                                                                                     <  4:252275/149‑1 (MQ=255)                                                                gTTGATGCGGGACAGTGCTGATGCGGCGTTGTGGCTGGTGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACGGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTGTGGATGCTGCTGACCGATGTGGTTg                                                                                   <  2:57281/149‑1 (MQ=255)                                                                     tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTg                                                                              >  7:76738/1‑149 (MQ=255)                                                                     tGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTg                                                                              >  6:253097/1‑149 (MQ=255)                                                                                gCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                 >  6:190162/1‑135 (MQ=255)                                                                                gCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAg                                                                                 <  5:190162/135‑1 (MQ=255)                                                                                   gATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATc                                                                 >  6:323905/1‑148 (MQ=255)                                                                                     tCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAAGCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTAt                                                              <  5:216395/149‑1 (MQ=255)                                                                                       cGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTg                                                            >  1:91269/1‑149 (MQ=255)                                                                                         gCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGc                                                          <  3:18091/149‑1 (MQ=255)                                                                                           gTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCtg                                                        >  5:240859/1‑149 (MQ=255)                                                                                                    ggAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGa                                                                                  <  7:2346/114‑1 (MQ=255)                                                                                                    ggAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGa                                                                                  >  8:2346/1‑114 (MQ=255)                                                                                                                cGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCacac                                    >  4:3196/1‑148 (MQ=255)                                                                                                                   gcaaTCGTGGCTGCACCCTCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTTGTGGGGAGCATAAGGCACACCAtt                                <  4:325088/147‑1 (MQ=255)                                                                                                                                gAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGccc                     <  2:91269/147‑1 (MQ=255)                                                                                                                                         cGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGc              <  5:253097/145‑1 (MQ=255)                                                                                                                                            gttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACgg               >  4:165978/3‑141 (MQ=255)                                                                                                                                            gttCACCA‑CGGGTCACCCGGTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACgg               <  3:165978/139‑1 (MQ=255)                                                                                                                                               caccacCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCgg    >  5:75304/1‑149 (MQ=255)                                                                                                                                                 ccaccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGt  >  6:104596/1‑149 (MQ=255)                                                                                                                                                  caccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGa                              <  2:135189/120‑1 (MQ=255)                                                                                                                                                  caccGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGa                              >  1:135189/1‑120 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTTCCAGCAGGCGGCGCGGGAGCTGGAAGCGATGCGCAAGCGCAATGGCGCATTGTCGTGGGGGTTGATGCGGGACAGTGCTGATCCGGCGTTGTGGCTGGAGTTCTTCTTCGAGGAATCGTGGCTGGAGCATCTGCGACATCACCACCGGGTGACCC‑GTGGTGAACTGAAGATCGAGGCGCGGGTCAGGATGCTGCAGACCGATGGGGTTGAGGTGCGGATTCGGCATCTATTGGCTGGTGGGGAGCATAAGGCACACCATTGATGCCGGCCCTTACGGCCCATTCGCGGGT  >  NC_002947/5994608‑5994900

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