Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 617,030 A→G L958L (TTA→TTG entF → apo‑serine activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913617,0300AG100.0% 74.9 / NA 23L958L (TTA→TTGentFapo‑serine activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/11);  total (12/11)

TTGGATACCAGCGCATTGCAGGCGCAGCTTCGTGAAACATTGCCACCACATATGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTACCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCTGGGGTGTGACGTGCAGGATGCCGATGCTGATTTCTTCGCG  >  NC_000913/616905‑617159
                                                                                                                             |                                                                                                                                 
ttGGATACCAGCGCATTGCAGGCGCAGCTTCGTGAAACACTGCCACCACATATGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACCTAGCGCCAGCGGAAAGATGGATCGCTAAGcgttgccgttgc                                                                                                                            >  2:144007/1‑133 (MQ=255)
  ggATACCAGCGCATTGCAGGCGCAGCTTCGTGAAACATTGCCACCACATATGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCct                                                                                                                          >  1:22369/1‑133 (MQ=255)
                     gcgcAGCTTCGTGAAACATTGTCACCACATATGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAgcgc                                                                                                       >  1:432620/1‑133 (MQ=255)
                      cgcAGCTTCGGGAAACATTGCCACCACATATGTTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGcc                                                                                                      <  2:171557/133‑1 (MQ=255)
                       gcAGCTTCGTGAAACATTGCCACCACATATGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCa                                                                                                     >  2:9688/1‑133 (MQ=255)
                                 gAAGCATTGCCACCATATACGGTAACGGAGGTCCTGCTGCAGCTGCCGCAGTTACCGCTTAGCGCCACCGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTgcgc                                                                                           <  2:196165/133‑1 (MQ=255)
                                      actGCCACCACATATCGTCCCGCTGGTTCTTCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGaaa                                                                                      <  1:144007/131‑1 (MQ=255)
                                               acaTATGATACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGt                                                                             >  1:507609/1‑133 (MQ=255)
                                                   atGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGCGCGCCGAAAGCGGTCAGTGaaa                                                                         >  1:225359/1‑133 (MQ=255)
                                                        aCCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGt                                                                             <  1:369267/124‑1 (MQ=255)
                                                        aCCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGt                                                                             >  2:369267/1‑124 (MQ=255)
                                                           ggtggtTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATc                                                                 >  1:341358/1‑133 (MQ=255)
                                                               gtTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATcgccg                                                             <  1:385725/133‑1 (MQ=255)
                                                                           cTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCtcgt                                                 <  1:9688/133‑1 (MQ=255)
                                                                             tCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCTGTGGAACGATTATCGCAGCGGGATTCtcgtcg                                               >  2:212348/1‑133 (MQ=255)
                                                                                  aGTTACCACTTAGCGCCAACGGGAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCt                                          <  1:212348/133‑1 (MQ=255)
                                                                                      aCCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCTgggg                                      >  1:432923/1‑133 (MQ=255)
                                                                                      aCCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATCATCGCCGAGGCATTCTCGTCGGTGCTgggg                                      >  1:164640/1‑133 (MQ=255)
                                                                                            tAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCTGGGGTGTGAc                                <  1:510360/133‑1 (MQ=255)
                                                                                             aGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCGTTGTCGGCGGTGCTGGTGTGTGACg                               >  1:165236/1‑133 (MQ=255)
                                                                                             aGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCt                                          >  2:292475/1‑122 (MQ=255)
                                                                                             aGCGCCAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCt                                          <  1:292475/122‑1 (MQ=255)
                                                                                                  cAACGGCAAGCTGGATCGCAAAGCCTTGCCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCTGGGGTGTGACGTGCAg                          <  2:507609/133‑1 (MQ=255)
                                                                                                      ggCAAGCTGGCTCGCGAAGCCTTGCCGTTGCCAGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCTGGGGTGTGACGTGCAGGATg                      <  2:164640/133‑1 (MQ=255)
                                                                                                                          cTTGCCGTTGCCTGATCTGAAGGCCCAAGCGCTAGGGCGTGAGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCGTTCTCGTCGTTGCTGGGGTGTGACGTGCAGGATGCCGATGCTGATTTCTTcgcg  <  2:165236/133‑1 (MQ=255)
                                                                                                                             |                                                                                                                                 
TTGGATACCAGCGCATTGCAGGCGCAGCTTCGTGAAACATTGCCACCACATATGGTACCGGTGGTTCTGCTGCAACTTCCACAGTTACCACTTAGCGCCAACGGCAAGCTGGATCGCAAAGCCTTACCGTTGCCTGAACTGAAGGCACAAGCGCCAGGGCGTGCGCCGAAAGCGGGCAGTGAAACGATTATCGCCGCGGCATTCTCGTCGTTGCTGGGGTGTGACGTGCAGGATGCCGATGCTGATTTCTTCGCG  >  NC_000913/616905‑617159

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: