Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,898,152 G→A A407A (GCG→GCA sdaA → L‑serine deaminase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,1520GA100.0% 122.2 / NA 30A407A (GCG→GCAsdaAL‑serine deaminase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/11);  total (19/11)

CATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGA  >  NC_000913/1898107‑1898196
                                             |                                            
cATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGa                                          >  1:1085624/1‑50 (MQ=255)
cATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGa                                          >  1:584198/1‑50 (MQ=255)
 aTTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGAt                                         <  1:2509142/50‑1 (MQ=255)
    gAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGc                                      >  1:358492/1‑50 (MQ=255)
     aGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGct                                     >  1:1746095/1‑50 (MQ=255)
     aGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGct                                     >  1:121579/1‑50 (MQ=255)
      gCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGctc                                    <  1:205489/50‑1 (MQ=255)
       cGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGctct                                   >  1:2292999/1‑50 (MQ=255)
          aaTGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTgcg                                >  1:98073/1‑50 (MQ=255)
           aTGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTgcgc                               >  1:1905327/1‑50 (MQ=255)
           aTGCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTgcgc                               >  1:2073428/1‑50 (MQ=255)
             gCCATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTGcgccg                             >  1:890115/1‑50 (MQ=255)
               cATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCa                           >  1:1488786/1‑50 (MQ=255)
               cATTGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCa                           >  1:2063754/1‑50 (MQ=255)
                 ttGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCAcc                         <  1:762723/50‑1 (MQ=255)
                 ttGCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCAcc                         <  1:1398200/50‑1 (MQ=255)
                   gCCTCTGTGAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAg                       >  1:2022728/1‑50 (MQ=255)
                     ctctGTGAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTg                     >  1:86824/1‑50 (MQ=255)
                        tgtgAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGc                    <  1:113490/48‑1 (MQ=255)
                        tgtgAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCAc                  >  1:975453/1‑50 (MQ=255)
                          tgAAGGCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCAc                  <  1:550433/48‑1 (MQ=255)
                               gCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGt           >  1:2349423/1‑50 (MQ=255)
                               gCGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGt           >  1:1729879/1‑50 (MQ=255)
                                cGATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtc          <  1:1706120/50‑1 (MQ=255)
                                 gATTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtct         >  1:2769864/1‑50 (MQ=255)
                                  aTTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc        <  1:429475/50‑1 (MQ=255)
                                  aTTAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGtctc        <  1:2399047/50‑1 (MQ=255)
                                   ttAACGCCGCACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCg       <  1:2033683/50‑1 (MQ=255)
                                       cgccgcACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTgg   <  1:2796750/50‑1 (MQ=255)
                                        gccgcACGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGa  >  1:581898/1‑50 (MQ=255)
                                             |                                            
CATTGAGCGTAATGCCATTGCCTCTGTGAAGGCGATTAACGCCGCGCGGATGGCTCTGCGCCGCACCAGTGCACCGCGCGTCTCGCTGGA  >  NC_000913/1898107‑1898196

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