Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,095,189	G→C	11.5%	intergenic (‑240/+77)	htpX ← / ← alaA	protease HtpX/glutamate‑pyruvate aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,095,189	0	G	C	11.5%	119.1 / 3.5	42	intergenic (‑240/+77)	htpX/alaA	protease HtpX/glutamate‑pyruvate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (14/23);  new base C (2/3);  total (16/26)
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

ACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCAC  >  NC_002947/2095048‑2095328                                                                                                                                              |                                                                                                                                            aCGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCgcgggcgg                                                                                                                                      >  1:295399/1‑149 (MQ=255)   ggcggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGAGGgcggcg                                                                                                                                    >  5:59618/1‑149 (MQ=255)    gcggAGGATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTGGCAGGATGTAAGAAACATTCGGCTGGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGGGGCGCAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                                   <  1:69640/149‑1 (MQ=37)     cggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                  <  2:189591/149‑1 (MQ=255)     cggAGTAGGGCAGAGGCTGGCGGGGGGGCGGTGGGATAGAGCGGTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCt                                                                                                                                  <  3:187989/149‑1 (MQ=255)      ggAGTATGGCAGCGGAGGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGGGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTc                                                                                                                                 <  1:73115/149‑1 (MQ=39)      ggAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTc                                                                                                                                 <  6:273575/149‑1 (MQ=255)      ggAGTATGGAAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGTAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCTCCGCGCGGGCGGCGCTc                                                                                                                                 <  4:264997/149‑1 (MQ=255)        aGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGa                                                                                                                               >  1:91148/1‑149 (MQ=255)            tGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                   >  5:326270/1‑141 (MQ=255)            tGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGgcgc                                                                                                                                   <  6:326270/141‑1 (MQ=255)                 gagGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGccc                                                                                                                      <  2:118919/149‑1 (MQ=255)                  agGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGcccc                                                                                                                     >  2:122619/1‑149 (MQ=255)                   ggATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCcca                                                                                                                    >  4:386372/1‑149 (MQ=255)                        gCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCccaccac                                                                                                                <  2:91148/148‑1 (MQ=255)                        gCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                               <  7:289339/133‑1 (MQ=255)                        gCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTCGa                                                                                                                               >  8:289339/1‑133 (MQ=255)                        gCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGCCTCCGGCAAGGCCGCTGGTTGTTTGCCCAAATTGTTTTGTTGCATTTGATATCAACCGCCGCCCGGCCGGCCCGCGATATGCACCcccccacc                                                                                                               >  1:137133/1‑149 (MQ=255)                                        aTAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCCCCCCCCCCCCCGCGGCGCAc                                                                                                >  4:18355/1‑148 (MQ=39)                                                    aGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGcccc                                                                                   <  6:175437/149‑1 (MQ=255)                                                       aTGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCCCCCCCCCGGGGAACCCCGCGGAGccccccg                                                                                >  6:253910/1‑147 (MQ=38)                                                         gTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCCCCACCACCCCCGCGGGGGACCTCTCGTAGCCCCCcgcg                                                                              >  6:291771/1‑149 (MQ=39)                                                            aGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCCCCACCCCCCCCCCGGGGAACCGCGAGGAGCCCCCGGCGccc                                                                           >  4:89422/1‑149 (MQ=38)                                                                  aTTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCCCCCCCCCGGGGCACCGGTCGGAGCCCCCCGCGccccccgcc                                                                     >  6:369785/1‑145 (MQ=37)                                                                  aTTCGGCTTGGGGGCGGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGGGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCg                                                                     <  5:253910/149‑1 (MQ=255)                                                                       gCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCCCCCCCCCCCCCCCGTGGGACCGCGCGTGGcccccggcgcccccggcggcccc                                                                >  4:128994/1‑144 (MQ=37)                                                                            ggggCCGCGGCAAGGCCGCTGGGTGTTTGCCGCAAGTGGGGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTg                                                           <  4:41595/149‑1 (MQ=255)                                                                               gCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGa                                                        <  3:200723/149‑1 (MQ=255)                                                                                    gTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAAc                                                   <  2:137133/149‑1 (MQ=255)                                                                                         ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCACCCCCGCGGGGAACCTCTCGGAGCCCCCGGCGCCCCCGGCGAACCCAACTGCGAAAAAcccat                                               >  6:168057/1‑146 (MQ=38)                                                                                         ggCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCCCCACCCCCCCGGCGACCCGCTCGTAGCCCCCGGGGCCCCCGGCGCACCCGACGGCGAATAAcccgt                                               >  4:55240/1‑146 (MQ=37)                                                                                            cGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAgg                                           <  7:6629/149‑1 (MQ=255)                                                                                               tGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCACCGTCGAACCCTACTGCGAATAAccc                                                 <  1:12245/140‑1 (MQ=255)                                                                                                 ggTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAAtt                                      <  3:18355/149‑1 (MQ=255)                                                                                                   tgtTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGc                                    <  5:400197/149‑1 (MQ=255)                                                                                                    gtTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCAACCCCACCCCCCCGGCGGCCCGGGGCTTGCCCCCCCCCCCCCCCGAGCACCCAACCGCGAATAACCCCTCACGAAAATggc                                   >  2:9505/1‑147 (MQ=12)                                                                                                     tttGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGcc                                   <  8:251462/148‑1 (MQ=255)                                                                                                              aGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCa                         <  5:71395/149‑1 (MQ=255)                                                                                                                             gtgAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAgg          <  1:184130/149‑1 (MQ=255)                                                                                                                                   tCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATcc    <  8:4001/149‑1 (MQ=255)                                                                                                                                   tCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATcc    <  1:122619/149‑1 (MQ=255)                                                                                                                                     gAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTa                                                              <  5:352763/89‑1 (MQ=255)                                                                                                                                     gAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCAc  <  2:295399/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                            ACGGCGGAGTATGGCAGAGGATGGCGTGGGGGCTGTGGGATAGAGCGTTGCAGGATGTAAGAAACATTCGGCTTGGGGGCTGCGTCAAGGCCGCTGGGTGTTTGCCGCAAGTGGTGTGGTGCCTGTGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTGCGCCCCACCACCACCCCCGCGTCGAACCTCTCGTAGCCCCCGGCGCCCCCGTCGAACCCTACTGCGAATAACCCTTCAGGAAATTGCCGATCCGGCCAATTGCCGCTTCCAGGTCATCCAC  >  NC_002947/2095048‑2095328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑