Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,125,056 (T)7→8 coding (691/1224 nt) wcaI ← putative colanic biosynthesis glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,125,0491.T100.0% 37.3 / NA 12I233N (ATC→AAC) wcaIputative colanic biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (6/6);  total (6/6)

CATCGCGCAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACATCGGCATCTGCAACATGCTGAAAACGG  >  NC_000913/2124972‑2125123
                                                                              |                                                                          
cATCGCGCAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAtt                                                                          <  1:334178/80‑1 (MQ=255)
 aTCGCGCAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAtttt                                                                        <  1:334180/81‑1 (MQ=255)
 aTCGCGCAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAttt                                                                         <  1:334179/80‑1 (MQ=255)
   cgcgCAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAtttttt                                                                      >  1:334220/1‑81 (MQ=255)
     cgcAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAtttttttt                                                                    <  2:334181/81‑2 (MQ=255)
      gcAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAttttttttg                                                                   <  2:699457/81‑3 (MQ=255)
        aGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGAttttttttgt                                                                  >  1:334212/1‑77 (MQ=255)
                           tAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTgg                                               >  2:334225/1‑81 (MQ=255)
                              cGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTa                                            <  2:334182/81‑1 (MQ=255)
                              cGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTa                                            >  2:334209/1‑81 (MQ=255)
                                  ttCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGa                                         <  1:334183/80‑1 (MQ=255)
                                  ttCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGa                                         <  2:334184/80‑1 (MQ=255)
                                          ccTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCa                                <  1:334185/81‑1 (MQ=255)
                                          ccTGCTTTTCACCAATATTGCCGGAGTAAAGAAAGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATc                                 >  1:334213/1‑80 (MQ=255)
                                              cTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACAt                            >  2:334217/1‑81 (MQ=255)
                                                   cACCAATATTGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACATCGGc                        <  1:334186/80‑1 (MQ=255)
                                                           ttGCCGGAGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACATCGGCATCTGCAAc               <  2:334187/81‑1 (MQ=255)
                                                                  aGTAAAGAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACATCGGCATCTGCAACATGCTg         >  1:334216/1‑80 (MQ=255)
                                                                        gAATGATTTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACATCGGCATCTGCAACATGCTGAAAACgg  >  1:334223/1‑81 (MQ=255)
                                                                              |                                                                          
CATCGCGCAGACGATCGGCAGCTTCAATAACGTTTTCCAGCCCCTGCTTTTCACCAATATTGCCGGAGTAAAGAATGATTTTTTTGTTATCCGGCAGGTCAAGCTGGTTACGAAGGGCATCAACATCGGCATCTGCAACATGCTGAAAACGG  >  NC_000913/2124972‑2125123

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