Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,977,784	T→C	T21A (ACC→GCC)	flhC ←	DNA‑binding transcriptional dual regulator FlhC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,977,784	0	T	C	100.0%	80.9 / NA	25	T21A (ACC→GCC)	flhC	DNA‑binding transcriptional dual regulator FlhC
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/15);  total (10/15)

AGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAG  >  NC_000913/1977658‑1977910                                                                                                                               |                                                                                                                               aGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCaa                                                                                                                          >  1:57553/1‑133 (MQ=255)  gTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTTCCGCGCCGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAAt                                                                                                                         >  2:378217/1‑133 (MQ=255)   tCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAAtt                                                                                                                        <  2:11902/133‑1 (MQ=255)            aaTGGCAGCATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAAGGCGATCAATTTCATTGCCa                                                                                                               >  1:513504/1‑133 (MQ=255)                  agcaTGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGaa                                                                                                         <  2:57553/133‑1 (MQ=255)                                   tGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGa                                                                                        >  1:315500/1‑133 (MQ=255)                                     ggCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAAc                                                                                      <  2:61458/133‑1 (MQ=255)                                      gCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACa                                                                                     <  2:478514/133‑1 (MQ=255)                                                  ttCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATCTCCCGCGGTTCTTGAACAATCCTTTTTTCa                                                                         >  1:198444/1‑133 (MQ=255)                                                  ttCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCa                                                                         <  2:309152/133‑1 (MQ=255)                                                   tCTTTATAAATTTTTATCAGGCGTCCGCCACTGAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCAc                                                                        <  2:198444/133‑1 (MQ=255)                                                            aGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAAGATCCCGCGATTCGTGAGCAATGCTTGTTTCCCTCATGATCa                                                               >  2:129825/1‑133 (MQ=255)                                                              ttttATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAgg                                                             <  2:488672/133‑1 (MQ=255)                                                                        cGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCCCTCGTGATCAGGCGCTTTTCtt                                                   >  1:239498/1‑133 (MQ=255)                                                                            cgcgACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTgcgc                                               <  1:36631/133‑1 (MQ=255)                                                                                aCTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAgcg                                           <  2:82044/133‑1 (MQ=255)                                                                                aCTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAgcg                                           <  2:284190/133‑1 (MQ=255)                                                                                          gtTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAgg                                 <  1:13553/133‑1 (MQ=255)                                                                                          gtTTCGCTCTCCAGCGTCAGCATACGAGCGTCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAgg                                 <  1:129825/133‑1 (MQ=255)                                                                                           tttCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGc                                >  1:67720/1‑133 (MQ=255)                                                                                           tttCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGc                                >  2:509381/1‑133 (MQ=255)                                                                                              cGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGa                             >  1:463488/1‑133 (MQ=255)                                                                                               gCTTTCCAGCATCTGCAAACGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGAt                            <  1:509381/133‑1 (MQ=255)                                                                                                                aaCGAGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAg           >  1:348191/1‑133 (MQ=255)                                                                                                                    aGCGCCCAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCgtg       <  1:310702/133‑1 (MQ=255)                                                                                                                         ccAGGGCGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTgag  <  2:463488/133‑1 (MQ=255)                                                                                                                               |                                                                                                                               AGTCGGTTGAGAATGGCAGCATGCCTTTCGGCGGTGGGCTTCCGCGCAGTTCTTTATAAAGTTTTATCAGGCGTCCGCGACTTAACTGTGTTTCGCTTTCCAGCATCTGCAAACGAGCGCCCAGGGTGATCAATTCCATTGCCAGCTGAATATCCCGCGCTTCCTGAACAATGCTTTTTTCACTCATGATCAGGCCCTTTTCTTGCGCAGCGCTTCTTCAGGCTGATTAACATCATTCAGCAAGCGTGTTGAG  >  NC_000913/1977658‑1977910

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