Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	1,843,140	G→A	Q566* (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,140	0	G	A	97.9%	143.9 / NA	47	Q566* (CAG→TAG)	gacS	sensor protein GacS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/21);  minor base C (1/0);  total (26/21)
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

TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGG  >  NC_002947/1842993‑1843283                                                                                                                                                    |                                                                                                                                                tCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTaca                                                                                                                                               <  1:500069/150‑1 (MQ=255)       acaGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACt                                                                                                                                         <  1:1329701/150‑1 (MQ=255)              tACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTa                                                                                                                                  >  2:1400815/1‑150 (MQ=255)              tACGCATTGGCAGTGGTAACGCTCGAGCTGCAGATGATAATGGCCCAGTTGCACGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGCTCGTa                                                                                                                                  >  1:708050/1‑150 (MQ=255)               aCGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTag                                                                                                                                 <  1:1400815/150‑1 (MQ=255)                 gCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTagga                                                                                                                               <  1:1462553/150‑1 (MQ=255)                        cAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAAc                                                                                                                        <  1:1094611/150‑1 (MQ=255)                                   cTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCGCTGCTGGCGGCCTGCACCCCCTACAGCAAATGGTCGTAGGAGGCGAACAGGCTCACGCt                                                                                                             >  1:1128778/1‑150 (MQ=255)                                        gctgctGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGc                                                                                                        <  1:843048/150‑1 (MQ=255)                                            ctgATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCa                                                                                                    >  1:1028247/1‑150 (MQ=255)                                             tgatgaTAATGGCCCAGTTGCTTGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCAAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCGCCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAg                                                                                                   <  1:906750/150‑1 (MQ=255)                                                  atAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTc                                                                                              >  1:240636/1‑150 (MQ=255)                                                              ttGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTc                                                                                  >  1:1205319/1‑150 (MQ=255)                                                                  tCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGt                                                                              >  1:891965/1‑150 (MQ=255)                                                                    gggCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGaa                                                                                                                         <  2:80807/105‑1 (MQ=255)                                                                    gggCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGaa                                                                                                                         >  1:80807/1‑105 (MQ=255)                                                                             gTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTCCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGgcgcgc                                                                   >  2:354405/1‑150 (MQ=255)                                                                             gTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGgcgcgc                                                                   <  1:199721/150‑1 (MQ=255)                                                                                aCCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTg                                                                                        <  1:412576/126‑1 (MQ=255)                                                                                aCCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTg                                                                                        >  2:412576/1‑126 (MQ=255)                                                                                  ccAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGtt                                                              >  1:17680/1‑150 (MQ=255)                                                                                  ccAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCGGCACCCCCTACAGCAACTGGTCGAAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGtt                                                              >  2:1204037/1‑150 (MQ=255)                                                                                    aGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGc                                                                                                     >  1:1470158/1‑109 (MQ=255)                                                                                    aGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGc                                                                                                     <  2:1470158/109‑1 (MQ=255)                                                                                    aGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCg                                                            >  1:136926/1‑150 (MQ=255)                                                                                      gTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCaa                                                                                                                                           >  1:1334657/1‑69 (MQ=255)                                                                                      gTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCaa                                                                                                                                           <  2:1334657/69‑1 (MQ=255)                                                                                      gTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTg                                                          <  2:1269863/150‑1 (MQ=255)                                                                                           cACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTc                                                                                              >  2:280881/1‑109 (MQ=255)                                                                                           cACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTc                                                                                              <  1:280881/109‑1 (MQ=255)                                                                                           cACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCAt                                                     <  2:1205319/150‑1 (MQ=255)                                                                                            aCCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTGGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCATTGCCTGGCGCGCCACTTCGTGGCCATc                                                    >  2:929157/1‑150 (MQ=255)                                                                                                     cgcgAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATggc                                           >  2:908634/1‑150 (MQ=255)                                                                                                            tCGAAGGGCAGGCCGGCCTGGCTGGCGTCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGAAGAGCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcgg                                    >  1:90181/1‑150 (MQ=255)                                                                                                              gAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcggcg                                  >  1:84433/1‑150 (MQ=255)                                                                                                              gAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCAcggcg                                  >  1:1199812/1‑150 (MQ=255)                                                                                                                                gctggcGGCCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCgg                <  2:240636/150‑1 (MQ=255)                                                                                                                                      ggcCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAg          <  2:136926/150‑1 (MQ=255)                                                                                                                                       gcCTGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTAGCaaa         <  2:891965/150‑3 (MQ=255)                                                                                                                                          tGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAgc         >  2:1270449/1‑147 (MQ=255)                                                                                                                                          tGCACCCCCTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAgc         <  1:1270449/147‑1 (MQ=255)                                                                                                                                              cccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCg                                                                                                      >  2:1018915/1‑50 (MQ=255)                                                                                                                                              cccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCg                                                                                                      <  1:1018915/50‑1 (MQ=255)                                                                                                                                              cccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCgg  >  1:1226214/1‑150 (MQ=255)                                                                                                                                               ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCa                                                                                                                 <  2:289959/38‑1 (MQ=255)                                                                                                                                               ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCa                                                                                                                 >  1:289959/1‑38 (MQ=255)                                                                                                                                               ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCa                             >  2:489184/1‑122 (MQ=255)                                                                                                                                               ccccTACAGCAACTGGTCGTAGGAGGCGAACAGGCGCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCa                             <  1:489184/122‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TCGGGCACAGCACTACGCATTGGCAGTGGTAACGCTCGAGCTGCTGATGATAATGGCCCAGTTGCTCGGGCGACAGGTTACCCAGGTTGGCACCGAGCACCGCGAACTCGAAGGGCAGGCCGGCCTGGCTGGCGGCCTGCACCCCCTGCAGCAACTGGTCGTAGGAGGCGAACAGGCTCACGCTGAGGCCGCAGTCTTCCAGCTGGTGCTCCAGTGCCTGGCGCGCCAGTTCGTGGCCATCGACAATGGCGGCACGGCGGCCCAGCAATGGCTGCGGTGGCAGCGCTTCGG  >  NC_002947/1842993‑1843283

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