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,6380TC92.3% 77.6 / ‑5.1 26intergenic (‑10/‑190)iclR/metHDNA‑binding transcriptional repressor IclR/cobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (6/18);  total (8/18)
Fisher's exact test for biased strand distribution p-value = 8.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTG  >  NC_000913/4223398‑4223878
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
gCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTa                                                                                                                                                                                                                                         <  2:104373/251‑1 (MQ=255)
    aCCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑TTGGCGCGGTTCGCGGGGAAGGGGGCGACCAAGA‑CAGGCT‑CCATTTTTTGTAttgt                                                                                                                                                                                                                                     >  1:120777/1‑248 (MQ=255)
                           ttGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAGAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTAtt                                                                                                                                                                                                              <  2:162189/250‑1 (MQ=255)
                            tGCAGCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTg                                                                                                                                                                                                             <  2:128279/250‑1 (MQ=255)
                                          cAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGGTTTGCGGGGATGGGGGCGACCATGAACA‑TCT‑CCTTTTTTTGTTTTGTGGGAAACA‑TTTTTTTTTTTTTTTTTATCTAATGCAAt                                                                                                                                                                                               >  1:117103/1‑250 (MQ=255)
                                                         ggAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATc                                                                                                                                                                               <  1:191284/251‑1 (MQ=255)
                                                              tGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGAt                                                                                                                                                                           <  2:220376/250‑1 (MQ=255)
                                                                             gttgCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑TTGCCCCGGTTCGCGGGGATGGGGGGGGCCACGA‑CAGTCT‑CCTTTTTTTTTATCGTGGGAATCA‑TTTTCTTTTTTTTTGTTTGCTTATGCAATTGGTACTGAACTGATCCGGTGAGTTAATGTTgtac                                                                                                                                                            >  1:144501/1‑247 (MQ=255)
                                                                                        tCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGt                                                                                                                                                 <  2:7007/250‑1 (MQ=255)
                                                                                         ccGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGtt                                                                                                                                                <  1:153959/250‑1 (MQ=255)
                                                                                                  ccACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTAAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑TTGCCCCGTTTCGCGGGAATGGGGGGGACAATTA‑CACATT‑CCCTTTTTTGTATCGTTGGAATCA‑TTTTTATTTTTATTTTTAGCTAATACCATAGTTAATTGACTGATCCCATGAATTAATGTTGAACAAATATCCTGTTTCgttggagg                                                                                                                                      >  1:30103/1‑245 (MQ=255)
                                                                                                          ccATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTa                                                                                                                              <  1:84513/251‑1 (MQ=255)
                                                                                                             gtGGAGTCGGCAATCCACTCCAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTTAACTGATCCGATGAGTTAATGTTGAACAAAGCTCATGTGTCGTGGTGGTctcatctacca                                                                                                                           <  2:13635/250‑11 (MQ=255)
                                                                                                                        aaTCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTTTCTCCCCCGTT‑CGCGGGGATGGGGGGGGCCAAGA‑CCCTCTTCTTTTTT‑TTTTTCGTGGGAATCA‑TTTTTATTTTTTTTGTTTGGTAATGCCATCGTTACTGTAATTATCCGATGCGTTAATGTTGAACAAAACTTATTTTGCGTCGTGGTCGCtttttatacagacgtggtt                                                                                                                >  2:86112/1‑236 (MQ=255)
                                                                                                                                ccAGTAAGTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCGGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTGTTACCACAGATGCGTTTATGCCa                                                                                                          <  2:81499/249‑1 (MQ=255)
                                                                                                                                     aaTTTCAGGCCACGCGTTAAAGACTGAACCTGTCAAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGGCGCGGTTCGCGGGGATGGGGGGGGCCATCA‑AATTCTACTTTTTT‑TGTATTGTGGGAATCAATTTTC‑TTTTTTTTGTTAGCTCATGGAAAAGTTTCTTAAATTATTCCATGGGTTAATGTTGTACAAATCTTCTTTTTTGTGGTGGTCGCTTTTTCCCCAGATACGtttttggccgtttggg                                                                                                   >  1:221035/1‑238 (MQ=255)
                                                                                                                                        ttCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGttt                                                                                                 <  1:126035/250‑1 (MQ=255)
                                                                                                                                         tCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑TTGCCGCGTTTCGCGGGGATGGGGGGGACCATGA‑CAGGCT‑CCTTTTTTTGTATCGTGGAAAACA‑TTTTTATTTTTTTTGTTAACTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAAAAAAATTCATGTTGCGGGGGGGTCGCTTTTTCCACAGATGCGGTTATTCCAGGATGGTttgt                                                                                               >  1:76808/1‑251 (MQ=255)
                                                                                                                                           aGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTttgttg                                                                                             <  1:33309/251‑1 (MQ=255)
                                                                                                                                                 cgcgTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGGTTT‑TTGCCGCGGTTTGCGGGGATGGGTGGGGCCACGA‑CAGCCT‑CCCTTTTCTGTTTTGTGGGAATCA‑TTTTTATTTTTTTTTTTTTCTAATGCAATAGTTACTGAACTGCTCCGATGAGTTAATTTTGTACAAAACTCATGTTGCGGGGTGGTGGCTTTTTCCACAGATGCGTTTTTTACAGTATGGGTTGTTTAAtttt                                                                                       >  2:16668/1‑251 (MQ=255)
                                                                                                                                                                aCCTGTCCAGTCGGTGGTGCGGTGGCAACGGCGGGGTTT‑CTGCCGGGTTTCGCGGGAATGGGTTCGATCAGTA‑CATGCT‑CCGTTGTCTGTATCGTGGAAATCC‑TTTTCATTTTTAGTTTTAGCTAATGCAATAGTTACTGAACTGAACCGATGAGTTAATGTTGTACAAAGCTCATGGTGTGTTGTGGTCGCGTTTACCGCAGATGCGTTTATGCCAGTATGGTTTGTTGAAGTTTGATTAAATCTGGGtt                                                                        <  2:221035/251‑1 (MQ=255)
                                                                                                                                                                         gTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTc                                                               <  1:187422/251‑1 (MQ=255)
                                                                                                                                                                                                           ccGCGTTTCGCGGGAAGGGGTGCGACCATGA‑CAGTCT‑CCTCTTTCTGTATCGTTGAAATCA‑TTATCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGAGCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTGACCACAGGTGCGTTGATGCCAGTAGGGTTTGGTGAATATGTAGTAAATCTGGGTTGAGCGTGTCGGGAGCAACTGTGAGCAGCAAAGGGGgacaaca                              <  2:158189/251‑2 (MQ=255)
                                                                                                                                                                                                                                    ccATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt     <  1:211315/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccATGA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTAtt     <  2:48476/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gA‑CAGTCT‑CCTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTg  <  1:190703/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
GCCAACCTGACGCACGAAACCCTGCTGTTGCATCGTGGTTAGCAGGCGGTGGGTCGTGGAATTGGGTAACCCGGCTTGTTGCGCCAGTTCCGTGAGTGCCACACTGCCATTGGATTCGGCAATCCACTCCAGTAATTTCAGGCCACGCGTTAAAGACTGAACCTGTCCAGTCGCTGGTGCGGTGGCAACGGCGGGTTTT‑CTGCCGCGTTTCGCGGGAATGGGTGCGACCATGA‑CAGTCT‑CTTTTTTCTGTATCGTGGAAATCA‑TTTTCATTTTTATTGTTAGCTAATGCAATAGTTACTGAACTGATCCGATGAGTTAATGTTGAACAAATCTCATGTTGCGTGGTGGTCGCTTTTACCACAGATGCGTTTATGCCAGTATGGTTTGTTGAATTTTTATTAAATCTGGGTTGAGCGTGTCGGGAGCAAGTGTGAGCAGCAAAGTGGAACAACTGCGTGCGCAGTTAAATGAACGTATTCTG  >  NC_000913/4223398‑4223878

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