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	87.0%	44.0 / ‑2.0	23	intergenic (‑10/‑190)	iclR/metH	DNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (4/16);  total (5/18)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.22e-01

CCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223399‑4223880                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    ccatccTGACGCCCGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGATTGCCACACTGCCATTTGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAATCTGTCCAGTCTCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCTTTTCGCGGGAATGGGTGAGACCATGACAGTCTCTTTTTTCTGt                                                                                                                                                                                                                                            <  1:247186/245‑1 (MQ=255)                   cccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAAtca‑ttttcat                                                                                                                                                                                                                       <  1:84229/251‑1 (MQ=255)                    ccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTCTCTGTATCGTGGAAAtca‑ttttcatt                                                                                                                                                                                                                      <  1:85387/251‑1 (MQ=255)                       gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATca‑ttttcatttt                                                                                                                                                                                                                    <  1:3616/250‑1 (MQ=255)                          gTTGCATCGTGGTTAGCATGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTTCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCTCTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCATTCGCTGGTGCGTTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATTACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTAt                                                                                                                                                                                                                 <  1:331083/250‑1 (MQ=255)                            tGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAAGCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTg                                                                                                                                                                                                               <  1:384299/250‑1 (MQ=255)                                 cGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTAGCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCt                                                                                                                                                                                                         <  1:267675/251‑1 (MQ=255)                                        aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTTTTAGCTAATGCaa                                                                                                                                                                                                  <  1:72762/251‑1 (MQ=255)                                                                ggTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAg                                                                                                                                                                          <  1:260093/251‑1 (MQ=255)                                                                                                       cTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTTCCCCGTTTCGCGGGAATGGGGGCGGCCCTGACAGTCTCTTTTTTTTTTTTCGTGGGAATCAATTTTC‑TTTTTTTTTTTTTCTTATGCAATAATTACTTTAATGATCCCATAGGTTTATGTTTAACAAATTTCATATTTCGTGGGGGGCGCtt                                                                                                                                   >  1:145468/1‑251 (MQ=255)                                                                                                         gCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGGCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCtt                                                                                                                                   <  1:109381/249‑1 (MQ=255)                                                                                                               ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTTGGGGGGATGGGTGGGGCCAAGACCGTCTCCTTTTTTTGTATAGTGGGAATCC‑TTTTTATTTTTTTTTTTTACTAAAGCAATAGTTTCTTAAATTATCCGATGAGTTTATGTTTTACAAAATTCATGTTTCGtggtgggtggtttttcccc                                                                                                                            >  1:340900/1‑237 (MQ=255)                                                                                                                  ttCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGGTTTGCGGGGATGGGTGGGGCCCTGACAGGCTCCTTTTTTTGTTTCGTGGGAATCA‑TTTTTATTTTTTTTGTTTACTAATTAAAAAGGTACTGAACTGATACGATGAGTTAATTTTTGACAAATTTCCTTTTGCGGGGGGGGCGATTTTACcac                                                                                                                            >  1:261159/1‑247 (MQ=255)                                                                                                                            cACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCAt                                                                                                                                                     <  1:176257/212‑1 (MQ=255)                                                                                                                                         tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttg                                                                                                  <  1:375208/250‑1 (MQ=255)                                                                                                                                           aGGCCACGCGTTAAAGGCGGAACCTGTACAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGATAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCTTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgttg                                                                                               <  2:261159/251‑1 (MQ=255)                                                                                                                                                 cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGTTTCGCGGGGATGGGGGCGACCCTGACAGTCTCCTTTTTCTTTATCGTGGGAAACC‑TTTTCATTTTTTTTGTTAGCTTATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAAAAAATCTCATGTTTCGTGGTGGTCGCTTTTACCACAGAATCGGTTATGGCAGGATGGGTTGTTTAAtttt                                                                                         >  1:126766/1‑251 (MQ=255)                                                                                                                                                                                                 ggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGt                                         <  1:41623/251‑1 (MQ=255)                                                                                                                                                                                                      tCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGaa                                     <  1:97931/250‑1 (MQ=255)                                                                                                                                                                                                                     gggAATTGGTGCGACCATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTCCTGATCTGTTCCGATGTGTTAATGTTGAACAAATCTCATGTTGCGTGTTGGTCGCTGTTACCACAGATGCGTTTATGCCAGTAGGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTTTCGGGAGCCAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                     <  2:340900/251‑1 (MQ=255)                                                                                                                                                                                                                                   ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgtt                                                                                                <  1:72317/162‑1 (MQ=255)                                                                                                                                                                                                                                   ccATGACAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgtt                                                                                                >  2:72317/1‑162 (MQ=255)                                                                                                                                                                                                                                        aCAGTCTCCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAGCAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGt  <  1:311940/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTCTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTGGT  >  NC_000913/4223399‑4223880

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