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	96.0%	77.3 / ‑5.8	25	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 (9/15);  total (10/15)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAC  >  NC_000913/4223394‑4223870                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         gTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTc                                                                                                                                                                                                                                    <  1:245495/251‑1 (MQ=255)        aaCCCGACGCGCGAAACCCTGCTGTTGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTGCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGGGCGGTGGCAACGGCGGGTTTGCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGTCAGTCTA‑CTTTTTCTGTATct                                                                                                                                                                                                                             <  2:26844/251‑2 (MQ=255)            ggACGCACGAAACCCTGCTGTTGCATCGGGGTTAGCAGGCGGTGGGTCGGGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCAGTCCAGTCGCTGGTGCGGTGTCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGGCTC‑CTTTTTCTGTATCGTGGa                                                                                                                                                                                                                         <  2:412631/250‑1 (MQ=255)                gCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAAtc                                                                                                                                                                                                                     <  1:2950/251‑1 (MQ=255)                        cccTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAAtcattttcat                                                                                                                                                                                                             <  1:327796/251‑1 (MQ=255)                                          gTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATGTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATg                                                                                                                                                                                           <  2:21047/251‑1 (MQ=255)                                             aGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTCGCGGGGATGGGGGCGGCCACGACAGTATC‑CCTTTTCTGTATCGTGGGAATCATTTTCATTTTTTTTGTTTGCTAATGCa                                                                                                                                                                                         >  1:296760/1‑250 (MQ=255)                                                gggCGGTGGGTTGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGGGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAg                                                                                                                                                                                     <  2:51120/250‑1 (MQ=255)                                                       gggTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGAGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGCCCCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTTCCGCGGTTCGCGGGAGTGGGTGCGACCACGACACTCTCCTTTTTTTTGT‑TCGGTGTAATCATTTTCATTTTTTTTGTTTGCTAAttcacaagatatta                                                                                                                                                                               >  1:290987/1‑237 (MQ=255)                                                                          cccGGCTTGGTGCGCCATTTCCGTGAGTGCCCCACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTTCTGCGGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTGAGCTAATGCAATAGTGACTTAACTGAGCCGATGAGTTAAt                                                                                                                                                           <  2:261970/251‑1 (MQ=255)                                                                                ttgttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGTTTCGCGGGGATGGGTGGGACCCTGACAGTCTC‑CTTTTTTTGTATTGGGGGAATACTTTTTATTTTTTTTTTTTGCTAATGCAATAATTACTGAACTTATCCGATGAGGTAATTTTgt                                                                                                                                                      >  1:340947/1‑249 (MQ=255)                                                                                                           acTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCt                                                                                                                          <  1:280599/251‑1 (MQ=255)                                                                                                                    ggATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACcaca                                                                                                                 <  2:369245/251‑1 (MQ=255)                                                                                                                      aTTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCCCGGTTCCCGGGGATGGGGGCGGCCCTGCAAGACTC‑CTTTTTTTTTATCGGGGGGATAATTTTCCTTTTTTTTGTTTGCTAATGCAATAAATACTGAAACGATCCGCTGTGGTTATTTTGAACAAAACTCATTTTTCGTGGTGGGGGCtttttcaccca                                                                                                                >  2:265061/1‑243 (MQ=255)                                                                                                                              aTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGTGCGGCCATGACAGTCTC‑CTTTTTTTTTATCGGGGGAATCATTTTTCTTTTTTTTTTTAGCTAATTCAATATTTACTGAACTGATTCGCTGAGTTAATGTTTTACAAATCTCATGTTTCGTGGTGGTCGCTTTTACCACAGATGTGtttt                                                                                                       >  2:138391/1‑250 (MQ=255)                                                                                                                               tCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTTGCCGCGGTTCGCGGGGATGGGGGCGGCCCTGACAGTCTC‑CTTTTTTTTTATCGTGGAAATCATTTTTCTTTTTTTTTTTTTCTAATGCAATAGTTACTGAAATGATCCGATGAGTTAATGTTGAACAAATCTCATGTTTCGTGGTGGTCGCTTTTACCACAAATGCGTTTAt                                                                                                      >  1:308201/1‑251 (MQ=255)                                                                                                                                 cACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGc                                                                                                    <  1:85466/251‑1 (MQ=255)                                                                                                                                        gTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATg                                                                                             <  1:341462/251‑1 (MQ=255)                                                                                                                                               cAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgtt                                                                                      <  1:164636/251‑1 (MQ=255)                                                                                                                                                                                                  ggcggGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGagcagc                                    >  1:71929/1‑250 (MQ=255)                                                                                                                                                                                                     aggTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGCTCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCaa                                  >  2:1762/2‑249 (MQ=255)                                                                                                                                                                                                                   gTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACtgcgt                  >  1:444897/1‑251 (MQ=255)                                                                                                                                                                                                                          gggAATGGGTGCGACCATGGTAGTCTC‑CTTTTTCTGTATCGTGGAAATCAGTTTCATTTGTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTAGCACAGATGCGTTTAGGGCAGGATGGTTTGTTGAATTGTGATTAAGGCTGGCTTGAGCGTGTCGAGAGCAAGGGTGAGCAGCAAAGTGGAACAACtgggggcgtagt           <  2:179368/251‑11 (MQ=255)                                                                                                                                                                                                                             aaTGGGTGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTaa        <  2:296760/251‑1 (MQ=255)                                                                                                                                                                                                                                   tGCGACCATGACAGTCTC‑CTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc  >  1:21151/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         GTTCGCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCGACCATGACAGTCTC‑TTTTTTCTGTATCGTGGAAATCATTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAC  >  NC_000913/4223394‑4223870

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