Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,001,486 2 bp→GC coding (316‑317/903 nt) rarD ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,001,4860CG100.0% 278.0 / NA 103R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (44/59);  total (44/59)
*NC_0009134,001,4870GC100.0% 280.1 / NA 103R106A (CGG→GCG) rarDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (44/59);  total (44/59)

Only 100 of 103 total aligned reads displayed.

CCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAG  >  NC_000913/4001420‑4001555
                                                                  ||                                                                    
cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                              <  1:925915/73‑1 (MQ=255)
cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                              <  2:967686/73‑1 (MQ=255)
cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                              >  1:4987431/4‑76 (MQ=255)
cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                              >  1:4921237/4‑76 (MQ=255)
cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                              <  1:4793014/73‑1 (MQ=255)
cgctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa                                                              >  1:4453968/4‑76 (MQ=255)
 gctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCGTCCAGCat                                                             <  2:4350867/74‑1 (MQ=255)
  ctCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata                                                            >  2:5096314/2‑76 (MQ=255)
   tCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCatat                                                           >  1:1033117/1‑76 (MQ=255)
    cGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATg                                                          <  2:4139313/76‑1 (MQ=255)
     gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATCTgg                                                         >  2:681988/1‑76 (MQ=255)
     gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                         <  1:3998079/76‑1 (MQ=255)
     gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                         >  1:5527241/1‑76 (MQ=255)
     gCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg                                                         <  1:5209465/76‑1 (MQ=255)
      ccGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGt                                                        <  1:5247023/76‑1 (MQ=255)
       cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTg                                                       <  2:3146218/76‑1 (MQ=255)
       cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTg                                                       <  2:5669952/76‑1 (MQ=255)
       cGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTg                                                       >  2:3355728/1‑76 (MQ=255)
        gAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa                                                      >  1:665782/1‑76 (MQ=255)
        gAGGAAAAACACCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa                                                      <  2:220730/76‑1 (MQ=255)
         agaaaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt                                                     <  1:2722164/73‑1 (MQ=255)
          ggAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt                                                    <  1:1913847/76‑1 (MQ=255)
            aaacTCATCCCCAGCACAATGTGCAGCGCGGGGTTCATAAAGTAACCACGGCTCGCTTCCAGCATGTGGTGATTGt                                                  <  2:2233867/76‑1 (MQ=255)
            aaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt                                                  >  1:2161594/1‑76 (MQ=255)
            aaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt                                                  >  2:901140/1‑76 (MQ=255)
            aaaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt                                                  >  2:3240097/1‑76 (MQ=255)
             aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt                                                 >  2:512244/1‑76 (MQ=255)
             aaaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGtt                                                 <  2:4160422/76‑1 (MQ=255)
              aaTCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                <  2:2583348/76‑1 (MQ=255)
                tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                              >  2:4809347/1‑76 (MQ=255)
                tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                              >  1:3460450/1‑76 (MQ=255)
                tcatcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAc                                              >  1:2702798/1‑76 (MQ=255)
                 catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                             <  2:3475399/76‑1 (MQ=255)
                 catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                             >  2:4384045/1‑76 (MQ=255)
                 catcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCAcc                                             <  2:3085990/76‑1 (MQ=255)
                   tcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                           <  1:2626226/76‑1 (MQ=255)
                   tcCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGc                                           <  2:3940061/76‑1 (MQ=255)
                     ctcAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                         <  1:1728856/74‑1 (MQ=255)
                     cccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGccc                                         <  2:3201076/76‑1 (MQ=255)
                      ccAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCa                                        <  1:4678607/76‑1 (MQ=255)
                        cgCACAATGGTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                      <  2:5982/75‑1 (MQ=255)
                        aGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                      >  2:2558834/1‑76 (MQ=255)
                        aGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGa                                      <  1:3916839/76‑1 (MQ=255)
                         gCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGAt                                     >  1:3801503/1‑76 (MQ=255)
                         gCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGACTGTTCACCGCCCAGAt                                     >  2:5608584/1‑76 (MQ=255)
                          cacaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATa                                    <  1:3006969/76‑1 (MQ=255)
                           acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                   >  1:3103191/1‑76 (MQ=255)
                           acaATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaa                                   <  2:758935/76‑1 (MQ=255)
                            taaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATaaa                                  >  2:928424/2‑76 (MQ=255)
                             aaTGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAAc                                 <  2:1506611/76‑1 (MQ=255)
                               tGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAg                               >  1:2770242/1‑76 (MQ=255)
                                gTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGt                              <  1:4847059/76‑1 (MQ=255)
                                 ttCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTa                             >  2:267641/1‑76 (MQ=255)
                                  tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                            <  1:993913/76‑1 (MQ=255)
                                  tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                            >  2:435993/1‑76 (MQ=255)
                                  tCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAg                            >  2:2887763/1‑76 (MQ=255)
                                   cACCAGCGGGTTAATAAAGGAACCAAGGCGCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGc                           >  1:3359825/1‑76 (MQ=255)
                                    aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCCCCGCCCAGATAAACACTAGcc                          >  2:5141813/1‑76 (MQ=255)
                                    aCCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGcc                          >  1:2513063/1‑76 (MQ=255)
                                      cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                <  2:3591342/52‑1 (MQ=38)
                                      cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTc                                                >  1:3591342/1‑52 (MQ=38)
                                      cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                        <  2:1858146/76‑1 (MQ=255)
                                      cAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAg                        <  1:3423180/76‑1 (MQ=255)
                                        gCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGtt                      >  2:2121595/1‑76 (MQ=255)
                                          gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                    >  1:5809430/1‑76 (MQ=255)
                                          gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                    >  1:1674551/1‑76 (MQ=255)
                                          gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                    >  2:5294107/1‑76 (MQ=255)
                                          gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                    <  1:4196620/76‑1 (MQ=255)
                                          gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                    <  1:3952030/76‑1 (MQ=255)
                                          gggTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGc                    >  1:3742139/1‑76 (MQ=255)
                                            gTTAATAAAGTAACCACGGCTCGCTTCCCGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcct                  <  2:2351581/76‑2 (MQ=255)
                                            gTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGcca                  >  1:4830580/1‑76 (MQ=255)
                                              taataaAGTATCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                <  1:2500527/76‑1 (MQ=255)
                                              taataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                >  1:1273799/1‑76 (MQ=255)
                                              taataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacc                <  1:4527708/76‑1 (MQ=255)
                                               aataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGccacca               <  2:2649233/76‑1 (MQ=255)
                                                ataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCCCCaa              >  2:5708772/1‑76 (MQ=255)
                                                ataaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCaa              <  1:4847256/76‑1 (MQ=255)
                                                 taaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt             <  1:435993/76‑1 (MQ=255)
                                                 taaAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAAt             <  1:4432095/76‑1 (MQ=255)
                                                  aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc            <  1:2319569/76‑1 (MQ=255)
                                                  aaaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATc            <  1:5573272/76‑1 (MQ=255)
                                                   aaGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGGTGCCACCAATCa           >  2:1305927/1‑76 (MQ=255)
                                                   aaGTAACCAAGGCTCGCTTCCAGCATATGGGGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCa           <  2:171010/76‑1 (MQ=255)
                                                    aGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAg          <  1:4809348/76‑1 (MQ=255)
                                                     gTACCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGc         <  2:401285/76‑1 (MQ=255)
                                                      tAACCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCa        <  1:924650/76‑1 (MQ=255)
                                                        aCCAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACg      <  2:5521799/76‑1 (MQ=255)
                                                         ccAGGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg     <  2:2155493/76‑1 (MQ=255)
                                                         ccAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg     >  2:5027221/1‑76 (MQ=255)
                                                         ccAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg     <  1:3355728/76‑1 (MQ=255)
                                                         ccAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg     >  1:5032095/1‑76 (MQ=255)
                                                         acAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACgg     <  2:2405608/75‑1 (MQ=255)
                                                          cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCCCGGc    >  2:625721/1‑76 (MQ=255)
                                                          cAAGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCAAGATAAACAGAACCCAGTTGCCACCAATCAGCACGGc    >  1:3351285/1‑76 (MQ=255)
                                                           aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa   <  1:1939493/76‑1 (MQ=255)
                                                           aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa   <  2:4806351/76‑1 (MQ=255)
                                                           aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa   <  1:3519567/76‑1 (MQ=255)
                                                           aaGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCa   <  2:5070360/76‑1 (MQ=255)
                                                            aGGCTCGCTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCag  <  1:5788049/76‑1 (MQ=255)
                                                                  ||                                                                    
CCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTTCACCGCCCAGATAAACAGTAGCCAGTTGCCACCAATCAGCACGGCAG  >  NC_000913/4001420‑4001555

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: