Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,610,589 A→G V31A (GTT→GCT)  uxaB ← tagaturonate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,610,5890AG100.0% 101.2 / NA 29V31A (GTT→GCT) uxaBtagaturonate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/19);  total (10/19)

TCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGC  >  NC_000913/1610517‑1610661
                                                                        |                                                                        
tCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGgcg                                                                  >  1:292752/1‑81 (MQ=255)
tCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGgc                                                                   <  1:292718/80‑1 (MQ=255)
  aaTCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGgcgc                                                                 <  2:292719/80‑1 (MQ=255)
    tCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCa                                                              <  2:292720/81‑1 (MQ=255)
      ggACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAg                                                             <  2:292722/80‑1 (MQ=255)
      ggACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAg                                                             <  2:292721/80‑1 (MQ=255)
              aacGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTAcc                                                    >  1:292760/1‑81 (MQ=255)
               acGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCt                                                   <  1:292723/81‑1 (MQ=255)
               acGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCt                                                   <  2:292724/81‑1 (MQ=255)
                cGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCtt                                                  <  1:292725/81‑1 (MQ=255)
                 gACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCtt                                                  >  2:292763/1‑80 (MQ=255)
                  aCCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCa                                                <  2:292726/81‑1 (MQ=255)
                     aCGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCa                                             <  1:292727/81‑1 (MQ=255)
                          aGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACt                                         >  1:292753/1‑80 (MQ=255)
                          aGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACt                                         >  2:292750/1‑80 (MQ=255)
                                cAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGa                                   <  2:292728/80‑1 (MQ=255)
                                  gATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATAc                                <  2:292729/81‑1 (MQ=255)
                                    tCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGt                              >  2:292783/1‑81 (MQ=255)
                                       gtgtGCTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCt                           >  2:292761/1‑81 (MQ=255)
                                            cTCATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGAt                       <  1:292730/80‑1 (MQ=255)
                                              cATTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACt                    >  1:292762/1‑81 (MQ=255)
                                               aTTCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATAc                     <  2:292731/79‑1 (MQ=255)
                                                ttCAGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACtgt                  <  1:292732/81‑1 (MQ=255)
                                                   aGGAGATCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGc                <  2:292733/80‑1 (MQ=255)
                                                        aTCGATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACCGAATGATACGTTCTGGATACTGTGCACCgg           >  2:292758/1‑80 (MQ=255)
                                                           gATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGaaa       <  1:292735/81‑1 (MQ=255)
                                                           gATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGaaa       <  2:292734/81‑1 (MQ=255)
                                                           gATTTGCCAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGaa        >  1:292766/1‑80 (MQ=255)
                                                                 ccAGTCAGCAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATcgc  <  1:292736/80‑1 (MQ=255)
                                                                        |                                                                        
TCAATCGGACGAACAACGACCACGCCAGAATTCAGATCGGTGTGCTCATTCAGGAGATCGATTTGCCAGTCAACAAAGGCGCGCAGGAAGTTACCTTCACCAAACTGAATGATACGTTCTGGATACTGTGCACCGGGAAAATCGC  >  NC_000913/1610517‑1610661

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