Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,099 C→T intergenic (‑130/‑223) yifB ← / → ilvL putative magnesium chelatase YifB/ilvXGMEDA operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,0990CT100.0% 58.9 / NA 19intergenic (‑130/‑223)yifB/ilvLputative magnesium chelatase YifB/ilvXGMEDA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/11);  total (8/11)

TCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAG  >  NC_000913/3949854‑3950340
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
tCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGc                                                                                                                                                                                                                                               <  2:37002/250‑1 (MQ=255)
            accaTCGTTAAGCCGGGAAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCgg                                                                                                                                                                                                                                  >  2:407838/1‑251 (MQ=255)
                        ccGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGaa                                                                                                                                                                                                                      <  1:195500/251‑1 (MQ=255)
                        ccGGGGAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGa                                                                                                                                                                                                                       >  1:349696/1‑250 (MQ=255)
                                      aCTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGa                                                                                                                                                                                                        <  2:61290/251‑1 (MQ=255)
                                                           gATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGATTTGTTCTTGCGGTGGAATTATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTg                                                                                                                                                                                   <  2:248750/251‑1 (MQ=255)
                                                                                 aGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGGACAACATTCATACTGAAATTGAATTTTTTTCACTCACtatttt                                                                                                                                                             >  2:433145/1‑251 (MQ=255)
                                                                                 aGGGCTGCGCGGGTATCAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACtatttt                                                                                                                                                             <  1:45838/251‑1 (MQ=255)
                                                                                                            aGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                                                                                                                                  <  1:433145/251‑1 (MQ=255)
                                                                                                                     aGCCCTCCTTGAGTCACCAGTATGTGCATAAGATATCGCTGCTGTAGCCAGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTGTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGTAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACACACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                                                                                                                         <  2:174067/251‑1 (MQ=255)
                                                                                                                                           tgtgCATAAGATATCGCTGCTGTAGCACGGTAATTCGTGACTTTTAGTGGCTGATTCCTGTTTCTTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGCATGATGTTCGCAACAATGCCGCGGACAAAGGATGAACTACGCGGAAGGGAACAACATTAATACTGCAATTGAATTTTTTTCACTCACTTTTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAACCGCATCATAAAAATCGGcc                                                                                                    >  2:426409/1‑250 (MQ=255)
                                                                                                                                                      tatCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCtt                                                                                        >  2:72720/1‑251 (MQ=255)
                                                                                                                                                                                   aTTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAActct                                                           <  2:260655/251‑1 (MQ=255)
                                                                                                                                                                                            gCTGATTCCTGTTTATTTGTGCAAGTGCAGTTGAGTTGTTCTGGCGGTGGCATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGACTTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCAtt                                                  >  2:179928/1‑251 (MQ=255)
                                                                                                                                                                                                     tGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACaaaa                                                                         <  2:335678/219‑1 (MQ=255)
                                                                                                                                                                                                     tGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACaaaa                                                                         >  1:335678/1‑219 (MQ=255)
                                                                                                                                                                                                      gTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATATATATTGAAATTAGTAAACTCATCATAAAAATCGGCCAAAAAATATCTTGTACGAGTTACAAAACCTATGTTAACTCTGTAGGCATTCCTGCGAAct                                        <  2:60404/251‑2 (MQ=255)
                                                                                                                                                                                                                                         ggTGGAATGATGTTCTCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTAAAAAATCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTAc     >  2:14429/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ggAATGATGTTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAg  <  1:120130/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
TCTGGTAAGCCCACCATCGTTAAGCCGGGTAGACCTTTACTGATATGTACCTCAACAGTGATCGGGGGCGCATTTACTCCCAGGGCTGCGCGGGTATGAACAATTGACAGTGACATAAGCCCTCCTTGAGTCACCATTATGTGCATAAGATATCGCTGCTGTAGCCCGCTAATTCGTGAATTTTAGTGGCTGATTCCTGTTTATTTGTGCAAGTGAAGTTGAGTTGTTCTGGCGGTGGAATGATGCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAAAAATATCTTGTACTATTTACAAAACCTATGGTAACTCTTTAGGCATTCCTTCGAACAAGATGCAAGAAAAGACAAAATGACAGCCCTTCTACGAG  >  NC_000913/3949854‑3950340

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