Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,742,755 G→C P163R (CCG→CGG)  ribC ← riboflavin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,742,7550GC100.0% 158.5 / NA 40P163R (CCG→CGG) ribCriboflavin synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (15/25);  total (15/25)

GCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742710‑1742800
                                             |                                             
gcgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:1522840/1‑50 (MQ=255)
 cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          <  1:3341440/50‑1 (MQ=255)
  gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          >  1:3980927/1‑49 (MQ=255)
    cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  1:32199/1‑50 (MQ=255)
     aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:89187/50‑1 (MQ=255)
        ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   <  1:714725/50‑1 (MQ=255)
        ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   >  1:2558745/1‑50 (MQ=255)
            ttcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGc                                <  1:1127320/49‑1 (MQ=255)
             tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               <  1:3218950/49‑1 (MQ=255)
               ttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             >  1:2053872/1‑49 (MQ=255)
                 cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           <  1:2106903/49‑1 (MQ=255)
                 cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaa                            <  1:1601376/48‑1 (MQ=255)
                 cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             <  1:4091051/47‑1 (MQ=255)
                  ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         <  1:3099459/50‑1 (MQ=255)
                  ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          <  1:2723690/49‑1 (MQ=255)
                   cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:3957559/50‑1 (MQ=255)
                   cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:4662990/50‑1 (MQ=255)
                     agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCCAAAACGCGt                      >  1:2555462/1‑50 (MQ=255)
                     agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      <  1:3534571/50‑1 (MQ=255)
                     agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      <  1:3194828/50‑1 (MQ=255)
                       agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    <  1:2998547/50‑1 (MQ=255)
                           gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:3652840/1‑50 (MQ=255)
                           gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:2386798/1‑50 (MQ=255)
                            tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:2485040/1‑49 (MQ=255)
                              cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             <  1:2938468/50‑1 (MQ=255)
                              cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             <  1:2674540/50‑1 (MQ=255)
                              cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             >  1:4350586/1‑50 (MQ=255)
                              cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             <  1:174044/50‑1 (MQ=255)
                               gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:1569505/1‑50 (MQ=255)
                               gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:2131568/1‑50 (MQ=255)
                                ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           <  1:4665200/50‑1 (MQ=255)
                                  ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:543302/49‑1 (MQ=255)
                                   cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         <  1:881459/49‑1 (MQ=255)
                                   cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGc        <  1:3894376/50‑1 (MQ=255)
                                    aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:1298871/1‑50 (MQ=255)
                                     gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:2945884/1‑49 (MQ=255)
                                     gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCg      <  1:3612190/50‑1 (MQ=255)
                                        tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:44741/49‑1 (MQ=255)
                                         ttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  >  1:867891/1‑50 (MQ=255)
                                          tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  <  1:2365468/49‑1 (MQ=255)
                                             |                                             
GCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742710‑1742800

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