Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,843,584	G→A	9.7%	Q404* (CAG→TAG)	gacS ←	sensor protein GacS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,843,584	0	G	A	9.7%	108.8 / 2.8	41	Q404* (CAG→TAG)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (19/18);  new base A (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

TCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAACGGTA  >  NC_002947/1843436‑1843724                                                                                                                                                     |                                                                                                                                              tCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTg                                                                                                                                               >  7:331711/1‑149 (MQ=255)   cTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCt                                                                                                                                             <  8:331711/149‑1 (MQ=255)    tGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGACCTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCtt                                                                                                                                            >  8:39581/1‑149 (MQ=255)                       cGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGcc                                                                                                                          >  4:235298/1‑148 (MQ=255)                              cTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTACTTGAGCCGCAACGGGTCGCCTACCAGc                                                                                                                   >  5:286207/1‑148 (MQ=255)                              cTGCACGCTGAAGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTACTTGAGCCGCAACGGGTCGCCTACCAGc                                                                                                                   >  1:128349/1‑148 (MQ=255)                                 cACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGaa                                                                                                                <  7:39581/148‑1 (MQ=255)                                        aTGCTCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                        <  3:235298/149‑1 (MQ=255)                                            gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAg                                                                                                            >  1:346702/1‑141 (MQ=255)                                            gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAg                                                                                                            <  2:346702/141‑1 (MQ=255)                                            gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                        <  5:81319/145‑1 (MQ=255)                                            gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGgtgt                                                                                                        >  6:81319/1‑145 (MQ=255)                                            gcAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTcg                                                                                                    >  7:76739/1‑149 (MQ=255)                                                               gTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGc                                                                                  >  6:342123/1‑148 (MQ=255)                                                                           cAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAgctgct                                                                     >  6:141249/1‑149 (MQ=255)                                                                            aGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCtt                                                                    >  1:276861/1‑149 (MQ=255)                                                                                  gcgcgcACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTg                                                              >  3:226787/1‑149 (MQ=255)                                                                                        acgacgATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGg‑gcggc                                                       >  8:50943/1‑149 (MQ=255)                                                                                            cgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTACTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGa                                                   <  2:128349/149‑1 (MQ=255)                                                                                            cgaTGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTACTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGa                                                   <  6:286207/149‑1 (MQ=255)                                                                                                ggTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                          <  6:102141/107‑1 (MQ=255)                                                                                                ggTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCa                                                                                          >  5:102141/1‑107 (MQ=255)                                                                                                   gCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCATCAATTCCAGCTGCTTg                                                                   >  5:349808/1‑127 (MQ=255)                                                                                                   gCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTg                                                                   <  6:349808/127‑1 (MQ=255)                                                                                                   gCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGa                                            >  1:206047/1‑149 (MQ=255)                                                                                                            ggTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGGTGCGGCCGGAGCGAGGATGGTCAgg                                    <  7:50943/149‑1 (MQ=255)                                                                                                              tGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGgtg                                 <  8:76739/149‑1 (MQ=255)                                                                                                               gAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGgtgt                                >  6:259439/1‑149 (MQ=255)                                                                                                                   ttGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGAAGGTCAGGGTGTCCTg                            >  5:346786/1‑149 (MQ=255)                                                                                                                      aTGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGAt                         >  3:289660/1‑149 (MQ=255)                                                                                                                       tGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTgg‑g                                                           <  2:51518/114‑1 (MQ=255)                                                                                                                         gCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCa                       <  2:276861/148‑1 (MQ=255)                                                                                                                           gTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGa                                               <  4:314039/122‑1 (MQ=255)                                                                                                                           gTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGa                                               >  3:314039/1‑122 (MQ=255)                                                                                                                              gCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTcc                  <  5:342123/148‑1 (MQ=255)                                                                                                                                  aCCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAgg                                    <  6:239309/126‑1 (MQ=255)                                                                                                                                   ccATGTTTGTGAGGATCTGCTTGAGCCGGAAGGGGTCGCCTACCAGCGAAGCAGGTGTGTAGCGGTAGATCAGGCTCAGCAATTCCAGCTGGTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAgg            <  4:226787/149‑1 (MQ=255)                                                                                                                                   ccAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAg                      <  4:19485/139‑1 (MQ=255)                                                                                                                                   ccAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAg                      >  3:19485/1‑139 (MQ=255)                                                                                                                                           gTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAACgg    >  5:195396/1‑149 (MQ=255)                                                                                                                                           gTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAACg     <  3:147116/148‑1 (MQ=255)                                                                                                                                              aGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAACGGTa  <  5:259439/148‑1 (MQ=255)                                                                                                                                                     gCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGa                                               <  6:225326/96‑1 (MQ=255)                                                                                                                                                     gCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGa                          >  6:4228/1‑117 (MQ=255)                                                                                                                                                     gCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGa                          <  5:4228/117‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              TCCTGAGGTGACAGGCCGATGCCGGTGTCCTGCACGCTGATGCGCAACTGCGCGCAGTCTTCGTGCTCGTCTTCCAGCATGGCGCGCACGACGATGGTGCCTTCGCGGGTGAACTTGATGGCGTTGCTGACCAGGTTGGTGAGGATCTGCTTGAGCCGCAACGGGTCGCCTACCAGCGAAGCAGGTGTGTCGCGGTAGATCAGGCTCAGCAATTCCAGCTGCTTGGCGTGG‑GCGGCCGGAGCGAGGATGGTCAGGGTGTCCTGGATCAGGTCCCGCAGGTTGAACGGTA  >  NC_002947/1843436‑1843724

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