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	95.7%	67.9 / ‑4.9	23	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (5/17);  total (6/17)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCTTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCT  >  NC_000913/4223397‑4223883                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       cGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGt                                                                                                                                                                                                                                               <  1:55387/251‑1 (MQ=255)                                                   tGGGTCGTGGAATTGGGTAACCCGGCTTGCTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTTCCGCGTTTCGCGGGAATGGGGGCGACCCTGACAGT‑CTCCTTTTTTTGTATCG‑TGGGAATCATTTTTATTTTTTTTTTTAGCTTAAGCAATaatttttt                                                                                                                                                                                            >  1:58058/1‑244 (MQ=255)                                                         gTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGTGTTTTTTGCCGCGTTTCCGGGGAATGGGGGGGGCCAGGGCGGG‑CTGCTTTTTCTGCATCGGTGGAA‑TCATTTTTTTTTTTATTGTTTGCCAATTCAAAAGTTaataaaatgt                                                                                                                                                                                      >  1:164046/1‑242 (MQ=255)                                                            gAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATcc                                                                                                                                                                                   <  2:277095/251‑1 (MQ=255)                                                                 gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCAGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATg                                                                                                                                                                               <  2:186575/250‑1 (MQ=255)                                                                 gggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGGGGCCATGACAGT‑CTCCTTTTTTTGTATCG‑TGGAAATCATTTTTATTTTTTTTTTTTGCTTATTCAATAGGTACTGAACTGATCCGATg                                                                                                                                                                               >  1:23386/1‑250 (MQ=255)                                                                   gTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAACCGATGAGt                                                                                                                                                                            <  2:26242/251‑1 (MQ=255)                                                                     aaCCCGGCTTGCTGCGCCAGTTCCGTGAGTGCCACACTGCCATTTTAGTCGGCAATCCACTCCAGTTATTTCAGGTCACGCGTTAAAGACTTAACCTGTCCAGTCGCTGGTGCGGTGGCAACGTCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATGTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCTATGAGTTa                                                                                                                                                                          <  2:264884/251‑1 (MQ=255)                                                                                     ccAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAAtctc                                                                                                                                                          <  1:24650/251‑1 (MQ=255)                                                                                      cAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCTGATGAGTTAATGTTGAACAGATCTCa                                                                                                                                                         <  2:195476/251‑1 (MQ=255)                                                                                       aGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAt                                                                                                                                                        <  1:67139/251‑1 (MQ=255)                                                                                            cGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGc                                                                                                                                                   <  2:47030/251‑1 (MQ=255)                                                                                             gTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCCCGGGGGTGGGGGGGGCAACGGCAGGGCTCTTTTTT‑TTTTTCG‑TGGGAATCATTTTTTTTTTTATTGTTAGCTAATTCAAAAATTATATTACTGGTTCGCTAGGTTAATGTTGAACAAATCTCATTTTGCg                                                                                                                                                  >  2:308632/1‑251 (MQ=255)                                                                                               gAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACGGATCCGATGAGTTAAGGTTGAACAAATCTCATGTTGCgtg                                                                                                                                                <  2:184155/251‑1 (MQ=255)                                                                                                   gCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGGTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTTCGtggtgg                                                                                                                                            <  1:59990/251‑1 (MQ=255)                                                                                                                                                         aaaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatta                                                                                       <  1:300675/250‑1 (MQ=255)                                                                                                                                                                                         ggCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGa                                                      <  1:158054/251‑1 (MQ=255)                                                                                                                                                                                            aaCGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagca                                                   <  1:227457/251‑1 (MQ=255)                                                                                                                                                                                                     ttttCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTgg                                          >  2:315404/1‑251 (MQ=255)                                                                                                                                                                                                                                     ccATGACAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt          <  2:246301/251‑1 (MQ=255)                                                                                                                                                                                                                                           cAGT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTg    <  1:66091/251‑1 (MQ=255)                                                                                                                                                                                                                                             gT‑CTCCTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCt  <  2:23386/251‑1 (MQ=255)                                                                                                                                                                                                                                             gT‑CTCCTTTTTATGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTACATGAACGTATTCTGGTGCt  >  2:246076/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                       CGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGT‑CTCTTTTTTCTGTATCG‑TGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGTGCT  >  NC_000913/4223397‑4223883

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