Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,223,638	T→C	intergenic (‑10/‑190)	iclR ← / → metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,223,638	0	T	C	94.7%	52.8 / ‑6.1	19	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (6/12);  minor base G (0/1);  total (6/13)
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

CATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTC  >  NC_000913/4223428‑4223876                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTa                                                                                                                                                                                                         <  1:156844/251‑1 (MQ=255)    cGTGGTTATCAGGCTGTGGGTCGTGGAATTGGGGAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACTCTGCCATTTGATTCTGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGTGTGCGTCCATGACTGTCTCCTTTTGCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCt                                                                                                                                                                                                      <  1:218415/251‑1 (MQ=255)           aGCAGGCGGTGGGTCGTGGAATTGGGTATCCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGAGTCGGCAATCCACTGCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCTCGGGAATGGGTGCGTGCAGGACAGGCTTGTTTTTCTGTAGCGTGGAAATCATTTTCATTTGTATTGGTAGCTGATGCaa                                                                                                                                                                                               <  2:22424/251‑1 (MQ=255)                                        cccGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCCCGGTTTGCGGGGATGGGGGGGGCCATGGCAGTCTCCCTTTTTTTTTTCGTGGGAATCATTTTTATTTTTTTTTTTTGCTTATTCAATAATTACT‑TAACTGATCCGGTGAGtttat                                                                                                                                                                 >  1:194149/1‑248 (MQ=255)                                              ttgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGAAGGGGGCGACCATGACAGTCTCCTTTTTTTGTTTTGTGGGAATTATTTTTATTTTTTTTTTTTTCTAATTCAATAGTTACT‑GCACTGATCCGATGAGTTAAtttttt                                                                                                                                                            >  1:76941/1‑245 (MQ=255)                                               tggtgCGACAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGGCAATCCGCTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTGCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCGCCTGTTGCTGTATCGTGGAAATCATGTTCATTTTTATTTTTAGCTAATGCAATAGTTACT‑GAACTGATCCGGTTAGTTAATGTTGaa                                                                                                                                                           <  2:212455/247‑1 (MQ=255)                                                gttgCGCCAGTTCCGTGAGTGCCACACTTCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACa                                                                                                                                                         <  2:157251/251‑1 (MQ=255)                                                                 ggTGCCACACTGCCATGGGATGCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTGTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTGTATTGTGAGCTAATGCAATAGTTACT‑GAACTGATCCGATGACTTAATGTTGAACATATCTCATTTTGCgtgg                                                                                                                                        <  2:238434/250‑1 (MQ=255)                                                                         acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTCGGGGGGCTGGGGGGGGCCATGACAGACTCCCTTTTTTTTTTTGTGGACATCATTTTTATTTTTTTTTTTTACTTATTCAATATTTTCTTGAA‑TGCTCCCATGGGTTAATGTTGAACAAAtctcctttttcgtggggggtggc                                                                                                                                >  2:219313/1‑232 (MQ=255)                                                                              cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtcattttc                                                                                                                                                                                                                      <  2:233658/161‑1 (MQ=255)                                                                              cATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAc                                                                                                                           <  2:96383/251‑1 (MQ=255)                                                                                       cGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGATCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGc                                                                                                                  <  2:113550/251‑1 (MQ=255)                                                                                        ggCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCg                                                                                                                 <  1:242863/251‑1 (MQ=255)                                                                                                             cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgt                                                                                             <  2:194149/250‑1 (MQ=255)                                                                                                                                                      cGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCTACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAGCATTTTCATTTTTATTGTTAGCTAATGCAATATGTACT‑GAACTGGTCCGATGAGTTAATTGGGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcaca                                                                                                                       <  1:250435/183‑1 (MQ=255)                                                                                                                                                                                       cgGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTgcgc                    >  1:240902/1‑249 (MQ=255)                                                                                                                                                                                                    gACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGt      >  1:175348/1‑250 (MQ=255)                                                                                                                                                                                                     aCCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCGAATGCAATAGTTACT‑GAACGGGTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCGAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGATCGTGTCGGGAGCAAGTGTGAGCAGCAATGTGGGACAACTGCTTGCGCATTTAAATGAACCTAt    <  2:231293/251‑1 (MQ=255)                                                                                                                                                                                                       cATGAAAGTCTCCTTTTTCTGTATAGTGGAACTCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAACTCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTc  >  2:68609/1‑251 (MQ=255)                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                CATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACT‑GAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTC  >  NC_000913/4223428‑4223876

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