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

CGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742711‑1742800
                                            |                                             
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          <  1:4806374/50‑1 (MQ=255)
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           <  1:2817949/49‑1 (MQ=255)
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:2314096/1‑49 (MQ=255)
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           >  1:1864851/1‑49 (MQ=255)
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           <  1:648844/49‑1 (MQ=255)
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           <  1:4217449/49‑1 (MQ=255)
cgcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAt                                           <  1:1441591/49‑1 (MQ=255)
 gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          >  1:266383/1‑49 (MQ=255)
 gcCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAAtt                                          >  1:1823490/1‑49 (MQ=255)
  ccAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaa                                        <  1:5492115/50‑1 (MQ=255)
  ccAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTa                                         >  1:5367507/1‑49 (MQ=255)
   cAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTaaa                                       <  1:3859550/50‑1 (MQ=255)
    aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:143962/50‑1 (MQ=255)
    aaGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAAt                                      <  1:6002661/50‑1 (MQ=255)
       ttttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   >  1:5963749/1‑50 (MQ=255)
       ttttttCTTCCCAAGAGACGTGCGTTCCAGTGTTTCCCGAATTAAATGGa                                   >  1:4361699/1‑50 (MQ=255)
        tttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  >  1:3844750/1‑50 (MQ=255)
         ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  <  1:6304546/49‑1 (MQ=255)
         ttttCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGAc                                  <  1:5399325/49‑1 (MQ=255)
           ttcttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGc                                <  1:4524195/49‑1 (MQ=255)
             cttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaa                              <  1:1953503/49‑1 (MQ=255)
              ttcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             >  1:5968265/1‑49 (MQ=255)
               tcCCAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           >  1:4755245/1‑50 (MQ=255)
                cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaaaa                           <  1:2904598/49‑1 (MQ=255)
                cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCaaa                             <  1:6782665/47‑1 (MQ=255)
                cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  1:3470022/1‑50 (MQ=255)
                cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  1:1355648/1‑50 (MQ=255)
                cccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          >  1:507722/1‑50 (MQ=255)
                 ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         <  1:2999839/50‑1 (MQ=255)
                 ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcg                         >  1:4695086/1‑50 (MQ=255)
                 ccAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAc                          <  1:4843226/49‑1 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:2390699/50‑1 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:7101697/50‑1 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:4815644/50‑1 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        >  1:4034635/1‑50 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        >  1:4030363/1‑50 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:3712495/50‑1 (MQ=255)
                  cAAGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgc                        <  1:2356584/50‑1 (MQ=255)
                   aaGAGTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAAcgcg                       >  1:2251075/1‑50 (MQ=255)
                    agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  1:1744400/1‑50 (MQ=255)
                    agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      <  1:4911103/50‑1 (MQ=255)
                    agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      >  1:2326741/1‑50 (MQ=255)
                    agagTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGt                      <  1:1560267/50‑1 (MQ=255)
                      agTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTg                    <  1:1720907/50‑1 (MQ=255)
                       gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   >  1:135141/1‑50 (MQ=255)
                       gTCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTgg                   <  1:248536/50‑1 (MQ=255)
                        tCGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGc                  <  1:1791879/50‑1 (MQ=255)
                         cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  1:855751/1‑50 (MQ=255)
                         cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  1:4100558/1‑50 (MQ=255)
                         cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  1:2976126/1‑50 (MQ=255)
                         cGTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 >  1:2826355/1‑50 (MQ=255)
                          gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCg                 <  1:3270953/49‑1 (MQ=255)
                          gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:1366458/1‑50 (MQ=255)
                          gTGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGt                >  1:609101/1‑50 (MQ=255)
                           tGCGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTg               <  1:711279/50‑1 (MQ=255)
                             cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             <  1:124258/50‑1 (MQ=255)
                             cGTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGAc             >  1:4875324/1‑50 (MQ=255)
                              gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:4617167/1‑50 (MQ=255)
                              gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:455066/1‑50 (MQ=255)
                              gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            <  1:4486683/50‑1 (MQ=255)
                              gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            <  1:2506373/50‑1 (MQ=255)
                              gTTCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACt            >  1:5540102/1‑50 (MQ=255)
                               ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           <  1:1380138/50‑1 (MQ=255)
                               ttCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACtt           <  1:3351225/50‑1 (MQ=255)
                                tCCAGTGTTTCCCGAATTAAATGGACGCACAAACGCGTTGGCGTGACTTc          >  1:1054255/1‑50 (MQ=255)
                                tCCAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          >  1:6737373/1‑50 (MQ=255)
                                 ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:5890402/49‑1 (MQ=255)
                                 ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTc          <  1:4119530/49‑1 (MQ=255)
                                 ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:5695716/1‑50 (MQ=255)
                                 ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:5120069/1‑50 (MQ=255)
                                 ccAGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         >  1:445107/1‑50 (MQ=255)
                                  cAGTGTTGCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCg         <  1:4530723/49‑1 (MQ=255)
                                   aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:1036368/1‑50 (MQ=255)
                                   aGTGTTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGcc       >  1:4476331/1‑50 (MQ=255)
                                    gtgtTTCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCg      <  1:3586611/50‑1 (MQ=255)
                                       tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:5544659/49‑1 (MQ=255)
                                       tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:202002/49‑1 (MQ=255)
                                       tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGAc    <  1:6681132/49‑1 (MQ=255)
                                       tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACg   >  1:13694/1‑50 (MQ=255)
                                       tttCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACg   <  1:6870438/50‑1 (MQ=255)
                                         tCCCGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACgg  <  1:6151463/49‑1 (MQ=255)
                                            |                                             
CGCCAAGTTTTTTCTTCCCAAGAGTCGTGCGTTCCAGTGTTTCCGGAATTAAATGGACGCAAAAACGCGTTGGCGTGACTTCGCCGACGG  >  NC_000913/1742711‑1742800

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