Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	244,829	A→G	16.2%	T169A (ACG→GCG)	algP →	transcriptional regulator AlgP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	244,829	0	A	G	16.2%	76.4 / 5.9	36	T169A (ACG→GCG)	algP	transcriptional regulator AlgP
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base G (3/3);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAAAGTACGCGAAGCTGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGC  >  NC_002947/244682‑244973                                                                                                                                                    |                                                                                                                                                 ggTAAAGTACGCGAAGCTGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAgacg                                                                                                                                                >  1:629052/1‑150 (MQ=255)   tAAAGTACGCGAAGCTGCAGGCAAGGCCCTTGATCAACGTGAGGCCAAACCGGCAGCCAAGGCCAACGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTCCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAgacgac                                                                                                                                              <  2:108939/150‑1 (MQ=255)         aCGCGAAGCTGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTgc                                                                                                                                        <  2:298875/150‑1 (MQ=255)           gcgAAGCTGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTgccg                                                                                                                                      <  1:586127/150‑1 (MQ=255)               aGCTGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCaa                                                                                                                                  <  1:490800/150‑1 (MQ=255)                  tGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCCGCTGCCAAGACGCCCACTGCCGCCAAgtg                                                                                                                               >  1:32127/1‑148 (MQ=255)                     aGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGc                                                                                                                            <  1:413681/150‑1 (MQ=255)                      ggCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCAATGCCGCCAAGCCTGCg                                                                                                                           >  1:97180/1‑150 (MQ=255)                               cTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCaa                                                                                                                                  >  2:405770/1‑134 (MQ=255)                               cTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCaa                                                                                                                                  <  1:405770/134‑1 (MQ=255)                                     cAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGATGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGcc                                                                                                            >  1:129747/1‑150 (MQ=255)                                              gCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCACGGCGACGGCTGCCGCCAAGCCTGCGGCCAAGCCGGCTGCGGAGGCGACCACTGCCGCAAAGGGTGCGGCCAAACAGGCTGCGCAGGCGAca                                                                                                   >  2:551719/1‑149 (MQ=14)                                                            cAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAg                                                                                                         >  2:621378/1‑130 (MQ=255)                                                            cAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAg                                                                                                         <  1:621378/130‑1 (MQ=255)                                                             aaGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAAc                                                                                                                    >  2:652936/1‑118 (MQ=255)                                                             aaGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAAc                                                                                                                    <  1:652936/118‑1 (MQ=255)                                                                    ccGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCa                                                                             >  2:674029/1‑150 (MQ=255)                                                                            gcCAAGCCTGCGGCCAAACCCCCTGCCAGGGCGACCGCCGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCt                                                                     <  2:129747/150‑1 (MQ=255)                                                                                aGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCa                                                                 >  1:646641/1‑150 (MQ=255)                                                                                 gCCTGCGGCCAAACCAGCAGCCAGGGCGCCCGCTGCCGCAAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCaa                                                                <  2:646641/150‑1 (MQ=255)                                                                                    tGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGc                                                             <  2:629052/150‑1 (MQ=255)                                                                                                     ccAGGGCTACCGGTGCCGCCCAGCTTGCGGCCAACCCAGCTGTCCTGGCGACCGCTGCCGCTAAGCCTGCGGCCCAACCGGCCGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGc                                            <  2:97180/150‑1 (MQ=11)                                                                                                        gggCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCt                                                                     >  1:514751/1‑122 (MQ=255)                                                                                                        gggCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCt                                                                     <  2:514751/122‑1 (MQ=255)                                                                                                                  tGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTCCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGTGGCCAAAcc                               <  1:342996/150‑1 (MQ=255)                                                                                                                         aaGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCg                                       <  1:389136/135‑1 (MQ=255)                                                                                                                         aaGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCg                                       >  2:389136/1‑135 (MQ=255)                                                                                                                                    cAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGc             >  1:213039/1‑150 (MQ=255)                                                                                                                                     aaaCCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTgc                                                                                              >  1:588049/1‑68 (MQ=255)                                                                                                                                     aaaCCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTgc                                                                                              <  2:588049/68‑1 (MQ=255)                                                                                                                                      aaCCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCACACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAAGCGGCTGCCAAGGCGAACGCTg           >  1:350665/1‑150 (MQ=255)                                                                                                                                             tGCCAAGGCGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCaa    <  1:673243/150‑1 (MQ=21)                                                                                                                                              gCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAg   >  2:673243/1‑150 (MQ=37)                                                                                                                                               ccAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGc                            <  2:187260/124‑1 (MQ=40)                                                                                                                                               ccAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGc                            >  1:187260/1‑124 (MQ=255)                                                                                                                                               ccAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGc  >  1:586283/1‑150 (MQ=255)                                                                                                                                                  aGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGc                                                       <  2:377372/94‑1 (MQ=255)                                                                                                                                                  aGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGc                                                       >  1:377372/1‑94 (MQ=255)                                                                                                                                                    |                                                                                                                                                 GGTAAAGTACGCGAAGCTGCAGGCAAGGCCCTTGATCAACGTGCGGCCAAACCGGCAGCCAAGGCAACCGCTGCTGCCAAGCCTGCGGCCAAACCAGCTGCCAGGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCAGCTGCCAAGACGACCACTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAACCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGCCTGCGGCCAAACCGGCTGCCAAGGCGACCGCTGCCGCCAAGC  >  NC_002947/244682‑244973

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑