Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 618,034 G→A R1293K (AGG→AAG)  entF → apo‑serine activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913618,0340GA100.0% 72.9 / NA 23R1293K (AGG→AAG) entFapo‑serine activating enzyme
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (11/12);  total (11/12)

GCAACTACGCTGATGCTGTGCGCCTGCTGACGACTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAGGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCATTTTCAGGTTATACCCCATGTCATCACT  >  NC_000913/617800‑618264
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
gCAACTACGCTGATGCTGTGCGCCTGCTGACGACTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATattatt                                                                                                                                                                                                                        <  1:254666/251‑1 (MQ=255)
   aCTACGCTGATGCTGTGCGCCTGCTGACGACTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATAttatttat                                                                                                                                                                                                                     >  1:346937/1‑251 (MQ=255)
     tACGCTGATGCTGTGCGCCTGCTGACGACTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGGATGAGTCCCGTACGCGCCTGGTCGACGTGGATAGCGGAGCTGGTTCTCTATCGCCAGGATTGTGCGCATGAGGATATTATCTCTCCAGGGACGTTTGAAACAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTtata                                                                                                                                                                                                                    >  2:364160/1‑250 (MQ=255)
                            tacgacTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAaataa                                                                                                                                                                                              <  2:202736/248‑1 (MQ=255)
                            gacgacTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAaataat                                                                                                                                                                                             <  1:151359/250‑1 (MQ=255)
                                                catttGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGCCGTTTGAAAAAATTTGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTAtttgtttg                                                                                                                                                                        <  1:17580/249‑1 (MQ=255)
                                                 gTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTAtttgtttg                                                                                                                                                                        >  2:342851/1‑250 (MQ=255)
                                                     gACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGAttt                                                                                                                                                                    <  2:381228/250‑1 (MQ=255)
                                                                            tgcggaatgCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGGGCGCATGTGGATATTATCTCTCCAGGGACGTGTGAAAAAAGTGGGCCGAGTATTCGCGCAACGCTAAACAAGTAAAGTAATATTATTTATAAACCCATAATTCCAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTctct                                                                                                                                             <  1:310803/242‑1 (MQ=255)
                                                                                 aaCGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTAAAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAAt                                                                                                                                        >  1:168329/1‑250 (MQ=255)
                                                                                                        aGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTgga                                                                                                                <  1:149876/251‑1 (MQ=255)
                                                                                                                 cgcnCCTGGTCGCCGTAGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATAGTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGGGGAAGGAATTTa                                                                                                       <  2:416413/251‑1 (MQ=255)
                                                                                                                              cGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCtt                                                                                          >  1:373025/1‑251 (MQ=255)
                                                                                                                                                  atctatcGTCCGGGTTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGACAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTAAAAACTTCATCTTCGGATACTTATCACCANCTg                                                                      >  2:410792/1‑251 (MQ=255)
                                                                                                                                                           cAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTg                                                              <  2:436060/250‑1 (MQ=255)
                                                                                                                                                                   atcgaatttgGATATTATCTCTACAGGGACGTTTGAAAAAATTGGGCCGATTTTTCGCGCAACGCTAAACAAGTAAATTATTTTTATTTATAAACGCATAATTACAGACAATAATTATGGGTTTTTTCTTTGCTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAATGCATTTACAAATTTCAGCTTCGGTTAATTATCAGCAACTGGTAAAGTGAGCGACGTg                                                     >  2:78223/9‑251 (MQ=255)
                                                                                                                                                                        aTGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTgg                                                    >  1:360555/1‑247 (MQ=255)
                                                                                                                                                                                          caaggACGTTTGAAAAAATTGGGCCTATTATTCGCGCAACGCTAAACAAGTAAAGTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGGTTGATTTATAGGGTTGATGAATATTTCTCTGAAATATAGTGAATGTTGCATGCAGTTGAGGGAAGTTACAAACTTCAGCTTCGGATAATTATCACCAACTCGTAAAGTGAGCGCCGTGGAGATTTCTCCCTATTGGACTCa                              <  2:433366/248‑1 (MQ=255)
                                                                                                                                                                                             ggACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAGCTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCAttt                           <  1:342851/251‑1 (MQ=255)
                                                                                                                                                                                                    tGAAAAGATTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCATTTTCAGGtt                    >  2:202495/1‑251 (MQ=255)
                                                                                                                                                                                                       aaaaaTTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATCGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCATTTTCAGGatat                  >  1:50533/1‑250 (MQ=255)
                                                                                                                                                                                                           aTTGGGCCGATTATTCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATCATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCATTTTCCGGTTCTAcccc             >  1:37619/1‑251 (MQ=255)
                                                                                                                                                                                                                       attCGCGCAACGCTAAACAAGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCATTTTCAGGTTATACCCCATGTCATCACt  <  1:405429/250‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      
GCAACTACGCTGATGCTGTGCGCCTGCTGACGACTGCTCATAGCGTACCGTTTGACGGTAAAGCGACGCTGTTTGTTGCTGAACGCACACTTCAGGAAGGTATGAGTCCCGAACGCGCCTGGTCGCCGTGGATAGCGGAGCTGGATATCTATCGTCAGGATTGTGCGCATGTGGATATTATCTCTCCAGGGACGTTTGAAAAAATTGGGCCGATTATTCGCGCAACGCTAAACAGGTAAATTAATATTATTTATAAACCCATAATTACAGAAAATAATTATGGGTTTTTTATTTGTTTGATTTATAGGTTTGATGAATATTTCTCTTAAATAGAGTGAATGTTGCATTCAGTGGAAGGAATTTACAAACTTCAGCTTCGGATAATTATCACCAACTGGTAAAGTGAGCGCCGTGGAGATTTCCCCCTATTGGACTCATTTTCAGGTTATACCCCATGTCATCACT  >  NC_000913/617800‑618264

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