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_0009134,223,6380TC96.5% 80.2 / ‑6.1 29intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (9/19);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 3.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

CCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CTTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTG  >  NC_000913/4223403‑4223878
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
ccTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGTGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTg                                                                                                                                                                                                                                     <  1:151489/251‑1 (MQ=255)
                  gctgTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCAttttt                                                                                                                                                                                                                   <  1:303855/251‑1 (MQ=255)
                         cATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGGCGCGGGTCGCGGGGGGGGGGGGG‑ACCATGGAAGTCT‑CCTTTTTTTTTTTC‑GTGGGAATCATTTTCATTTTTAtttttt                                                                                                                                                                                                            >  1:153605/1‑247 (MQ=255)
                                     cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGGTGCTCTAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACGCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTTTCCAGGCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATa                                                                                                                                                                                               <  1:34231/251‑1 (MQ=255)
                                         gggTGGGTCGTGGAATTGGGTAACCCGGCTTTTTGCGCCAGTTCCGTGAGTGCCACACTGCCAGTGGATTCGGCAATCCACTCCAGTAATTGCAGGCCACGCGTGAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGTGTGCG‑ACCATGACAGTCT‑CCTGTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTa                                                                                                                                                                                           <  2:316822/250‑1 (MQ=255)
                                                         tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTCGCGGGGATGGGTGCG‑GCCATGACAGGCT‑CCTTTTTTTTTATT‑GTGGGAATCAATTTTATTTTTTTTTTTAGCTTAAG‑CAAAAGTTACTGAACTGAttccatg                                                                                                                                                                           >  1:342540/1‑245 (MQ=255)
                                                         tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGATTTTTTCCCCCTTTTGCGGGGAGGGGGGGG‑GCCAAGGCAGTATACTTTTTT‑TGTTTC‑TTGGGGATCATTTTTATTTTTTTTGTTTTCTTATG‑CAATAGTTaatgaacttgatcggtg                                                                                                                                                                           >  2:315263/1‑242 (MQ=255)
                                                              aaCCCGGCTTGTTGCGCCAGTTCCGGGAGTGCCACACTGCCTTTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGtt                                                                                                                                                                       <  2:342540/250‑1 (MQ=255)
                                                                          tgCGCCAGTTCCGTGAGTTCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCGGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGTGTGCG‑ACCATGACAGTCT‑CCTTTTTCTTTATC‑GTGTAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACac                                                                                                                                                          <  2:126297/250‑2 (MQ=255)
                                                                           gcgcCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAAt                                                                                                                                                        <  1:294990/251‑1 (MQ=255)
                                                                              ccccttccgatcttGCCACACTGCCATTGGATTCGGCAACCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAAtctc                                                                                                                                                     <  2:129824/238‑1 (MQ=255)
                                                                                           tGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTTTGCCGCGTTTTGCGGGGATGGGGGCG‑GCCATGACAGTCT‑CCTTTTTCTGTTTC‑GTGGGAATCATTTTTATTTTTTTTTTTTGCTAATG‑CAATAATTACTGAAATTATTCGCTGAGTTAAT‑ATTGAACAAAactccgcacggaagagg                                                                                                                                        >  1:129824/1‑238 (MQ=255)
                                                                                                        ttGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTACAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGATTCGGGGGGAGGGGGGCGGACC‑TGACAGGTT‑CCTTTTTTTGTATC‑GTGGGAAACATTTTTATTTTTATTGTTAGTTAATGGCCA‑AGATAATTCACTGATTCGCTGAGGTAAT‑GTTGAACAAAACTAATGTTGCgtggggggcgcctttttc                                                                                                                            >  1:292730/1‑236 (MQ=255)
                                                                                                           gATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAg                                                                                                                        <  1:36992/251‑1 (MQ=255)
                                                                                                                ggCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGAGCGGTGGCAACGGCGGGCTTTATCCCCGGTTTCGCGCGGCTGGGGGCG‑TCCATGGCCGTCT‑CCTTTTTTTGTGTT‑CTGGGAATCATTTTTATTTTTTTTGTTTGCTAATG‑ACATAGTTACTGAACTTATACCATGGGGTAATAGTTGAA‑AAATCTCATGTTGTGGGGGGGTCGCATTTCCCCCAGATGCg                                                                                                                   >  1:250916/1‑251 (MQ=255)
                                                                                                                      ccACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGGTTTGCGGGGATGGGGGCG‑GGCATGAAAGTCT‑CCTTTTTTTGTATC‑GTGGAAATCCTTTTTCTTTTTTTTTTTAGCTAATT‑GAATAGTTACTGAAATGATCCGGTGAGGTTAT‑GTTGAACAAAACTCATGTTGCGTGGGGGTGGCGTTTTCcacaacagcggttttt                                                                                                             >  1:283395/1‑239 (MQ=255)
                                                                                                                          cccAGTAATTTCAGGCCACGCGTTAAAGACTGACCCTGTCCTGTCGCAGGGGCGGTGGCAACGGCTGGGTTTCTTCCGCGTTTTGCGGGAATGGGTGCG‑ACCATTACAGTCT‑CCGTTTTCGGTATC‑GTGGAAATCATTTTCATTTTTAGTGTTAGCTAATG‑CAATAGTGACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCTCTTTTACCACAGATTCGTTTATGCCAg                                                                                                         <  2:250916/250‑1 (MQ=255)
                                                                                                                                                 tAAAGACTGAACCTGTCCAGTCGCTGATGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTttatt                                                                                   <  2:283395/250‑1 (MQ=255)
                                                                                                                                                   aaGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTTCTGCCGCGGTTTGCGGGAATGGGGGGG‑ACCATTGCAGTCT‑CCTTTTTCTGTATCCGGGGAAA‑CATTTTTCTTTTTATTGTTATGTAATA‑CAATAGTTAATGAACTGATTCGCTGAGTTAAT‑GTTTAACAAATCTCATGGTGCGTGGGGGGCGCTTTTTACACAGATGCATTTATGTCAGTCTGGTTTGTTTCATTTTTTTTaa                                                                                 >  1:417997/1‑250 (MQ=255)
                                                                                                                                                    aGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTAGTGTTAGCTAATG‑CAATAGTTACTGAACTGAGCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTc                                                                                                                                     <  2:40305/197‑1 (MQ=255)
                                                                                                                                                               ggtcAGTCGCGGGCGCGGTGCCAACGTCGGGTGTTCTGCCTCGTTGCGCTGGAATGGGTGCG‑AACATGCTAATCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTGTATTCTTTGTGAATG‑CAATATCTTCTGAACTGATTCGATTAGTTAAT‑TTTGCACATATCTCAGGGTGCGTGGTTTTCGCTTTTCCCACAGATGTGTCGACGCCAGTAGGGTTTGTTTAATTTTTATTAAAttttgttggagc                                                                    <  2:877/248‑12 (MQ=255)
                                                                                                                                                                 ccAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCg                                                                   <  2:338833/250‑1 (MQ=255)
                                                                                                                                                                               ggTGGCAACGGCGGGTTTTCTTCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCgttatcaagtg                                                    <  2:394545/251‑11 (MQ=255)
                                                                                                                                                                                                                gggAATGGGTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGt                   >  1:44165/1‑251 (MQ=255)
                                                                                                                                                                                                                 ggAATGGGTGCG‑ACCTTGACATGCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTCTATTGTTAGCTAATG‑TAATAGTTACTGAGCTGAGCCGATGAGTGATT‑GTTGAACAAATCACATGTTGGGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAACGTTTTTTAAATCTGGGCTGAGCTTGTCGGGAGTAAGTGTGAGCAGCAAAGTGGAACTATTGCGTGCGCAGtt                  <  2:254903/251‑1 (MQ=255)
                                                                                                                                                                                                                      gggTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATCTTCATTTTTATTGTTAGCGAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGTTTTTACCACAGATGCGTTTATGCCAGTATGGTTTTTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGATCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATg             <  2:272973/251‑1 (MQ=255)
                                                                                                                                                                                                                       ggTGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGa            <  1:26125/251‑1 (MQ=255)
                                                                                                                                                                                                                         tGCG‑ACCATGACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATACGATGAGTTAAT‑GTTGAACAAATATCATGTTGCGTGGTGGTCGCTTTTACCACAGCTGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAAc          >  1:279854/1‑251 (MQ=255)
                                                                                                                                                                                                                                   gACAGTCT‑CCTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg  <  1:237748/250‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
CCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTCTGCCGCGTTTCGCGGGAATGGGTGCG‑ACCATGACAGTCT‑CTTTTTTCTGTATC‑GTGGAAATCATTTTCATTTTTATTGTTAGCTAATG‑CAATAGTTACTGAACTGATCCGATGAGTTAAT‑GTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTG  >  NC_000913/4223403‑4223878

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