Predicted mutation
evidence seq id position mutation annotation gene description
RA NODE_4 17,634 Δ1 bp intergenic (+1420/+310) 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,6320A.94.6% 87.3 / 4.1 37intergenic (+1418/+312)chaA/kdsAsodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (19/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

ATTGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17483‑17662
                                                                                                                                                     |                               
aTTGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGa                                 <  2:21156/150‑1 (MQ=37)
  tGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaac                               >  2:2270265/1‑149 (MQ=37)
   gAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaacc                              <  2:3663179/150‑3 (MQ=37)
   gAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaacc                              <  2:2154662/150‑3 (MQ=37)
    aGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaacca                             >  2:1965790/1‑147 (MQ=255)
     gTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGTGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCTGCACGCTGAGAGGTACAACCTGACAACACAACGGGGGCTAATCTTGaaccag                            >  1:2370844/1‑146 (MQ=11)
      tCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑a                         <  1:944600/150‑1 (MQ=255)
      tCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑a                         <  2:3936228/150‑1 (MQ=255)
       cGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aa                        >  1:1078749/1‑149 (MQ=255)
        gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAGCAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:306855/150‑3 (MQ=255)
        gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       >  1:3647451/1‑148 (MQ=255)
        gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:2555976/150‑3 (MQ=255)
        gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ac                        >  2:4025794/1‑150 (MQ=34)
           gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    >  1:3280517/1‑145 (MQ=255)
           gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    >  2:3486419/1‑145 (MQ=255)
            taTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaatt                   <  2:1703440/150‑7 (MQ=255)
            taTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaatt                   <  2:3195732/150‑7 (MQ=255)
             aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  <  2:4035274/150‑8 (MQ=11)
             aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  >  2:3527877/1‑143 (MQ=11)
               ccACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCaa                 <  1:3457915/150‑1 (MQ=11)
                cacaCCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAg                <  2:2728369/150‑1 (MQ=11)
                  cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              >  2:1358555/1‑150 (MQ=11)
                  cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              >  1:559896/1‑150 (MQ=255)
                   acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAAt             <  1:1286041/150‑1 (MQ=11)
                    ccAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATa            <  2:1849679/150‑1 (MQ=11)
                     cAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAg           >  1:1446893/1‑150 (MQ=11)
                      aacccGGCGTTATCCGGTGAGGCGCTATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          <  2:264901/147‑1 (MQ=11)
                      aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          >  2:103813/1‑150 (MQ=11)
                      aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          >  2:2062752/1‑150 (MQ=11)
                       aGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGcc         <  1:2585235/150‑1 (MQ=11)
                       aGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGcc         >  1:3921703/1‑150 (MQ=11)
                        gCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCg        <  2:137627/150‑1 (MQ=11)
                         ccGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGc       <  2:395294/150‑1 (MQ=11)
                          cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTTAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCa      >  1:889873/1‑150 (MQ=11)
                          cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGGTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑TACCAGAACAATTCAAGAATAACCGCa      >  1:1136621/1‑150 (MQ=11)
                           ggCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑CACCAGAACAATTCAAGAATAGCCGCaa     >  1:617605/1‑150 (MQ=11)
                            gCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaaa    >  2:2476387/1‑150 (MQ=11)
                             cGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGGGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaaaa   >  1:2842371/1‑150 (MQ=11)
                              gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  2:317786/1‑150 (MQ=11)
                              gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  1:73119/1‑150 (MQ=11)
                              gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  <  2:682607/150‑1 (MQ=11)
                              gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAt  <  2:1042174/150‑1 (MQ=11)
                              gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCCCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑TACCAGAACAATTCAAGAATAGCCGCAAAAt  >  1:2447925/1‑150 (MQ=11)
                                                                                                                                                     |                               
ATTGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17483‑17662

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