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% 328.9 / NA 81P163R (CCG→CGG) ribCriboflavin synthase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (43/38);  total (43/38)

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800
                                               |                                             
gTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGa                                             <  1:1655718/50‑1 (MQ=255)
  gcgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:3182575/1‑50 (MQ=255)
   cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:4928427/1‑49 (MQ=255)
   cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:1472528/1‑49 (MQ=255)
    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                         <  1:206902/50‑1 (MQ=255)
      cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  1:1139137/1‑50 (MQ=255)
      cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  1:5249977/1‑50 (MQ=255)
       aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      >  1:3029442/1‑50 (MQ=255)
       aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:6012275/50‑1 (MQ=255)
       aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:1767831/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     <  1:2316586/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     <  1:5187781/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     <  1:4616581/50‑1 (MQ=255)
         gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATgg                                    <  1:1210841/50‑1 (MQ=255)
         gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATgg                                    <  1:3331733/50‑1 (MQ=255)
          ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   <  1:4038466/50‑1 (MQ=255)
           tttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  >  1:3198911/1‑50 (MQ=255)
            ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  <  1:4118673/49‑1 (MQ=255)
            ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  <  1:1691796/49‑1 (MQ=255)
             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 <  1:3264424/49‑1 (MQ=255)
             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 <  1:201785/49‑1 (MQ=255)
             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 <  1:4150716/49‑1 (MQ=255)
               tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               <  1:4291698/49‑1 (MQ=255)
                cttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaa                              <  1:1010278/49‑1 (MQ=255)
                cttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaa                              <  1:5478374/49‑1 (MQ=255)
                  tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           >  1:5441981/1‑50 (MQ=255)
                  tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           >  1:202584/1‑50 (MQ=255)
                  tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           >  1:6026765/1‑50 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           <  1:4792170/49‑1 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         >  1:355493/1‑50 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         >  1:6746223/1‑50 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          <  1:6355150/49‑1 (MQ=255)
                      aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       >  1:2986382/1‑50 (MQ=255)
                      aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       <  1:43552/50‑1 (MQ=255)
                       agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  1:3685549/1‑50 (MQ=255)
                       agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  1:2410301/1‑50 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     <  1:1905206/50‑1 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     >  1:1341442/1‑50 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     >  1:799824/1‑50 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     >  1:1204830/1‑50 (MQ=255)
                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    <  1:2744096/50‑1 (MQ=255)
                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    <  1:2136298/50‑1 (MQ=255)
                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    <  1:423719/50‑1 (MQ=255)
                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   <  1:3042902/50‑1 (MQ=255)
                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   <  1:560499/50‑1 (MQ=255)
                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   >  1:4378719/1‑50 (MQ=255)
                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   >  1:231585/1‑50 (MQ=255)
                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   >  1:3497195/1‑50 (MQ=255)
                           tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  >  1:6398827/1‑50 (MQ=255)
                           tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  <  1:1365630/50‑1 (MQ=255)
                           tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  <  1:5107031/50‑1 (MQ=255)
                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:4256872/1‑50 (MQ=255)
                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:6438680/1‑50 (MQ=255)
                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                <  1:3180413/50‑1 (MQ=255)
                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:2640613/1‑50 (MQ=255)
                              tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:2162903/1‑49 (MQ=255)
                              tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:706494/1‑49 (MQ=255)
                              tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:1357972/1‑49 (MQ=255)
                                cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGa              <  1:820734/49‑1 (MQ=255)
                                cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGa              <  1:2021136/49‑1 (MQ=255)
                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:4587984/1‑50 (MQ=255)
                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:5995071/1‑50 (MQ=255)
                                  ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           >  1:3524434/1‑50 (MQ=255)
                                   tCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          >  1:4947928/1‑50 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:1884921/49‑1 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:3387202/49‑1 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:5903140/1‑50 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:1506420/1‑50 (MQ=255)
                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGc        <  1:3719387/50‑1 (MQ=255)
                                      aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:4083100/1‑50 (MQ=255)
                                      aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       <  1:6570009/50‑1 (MQ=255)
                                      aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       <  1:1186161/50‑1 (MQ=255)
                                      aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       <  1:339809/50‑1 (MQ=255)
                                       gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:4717602/1‑49 (MQ=255)
                                        tgtTTCCCGAATTAAATGGACGCAACAACGCGTTGGCGTGACTTCGCCGa     >  1:6596114/1‑50 (MQ=255)
                                        tgtTTCCCGAATTAAATGGACGCAAAAACGTGTTGGCGTGACTTCGCCGa     >  1:2634082/1‑50 (MQ=255)
                                        tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     >  1:5849393/1‑50 (MQ=255)
                                        tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     >  1:3963067/1‑50 (MQ=255)
                                         gtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    >  1:6383406/1‑50 (MQ=255)
                                          tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACg   >  1:6020994/1‑50 (MQ=255)
                                           ttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  >  1:2019341/1‑50 (MQ=255)
                                               |                                             
GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800

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