Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	385,476	C→G	6.8%	intergenic (+16/+19)	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,476	0	C	G	6.8%	167.3 / 2.2	54	intergenic (+16/+19)	PP_0320/ltaE	methyl‑accepting chemotaxis transducer/low specificity l‑threonine aldolase
Reads supporting (aligned to +/- strand):  ref base C (30/20);  new base G (2/2);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCC  >  NC_002947/385343‑385621                                                                                                                                       |                                                                                                                                                    gATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGTGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGTctct                                                                                                                                     <  7:329239/149‑1 (MQ=255) gATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGctct                                                                                                                                     >  5:76090/1‑149 (MQ=255)        gATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAg                                                                                                                              >  7:207519/1‑149 (MQ=255)             ggCGCGTGAGCCGGGGGCGGCGGTGCACGGGGCCCGCTCGGTGGGCCCGGGGTGGGGGGGGTGGGCTTATCACGTGGCGCGGCATTTGGGGGGGTGCC‑GGCTCTATAGCGCTGGGGCCGCAGG‑GCGGCCCCGGCTCTCGCTCAGccac                                                                                                                          <  4:159228/148‑1 (MQ=255)             ggCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGGGCGGCAGATGGTTGGGTTCA‑GGCTCTAGATCGCAGCGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGccac                                                                                                                          <  5:265094/148‑1 (MQ=255)                gcgTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAAt                                                                                                                      >  8:49779/1‑149 (MQ=255)                      gaCGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGt                                                                                                                <  1:329588/149‑1 (MQ=255)                               ggAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATc                                                                                                       >  4:212363/1‑149 (MQ=255)                                 agatGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATccg                                                                                                     >  5:97141/4‑149 (MQ=255)                                     gtaaGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                 >  8:178484/3‑149 (MQ=255)                                     gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCCGT‑GCGGccc                                                                                                                                            >  4:46636/1‑106 (MQ=255)                                     gCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCCGT‑GCGGccc                                                                                                                                            <  3:46636/106‑1 (MQ=255)                                      cAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCcgccgc                                                                                                 >  5:160467/1‑148 (MQ=255)                                            tGCCGCGAGGATGGCCCGGCGTGGCGTGGGTTGGCGGCTGACCGGGCGCGGCAGATGGGGGGGGGCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCa                                                                                          <  4:292926/149‑1 (MQ=255)                                              ccGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACg                                                                                        >  6:10381/1‑149 (MQ=255)                                               cgcgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGc                                                                                       >  4:247715/1‑149 (MQ=255)                                                gcgAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCa                                                                                      >  5:321200/1‑149 (MQ=255)                                                    ggagTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCcg                                                                                  >  2:71943/1‑149 (MQ=255)                                                       gTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCgcgc                                                                               >  3:225926/1‑149 (MQ=255)                                                        tGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCgcgc                                                                               <  3:212363/148‑1 (MQ=255)                                                          ccAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGAt                                                                                                                   <  1:105303/110‑1 (MQ=255)                                                          ccAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGAt                                                                                                                   >  2:105303/1‑110 (MQ=255)                                                             ggCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCttcttc                                                                         <  5:10381/149‑1 (MQ=255)                                                             ggCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCttcttc                                                                         >  4:120143/1‑149 (MQ=255)                                                              gCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                        <  3:120143/149‑1 (MQ=255)                                                              gCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCg                                                                        <  3:247715/149‑1 (MQ=255)                                                                   gCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGccac                                                                                                                          >  5:5376/1‑94 (MQ=255)                                                                          ttGGCTTATGACCTGGCGGGGCAGATGGGGGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAg                                                           <  6:76090/149‑1 (MQ=255)                                                                                gatgaCCTGGCGCGGCAGATGGGTGGGTTCAGGGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACAt                                                      >  1:44585/1‑149 (MQ=255)                                                                                    aCCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAgc                                                 >  4:39720/1‑149 (MQ=255)                                                                                     ccTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAgc                                                 <  7:178484/148‑1 (MQ=255)                                                                                      cTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAgcgc                                               >  6:88830/1‑149 (MQ=255)                                                                                        ggCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAgcgcgc                                             >  4:351400/1‑149 (MQ=255)                                                                                           gcgGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACGAGCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAgcgcgcgc                                           >  1:76645/1‑149 (MQ=255)                                                                                                  aTGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGaa                                                   <  4:310942/133‑1 (MQ=255)                                                                                                  aTGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGaa                                                   >  3:310942/1‑133 (MQ=255)                                                                                                   tgggtgggGTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATCTCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGc                                  <  8:207519/149‑1 (MQ=255)                                                                                                    gggtgggtTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGc                                  <  6:160467/148‑1 (MQ=255)                                                                                                      gtgggtTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGa                               >  4:344271/1‑149 (MQ=255)                                                                                                          gtTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATaaa                           <  6:97141/149‑1 (MQ=255)                                                                                                           ttCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAg                          <  4:225926/149‑1 (MQ=255)                                                                                                           ttCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAg                          >  3:171239/1‑149 (MQ=255)                                                                                                                ggCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGta                       >  1:240511/1‑148 (MQ=255)                                                                                                                 gCTCTAGATCGCCGGGGCCGCACG‑GCGGCCCCGGCGGTCGCTGAGCCCCCAATGCTCGTGCGGGTATGCGCCGCCCGCTCACGCCCCCGCGCGTGTTGGGTGT‑CCCACGCGCACCTGAAGCGCGCGCCCCCGCGGCCGCTAAAGGta                       >  2:345585/1‑147 (MQ=255)                                                                                                                 gCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTATCCCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGtat                      <  1:71943/149‑1 (MQ=255)                                                                                                                  ctctAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGta                       >  4:301143/1‑146 (MQ=255)                                                                                                                   tctAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCCCGCACCCGCGCTTCTTCGGTAT‑CCCCCGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGGATAg                    >  7:234715/1‑148 (MQ=255)                                                                                                                    ctAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAg                    >  5:379133/1‑147 (MQ=255)                                                                                                                         tCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCTCCCTCGCGTCTTCGGTGT‑CCCCCGAGCACATGAAGCGCGCGCCCCCGCTGCCGAGATTGGTATAGactcgg              >  8:319583/1‑143 (MQ=255)                                                                                                                           gccggggCCGCCGT‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAg           >  6:88964/4‑149 (MQ=255)                                                                                                                                ggCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCtt      <  2:44585/149‑1 (MQ=255)                                                                                                                                ggCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCtt      <  6:321200/149‑1 (MQ=255)                                                                                                                                  ccGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGc    >  3:134534/1‑149 (MQ=255)                                                                                                                                    gCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGccc  <  6:379133/149‑1 (MQ=255)                                                                                                                                       |                                                                                                                                                    GATTTCCGATTTGGCGCGTGAGACGGCGGCGGAGGTGCAAGGGTGCCGCGAGGAGTGCCAGGCGTTGCGTGGGTTGGCTGATGACCTGGCGCGGCAGATGGGTGGGTTCA‑GGCTCTAGATCGCAGGGGCCGCACG‑GCGGCCCCGGCTCTCGCTCAGCCACCAATGATCGTGCGGATATCCGCCGCCAGCTCACGCACCCGCGCTTCTTCGGTAT‑CCCACGAGCACATGAAGCGCGCGCCACCGCTGCCGATAAAGGTATAGAAGCGCCAGCCCTTGCCC  >  NC_002947/385343‑385621

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