Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	385,475	A→C	8.4%	intergenic (+15/+20)	PP_0320 → / ← ltaE	methyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	385,475	0	A	C	8.4%	119.5 / 2.5	45	intergenic (+15/+20)	PP_0320/ltaE	methyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base A (26/15);  new base C (2/2);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 6.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

CAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCG  >  NC_002947/385338‑385623                                                                                                                                            |                                                                                                                                                     cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCg                                                                                                                                           >  7:147876/1‑149 (MQ=255) cAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCg                                                                                                                                           >  8:39454/1‑149 (MQ=255)  aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCgg                                                                                                                                          <  7:39454/149‑1 (MQ=255)  aGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGAAGGGGCCGCACGGCGGCCCCgg                                                                                                                                          >  3:6562/1‑149 (MQ=255)   ggggAGTTGCGCTTTGGCGCGGGAGTCGTCGGCGGAGGTGTAAGGGTGCCGCGAGGAGGGACAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGGGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCCGTGCGGCCCCTGc                                                                                                                                         <  8:265269/149‑1 (MQ=255)    gggATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGct                                                                                                                                        >  5:158197/1‑149 (MQ=255)    gggATTTCCGAGTTGGCGCGTGAGACGGCGGCGTAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TGCAGGCTCTAGATCGCAGGGGCCGCCGTGCGGCCCCTGc                                                                                                                                         <  6:253319/148‑1 (MQ=255)           ccGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCt                                                                                                                                  <  6:150910/148‑1 (MQ=255)            cGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCTGCTCTCGCTc                                                                                                                                 <  4:113969/148‑1 (MQ=255)                        tGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGAt                                                                                                                    >  1:291667/1‑149 (MQ=255)                           gaCGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCTGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGt                                                                                                                 >  5:144538/1‑149 (MQ=255)                                         tGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccgccg                                                                                                   <  6:158197/149‑1 (MQ=255)                                          gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccgccg                                                                                                   >  1:94036/1‑149 (MQ=255)                                          gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTGGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑GTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGGCACCAATGATCGTGCGGATATCcgccgc                                                                                                  <  5:293407/149‑1 (MQ=255)                                                gggCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTc                                                                                            <  5:388785/147‑1 (MQ=255)                                                 tGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCa                                                                                           <  8:284503/149‑1 (MQ=255)                                                         ggagTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGGGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCcg                                                                                   <  1:198792/149‑1 (MQ=255)                                                                 aGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGGCTTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGcttct                                                                            >  1:84025/1‑149 (MQ=255)                                                                        gCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATc                                                                    >  1:13873/1‑149 (MQ=255)                                                                                    tgatgaCCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACAt                                                        >  8:242122/1‑149 (MQ=255)                                                                                    tgatgaCCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACAt                                                        >  8:242123/1‑149 (MQ=255)                                                                                          ccT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcg                                                  <  2:13873/149‑1 (MQ=255)                                                                                           cT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAgcgc                                                 >  1:163529/1‑149 (MQ=255)                                                                                                      aGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGc                                       >  5:349957/1‑149 (MQ=255)                                                                                                         tgggtggg‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                         <  5:369332/112‑1 (MQ=255)                                                                                                         tgggtggg‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                         >  6:369332/1‑112 (MQ=255)                                                                                                          gggtggg‑tTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGTGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGcc                                   >  4:128908/1‑149 (MQ=255)                                                                                                           ggtggg‑tTCAGTCTATATATCGCAGCGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGcc                                   <  4:6562/148‑1 (MQ=255)                                                                                                                  ttCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAg                            >  7:123412/1‑149 (MQ=255)                                                                                                                    cAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCg                                        <  2:35673/135‑1 (MQ=255)                                                                                                                    cAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCg                                        >  1:35673/1‑135 (MQ=255)                                                                                                                         tctAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAg                      >  8:183722/1‑148 (MQ=255)                                                                                                                         tctAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAg                      >  7:65923/1‑148 (MQ=255)                                                                                                                           tAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAg                   >  3:98252/1‑149 (MQ=255)                                                                                                                             gATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGTGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGaa                    <  3:128908/146‑1 (MQ=255)                                                                                                                             gATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCAc                                          >  6:50536/1‑124 (MQ=255)                                                                                                                                agccggggCCGCCGTGCGGCCCCTGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAgcgc                >  1:236319/5‑147 (MQ=255)                                                                                                                                 gCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAg             >  5:192400/1‑149 (MQ=255)                                                                                                                                 gCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTACTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAg             >  5:187006/1‑149 (MQ=255)                                                                                                                                   cggggCCGCCGTGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGcc           >  6:363793/2‑149 (MQ=255)                                                                                                                                   aGGGGCCTCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGcc           <  8:58580/149‑1 (MQ=255)                                                                                                                                   aGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGcc           >  6:113641/1‑149 (MQ=255)                                                                                                                                       gtcGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTg       >  3:165349/3‑149 (MQ=255)                                                                                                                                         cGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCa              >  7:290079/1‑140 (MQ=255)                                                                                                                                         cGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCa              <  8:290079/140‑1 (MQ=255)                                                                                                                                            aCGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCg  >  5:129826/1‑149 (MQ=255)                                                                                                                                            |                                                                                                                                                     CAGGGGATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCT‑GGCGCGGCAGATGGGTGGG‑TTCAGGCTCTAGATCGCAGGGGCCGCACGGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCCCG  >  NC_002947/385338‑385623

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