Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 698,541 G→T R213S (CGT→AGT)  asnB ← asparagine synthetase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913698,5410GT100.0% 110.5 / NA 31R213S (CGT→AGT) asnBasparagine synthetase B
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (10/21);  total (10/21)

AGACCACCAGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGAT  >  NC_000913/698465‑698616
                                                                            |                                                                           
aGACCACCAGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGc                                                                         >  2:135866/1‑81 (MQ=255)
  accaccAGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTc                                                                       <  1:135826/81‑1 (MQ=255)
    caccaGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGt                                                                     <  2:135827/81‑1 (MQ=255)
    caccaGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGt                                                                     >  1:135859/1‑81 (MQ=255)
    caccaGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGt                                                                     >  2:135869/1‑81 (MQ=255)
     accaGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGtt                                                                    <  1:135828/81‑1 (MQ=255)
     accaGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGtt                                                                    <  2:135830/81‑1 (MQ=255)
     accaGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGtt                                                                    <  2:135829/81‑1 (MQ=255)
        aGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGttttt                                                                 <  1:135831/81‑1 (MQ=255)
              cagcaCACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGt                                                           <  2:135833/81‑1 (MQ=255)
              cagcaCACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGt                                                           <  1:135832/81‑1 (MQ=255)
               agcaCACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGt                                                           >  2:135863/1‑80 (MQ=255)
                  acacCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACg                                                       >  1:135867/1‑81 (MQ=255)
                    acCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGtt                                                     <  2:135834/81‑1 (MQ=255)
                          aGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCttt                                               <  1:135835/81‑1 (MQ=255)
                             cacaTCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCAc                                            <  1:135836/81‑1 (MQ=255)
                              acatcagacatcagATGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCAcc                                           <  1:135837/81‑1 (MQ=255)
                                atcagacatcagaTGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCg                                          >  2:135868/1‑80 (MQ=255)
                                  cagacatcagaTGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCAt                                       >  2:135871/1‑81 (MQ=255)
                                   agtcatcagaTGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATc                                      <  2:135838/78‑1 (MQ=255)
                                    gacatcagaTGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCg                                     >  2:135865/1‑81 (MQ=255)
                                        tcagaTGGCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGt                                 <  2:135839/81‑1 (MQ=255)
                                           gaTGGCTTTTAACTGAATTTTCTAGTGCCTGAGTAAGCGCGTTTTTGTCGGTAAGGTTATCTTTAACCGCATCGTAGCCGa                              <  2:135840/81‑1 (MQ=255)
                                               gCTTTTAACTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAAcc                           <  1:135841/80‑1 (MQ=255)
                                                      aCTGAATCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCga                   >  1:135878/1‑81 (MQ=255)
                                                          aaTCTTCCAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCgatg                 <  1:135842/79‑1 (MQ=255)
                                                                ccAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTaa          <  2:135844/80‑1 (MQ=255)
                                                                ccAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTaa          <  1:135843/80‑1 (MQ=255)
                                                                 cAGTGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAg         <  1:135845/80‑1 (MQ=255)
                                                                    tGCCTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACg     <  1:135846/81‑1 (MQ=255)
                                                                       cTGACTCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGAt  >  1:135874/1‑81 (MQ=255)
                                                                            |                                                                           
AGACCACCAGAAAGCAGCACACCGTAAGGCACATCAGACATCAGATGGCTTTTAACTGAATCTTCCAGTGCCTGACGCAGCTCGTTTTTGTCGGTCACGTTATCTTTCACCGCATCGTAGTCGAACCAGTCGCGATGATAGTAAGAACGGAT  >  NC_000913/698465‑698616

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