Predicted mutation
evidence seq id position mutation annotation gene description
RA NODE_4 17,639 +A intergenic (+1425/+305) 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,6381.A95.7% 31.7 / ‑2.1 23intergenic (+1424/+306)chaA/kdsAsodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base A (8/14);  total (9/14)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17491‑17662
                                                                                                                                                    |                        
gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:2267059/150‑3 (MQ=255)
gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       >  2:3961047/1‑148 (MQ=255)
   gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGGTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGGACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    >  1:1170490/1‑145 (MQ=11)
   gtaTCCACACCAACTCGGCGTTATGCGGTGAGGCGCAATGTTGCGGGGGCTTTATCACTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    <  1:58847/150‑6 (MQ=11)
     aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  <  1:3650706/150‑8 (MQ=11)
     aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  >  1:4051866/1‑143 (MQ=11)
     aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGGTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  >  2:2166753/1‑143 (MQ=11)
       ccACACCAAGCCGGCGTTATCCGGTGAGGTGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCaa                 <  1:3205228/150‑1 (MQ=11)
        cccacCAAGCCGGCGTTATGGGGTGAGGCGCAATGTTGCGGGGGCATTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAg                <  1:444518/148‑1 (MQ=11)
        cacaCCAAGCCGGCGTTATCCGGTCAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGATGAGAGGTACAACCTGACGACAGAACGGGGGCTAATCTTGAACCAAA‑ACAAATCAAg                >  1:2383337/1‑150 (MQ=11)
          cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              <  1:955401/150‑1 (MQ=11)
          cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              <  1:3473510/150‑1 (MQ=11)
           acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAAt             >  1:600592/1‑150 (MQ=11)
           acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAAt             <  2:3361245/150‑1 (MQ=11)
            ccAAGCCGGCGTTATCCGGTGGGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATa            <  2:116946/150‑1 (MQ=11)
            ccAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATa            <  2:1282743/150‑1 (MQ=11)
             cAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAg           <  1:1504070/150‑1 (MQ=11)
             cAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGGTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAAATCAAGAATAg           >  1:2925825/1‑150 (MQ=11)
              aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          <  1:1860323/150‑1 (MQ=11)
              aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          <  2:2753101/150‑1 (MQ=11)
               ccgacttgTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGcc         <  1:2232093/143‑1 (MQ=11)
                gCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCg        <  2:580379/150‑1 (MQ=11)
                  cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCa      <  1:486615/150‑1 (MQ=11)
                  cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCa      >  1:2399348/1‑150 (MQ=11)
                  cGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCCCa      >  1:824617/1‑150 (MQ=11)
                      gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGccaagtcgg  >  2:2232093/1‑143 (MQ=11)
                      gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  2:87288/1‑150 (MQ=11)
                      gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  2:1042486/1‑150 (MQ=11)
                      gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  <  1:2787782/150‑1 (MQ=11)
                      gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCCCACGTTGGGGGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑GTCCAGAACAATTCAAGAATAGCCGCAAAAt  >  1:599727/1‑150 (MQ=11)
                                                                                                                                                    |                        
GTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17491‑17662

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