Predicted mutation
evidence seq id position mutation annotation gene description
RA NODE_4 17,237 A→G intergenic (+1023/+707) chaA → / ← kdsA sodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NODE_417,2370AG100.0% 102.4 / NA 48intergenic (+1023/+707)chaA/kdsAsodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/23);  total (25/23)

CGATCCTGGCGGGAATTATTACCGCAGCGATTGTCAGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCGAGC  >  NODE_4/17202‑17381
                                   |                                                                                                                                                
cGATCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGGAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGCCCATTTTTGTCTTTATCa                                >  1:1887473/1‑150 (MQ=2)
cGATCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCa                                >  1:342754/1‑150 (MQ=11)
 gATCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCAc                               >  2:3055049/1‑150 (MQ=14)
 gATCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCAc                               >  1:2583035/1‑150 (MQ=14)
  aTCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACt                              <  2:880685/150‑1 (MQ=17)
  aTCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACt                              >  2:3608454/1‑150 (MQ=17)
  aTCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTCATCACt                              >  1:2985376/1‑150 (MQ=17)
  aTCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCCGCCATTTTTGTCTTTATCACt                              >  2:1729759/1‑150 (MQ=17)
   tCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACtt                             <  2:3846619/150‑1 (MQ=17)
   tCCTGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACtt                             >  2:3732752/1‑150 (MQ=17)
    ccTGGCGGGATTTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTc                            <  1:2308558/150‑1 (MQ=17)
     cTGGCGGGAATTATTACAGCAGCGATTGTCGGGTGGTGCCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGACCATTTTTGTCTTTATCACTTcc                           <  1:522962/150‑1 (MQ=11)
      tGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGTAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCca                          >  1:2154662/1‑150 (MQ=17)
      tGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCca                          <  2:1206632/150‑1 (MQ=17)
      tGGCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCca                          >  1:217342/1‑150 (MQ=17)
       ggCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCcac                         <  2:2063625/150‑1 (MQ=17)
       ggCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCcac                         >  1:2153201/1‑150 (MQ=17)
        gCGGGAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCcaca                        <  2:3602027/150‑1 (MQ=17)
          gggAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAg                      >  1:2468681/1‑150 (MQ=17)
           ggAATTATTGCCGCGGCGATTGTCGGCGGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGt                     <  1:3338997/150‑1 (MQ=255)
           ggAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGt                     <  1:1301975/150‑1 (MQ=17)
            gAATTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTa                    >  1:3969491/1‑150 (MQ=17)
             aaTTATTACCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTaa                   <  1:2259764/150‑1 (MQ=17)
              attattACTGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTaaa                  <  1:2123823/150‑1 (MQ=17)
               ttattaCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAat                 <  1:1623588/150‑1 (MQ=17)
                tattaCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAata                <  2:3354069/150‑1 (MQ=17)
                 attaCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAatat               <  2:1503689/150‑1 (MQ=17)
                 attaCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGAAAatat               >  1:3195732/1‑150 (MQ=17)
                  ttaCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATAtt              <  1:1301254/150‑1 (MQ=17)
                    aCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTgg            >  2:1840515/1‑150 (MQ=17)
                    aCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTgg            >  1:125165/1‑150 (MQ=17)
                    aCCGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGCCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTgg            <  2:247539/150‑1 (MQ=17)
                     ccGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACGGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGc           <  1:4025794/150‑1 (MQ=17)
                     ccGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGc           >  1:3663179/1‑150 (MQ=17)
                      cGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAAGCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCa          <  1:244497/150‑1 (MQ=17)
                      cGCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCa          <  2:1476096/150‑1 (MQ=17)
                       gCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCaa         <  1:1417046/150‑1 (MQ=17)
                       gCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCACGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCaa         >  2:1052829/1‑150 (MQ=17)
                       gCAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCCGCCGCCGGGGGAACAGCCTGCGGCCATTTTTGGCTTTTTCACTTCCACAAGGAAATATTGGCaa         >  1:2082513/1‑150 (MQ=11)
                        cAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAAc        <  1:648356/150‑1 (MQ=17)
                        cAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAAc        >  1:2728369/1‑150 (MQ=17)
                        cAGCGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAAc        <  1:2034033/150‑1 (MQ=17)
                           cGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCg     >  1:3936228/1‑150 (MQ=17)
                           cGATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCg     >  2:749267/1‑150 (MQ=17)
                            gATTGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCGa    >  1:1147272/1‑150 (MQ=17)
                             aTTGTCGGCTGGTGGCGTACCCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGCAAATATTGGCAACCCGAg   >  2:2729573/1‑150 (MQ=12)
                              ttGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCGAGc  <  1:3160306/150‑1 (MQ=17)
                              ttGTCGGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCGAGc  >  2:142717/1‑150 (MQ=17)
                                   |                                                                                                                                                
CGATCCTGGCGGGAATTATTACCGCAGCGATTGTCAGCTGGTGGCGTAACCGGAAGTAAGGTGAAAACGGGCGTCAGGCTGCCGTGATGGCAATATGCGCCCGAAACCGGCCGCAGGGGTAACAGCCTGCGGCCATTTTTGTCTTTATCACTTCCACAAGTAAATATTGGCAACCCGAGC  >  NODE_4/17202‑17381

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 28 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑