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

TGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGT  >  NC_000913/1742709‑1742801
                                              |                                              
tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                             >  1:4906603/1‑50 (MQ=255)
tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                             <  1:2055242/50‑1 (MQ=255)
tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                             <  1:3282568/50‑1 (MQ=255)
  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                           <  1:748527/50‑1 (MQ=255)
  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                           <  1:1960957/50‑1 (MQ=255)
  cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                            >  1:3793889/1‑49 (MQ=255)
   gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                           >  1:4936557/1‑49 (MQ=255)
   gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                           >  1:613265/1‑49 (MQ=255)
    ccAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                          >  1:4423448/1‑49 (MQ=255)
     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                        <  1:4290851/50‑1 (MQ=255)
     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                        <  1:1899073/50‑1 (MQ=255)
     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                        <  1:4043105/50‑1 (MQ=255)
     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                        >  1:797135/1‑50 (MQ=255)
     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                        >  1:3896540/1‑50 (MQ=255)
     cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                        >  1:1010221/1‑50 (MQ=255)
      aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                       <  1:315965/50‑1 (MQ=255)
      aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                       >  1:5291964/1‑50 (MQ=255)
         ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                    <  1:2601220/50‑1 (MQ=255)
          tttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                   >  1:3842288/1‑50 (MQ=255)
           ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                   <  1:3149215/49‑1 (MQ=255)
           ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                  >  1:1210651/1‑50 (MQ=255)
            tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                  <  1:791764/49‑1 (MQ=255)
             ttcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGc                                 >  1:2779497/1‑49 (MQ=255)
              tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                                <  1:235698/49‑1 (MQ=255)
              tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                                <  1:3825030/49‑1 (MQ=255)
              tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                                >  1:3443949/1‑49 (MQ=255)
              tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                                <  1:5568112/49‑1 (MQ=255)
              tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                                <  1:4345592/49‑1 (MQ=255)
                ttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                              >  1:1349990/1‑49 (MQ=255)
                 tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                            >  1:2877217/1‑50 (MQ=255)
                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                              <  1:3967197/47‑1 (MQ=255)
                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           >  1:1954659/1‑50 (MQ=255)
                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           >  1:3344635/1‑50 (MQ=255)
                  cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           >  1:3622867/1‑50 (MQ=255)
                   ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                          >  1:3401704/1‑50 (MQ=255)
                   ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                          >  1:4474626/1‑50 (MQ=255)
                   ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                           <  1:411049/49‑1 (MQ=255)
                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                        >  1:1351788/1‑50 (MQ=255)
                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                        <  1:2786672/50‑1 (MQ=255)
                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                        >  1:2739240/1‑50 (MQ=255)
                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                        >  1:4831145/1‑50 (MQ=255)
                     aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                        >  1:2393923/1‑50 (MQ=255)
                      agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                       >  1:5549205/1‑50 (MQ=255)
                      agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                       <  1:402757/50‑1 (MQ=255)
                      agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                       >  1:5184213/1‑50 (MQ=255)
                       gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                      <  1:963058/50‑1 (MQ=255)
                       gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                      <  1:1283161/50‑1 (MQ=255)
                       gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                      >  1:1639566/1‑50 (MQ=255)
                        agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                     >  1:1324408/1‑50 (MQ=255)
                        agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                     >  1:2814094/1‑50 (MQ=255)
                        agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                     <  1:2313785/50‑1 (MQ=255)
                         gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                    <  1:3960478/50‑1 (MQ=255)
                          tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                   <  1:4928162/50‑1 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                  <  1:1228594/49‑1 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                  <  1:4034907/49‑1 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:3470640/1‑50 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:1699385/1‑50 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:239281/1‑50 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:2562662/1‑50 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:4497817/1‑50 (MQ=255)
                            gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                 >  1:2243760/1‑50 (MQ=255)
                             tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTg                <  1:4438634/50‑1 (MQ=255)
                             tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTg                <  1:4446850/50‑1 (MQ=255)
                               cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGa               <  1:2046045/49‑1 (MQ=255)
                                gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt             >  1:5473885/1‑50 (MQ=255)
                                gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt             >  1:1076791/1‑50 (MQ=255)
                                 ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt            <  1:4822440/50‑1 (MQ=255)
                                 ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt            <  1:3883659/50‑1 (MQ=255)
                                  tCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           >  1:4936642/1‑50 (MQ=255)
                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           <  1:5140413/49‑1 (MQ=255)
                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           <  1:4685068/49‑1 (MQ=255)
                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           <  1:4659839/49‑1 (MQ=255)
                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc           <  1:3071989/49‑1 (MQ=255)
                                   ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg          >  1:4315830/1‑50 (MQ=255)
                                     aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc        <  1:2202369/50‑1 (MQ=255)
                                     aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc        >  1:4003467/1‑50 (MQ=255)
                                         tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc     <  1:3226417/49‑1 (MQ=255)
                                         tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACg    <  1:3898654/50‑1 (MQ=255)
                                           tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGt  >  1:2496269/1‑50 (MQ=255)
                                           tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGt  >  1:457432/1‑50 (MQ=255)
                                              |                                              
TGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGGT  >  NC_000913/1742709‑1742801

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