Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,539,484 T→A N34Y (AAT→TAT)  cysM ← cysteine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,539,4840TA100.0% 104.2 / NA 29N34Y (AAT→TAT) cysMcysteine synthase B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (11/18);  total (11/18)

TTTAATTTCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATTGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCACCAGAGGCGTA  >  NC_000913/2539407‑2539554
                                                                             |                                                                      
tttAATTTCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGcc                                                                     <  1:462458/81‑1 (MQ=255)
   aaTTTCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCtt                                                                   <  1:462459/80‑1 (MQ=255)
   aaTTTCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTc                                                                  >  2:462502/1‑81 (MQ=255)
       tCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAg                                                               <  1:462461/80‑1 (MQ=255)
       tCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAg                                                               <  2:462460/80‑1 (MQ=255)
        ccccGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGtt                                                             >  1:462507/1‑81 (MQ=255)
        ccccGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGt                                                              <  2:462462/80‑1 (MQ=255)
          ccGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGtttt                                                           <  2:462464/81‑1 (MQ=255)
          ccGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGtttt                                                           <  2:462463/81‑1 (MQ=255)
           cgcgCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGtttt                                                           >  2:462496/1‑80 (MQ=255)
            gcgcTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTaa                                                         >  2:462504/1‑81 (MQ=255)
              gcTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAAc                                                        >  1:462503/1‑80 (MQ=255)
                  ttCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACcac                                                     <  2:462465/79‑1 (MQ=255)
                     cGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACtt                                                 <  1:462466/80‑1 (MQ=255)
                     cGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTc                                                <  1:462467/81‑1 (MQ=255)
                       ccTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCAc                                              <  2:462468/81‑1 (MQ=255)
                           gacgaTCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGcc                                          >  1:462499/1‑81 (MQ=255)
                           gacgaTCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGcc                                          <  2:462469/81‑1 (MQ=255)
                            acgaTCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCg                                         <  2:462470/81‑1 (MQ=255)
                                            ccGCACGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTc                         >  1:462509/1‑81 (MQ=255)
                                                aCGATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTg                     <  1:462471/81‑1 (MQ=255)
                                                  gATCTTTCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGc                    >  1:462517/1‑80 (MQ=255)
                                                        tCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACtt               <  2:462473/79‑1 (MQ=255)
                                                        tCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTc              <  2:462472/80‑1 (MQ=255)
                                                        tCACCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCa             >  2:462515/1‑81 (MQ=255)
                                                          aCCGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCAcc           <  2:462474/81‑1 (MQ=255)
                                                           ccGAACCTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCACCa          <  1:462475/81‑1 (MQ=255)
                                                                 cTGCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCACCAGAGGCg    >  1:462518/1‑81 (MQ=255)
                                                                   gCCGGGTTATAGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCACCAGAGGCGTa  >  2:462511/1‑81 (MQ=255)
                                                                             |                                                                      
TTTAATTTCCCCGCGCTTTTCCGCCTCGACGATCATCGAAAGTGCCGCACGATCTTTCACCGAACCTGCCGGGTTATTGCCTTCCAGTTTTAACCACACTTCACTGCCGTTATCCGGCCCCATTCGCTGCAACTTCACCAGAGGCGTA  >  NC_000913/2539407‑2539554

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