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

GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800
                                               |                                             
gTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGa                                             <  1:4245408/50‑1 (MQ=255)
 tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            >  1:4627671/1‑50 (MQ=255)
 tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            >  1:662534/1‑50 (MQ=255)
 tGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGaa                                            >  1:3380369/1‑50 (MQ=255)
  gcgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:640881/1‑50 (MQ=255)
   cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          <  1:41908/50‑1 (MQ=255)
   cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           <  1:4733282/49‑1 (MQ=255)
    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          >  1:2095225/1‑49 (MQ=255)
    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                         <  1:1154851/50‑1 (MQ=255)
    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                         <  1:2857224/50‑1 (MQ=255)
    gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                         <  1:742659/50‑1 (MQ=255)
      cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       >  1:5063736/1‑50 (MQ=255)
       aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:5255731/50‑1 (MQ=255)
       aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      >  1:4810198/1‑50 (MQ=255)
       aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:2547/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGTGTTCCAGTGTTTCCCGAATTAAATg                                     <  1:5605245/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     <  1:5222721/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     <  1:3193866/50‑1 (MQ=255)
        aGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     >  1:2600604/1‑50 (MQ=255)
         gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATgg                                    <  1:3738607/50‑1 (MQ=255)
         gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     >  1:1982160/1‑49 (MQ=255)
         gTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATg                                     >  1:1509530/1‑49 (MQ=255)
          ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   <  1:3192711/50‑1 (MQ=255)
           tttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  >  1:4945133/1‑50 (MQ=255)
            ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 >  1:3200854/1‑50 (MQ=255)
             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 >  1:4890893/1‑49 (MQ=255)
             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 >  1:3833716/1‑49 (MQ=255)
             tttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACg                                 >  1:412727/1‑49 (MQ=255)
              ttcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGc                                <  1:2798235/49‑1 (MQ=255)
               tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               <  1:360849/49‑1 (MQ=255)
               tcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCa                               >  1:492269/1‑49 (MQ=255)
                cttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaa                              <  1:727455/49‑1 (MQ=255)
                  tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           >  1:1406202/1‑50 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           <  1:4397963/49‑1 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaa                            <  1:2715537/48‑1 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaa                            <  1:1852890/48‑1 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             <  1:4493672/47‑1 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             <  1:2666467/47‑1 (MQ=255)
                   cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  1:771247/1‑50 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         >  1:5456360/1‑50 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         >  1:3872510/1‑50 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         <  1:3802574/50‑1 (MQ=255)
                    ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          <  1:2424963/49‑1 (MQ=255)
                     cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:776137/50‑1 (MQ=255)
                      aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       <  1:2037908/50‑1 (MQ=255)
                       agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  1:2249626/1‑50 (MQ=255)
                       agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  1:1842714/1‑50 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     <  1:1101317/50‑1 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     <  1:2312603/50‑1 (MQ=255)
                        gagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGtt                     >  1:3048047/1‑50 (MQ=255)
                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    >  1:1203131/1‑50 (MQ=255)
                         agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    >  1:2369178/1‑50 (MQ=255)
                          gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   >  1:4204014/1‑50 (MQ=255)
                            cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  1:5336910/1‑50 (MQ=255)
                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGTCGt                >  1:4717399/1‑50 (MQ=255)
                             gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:5716978/1‑50 (MQ=255)
                               gCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGa              >  1:3508537/1‑50 (MQ=255)
                                cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             >  1:4924031/1‑50 (MQ=255)
                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:5401343/1‑50 (MQ=255)
                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:3968088/1‑50 (MQ=255)
                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:3061301/1‑50 (MQ=255)
                                 gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:1199397/1‑50 (MQ=255)
                                  ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           >  1:3002199/1‑50 (MQ=255)
                                  ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           <  1:2691105/50‑1 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:1907899/49‑1 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:2849205/1‑50 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:4388873/1‑50 (MQ=255)
                                    ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:3139467/1‑50 (MQ=255)
                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         <  1:2920613/49‑1 (MQ=255)
                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         <  1:4229188/49‑1 (MQ=255)
                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGc        <  1:2684642/50‑1 (MQ=255)
                                     cAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGc        >  1:2405207/1‑50 (MQ=255)
                                      aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       <  1:5328207/50‑1 (MQ=255)
                                        tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     >  1:4649248/1‑50 (MQ=255)
                                        tgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGa     >  1:764499/1‑50 (MQ=255)
                                         gtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    >  1:3733289/1‑50 (MQ=255)
                                         gtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    >  1:4093748/1‑50 (MQ=255)
                                          tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:3675520/49‑1 (MQ=255)
                                           ttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  >  1:1304337/1‑50 (MQ=255)
                                           ttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGGGACTTCGCCGACgg  >  1:3260326/1‑50 (MQ=255)
                                            tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  <  1:1639659/49‑1 (MQ=255)
                                               |                                             
GTGCGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742708‑1742800

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: