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.8% 33.5 / ‑2.1 24intergenic (+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 (13/10);  total (14/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17490‑17662
                                                                                                                                                     |                        
cGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aa                        <  1:2010962/150‑2 (MQ=255)
cGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aa                        >  2:2005677/1‑149 (MQ=255)
cGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aa                        <  1:3091511/150‑2 (MQ=255)
 gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:644711/150‑3 (MQ=255)
 gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTCTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:2834397/150‑3 (MQ=255)
   agtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaa                     >  2:2583821/1‑146 (MQ=255)
    gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    <  1:1151293/150‑6 (MQ=255)
    gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    <  2:4090401/150‑6 (MQ=255)
    gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    <  1:3426887/150‑6 (MQ=255)
     taTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaatt                   >  2:2413055/1‑144 (MQ=255)
     taTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaatt                   <  1:4113144/150‑7 (MQ=255)
      aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  <  1:971736/150‑8 (MQ=11)
      aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  <  2:1773792/150‑8 (MQ=11)
       tCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAAAGGGGGCTAATCTTGAACCCGA‑ACAATTCa                  >  1:3065397/1‑150 (MQ=11)
        ccACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCaa                 <  1:1130503/150‑1 (MQ=11)
         cacaCCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAg                <  2:578117/150‑1 (MQ=11)
         cacaCCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAg                <  2:1512866/150‑1 (MQ=11)
           cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑GACCAGAACAATTCAAGaa              >  1:3533714/1‑150 (MQ=11)
           cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGCTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              >  2:819095/1‑150 (MQ=11)
            acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑ATCCAGAACAATTCAAGAAt             >  2:2170604/1‑150 (MQ=11)
            acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAAt             <  1:867584/150‑1 (MQ=11)
            acCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAAt             <  1:1283116/150‑1 (MQ=11)
             ccAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATa            <  2:2847464/150‑1 (MQ=11)
             ccAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCCCACGCTGAGAGGTACAACCTGACAACAGAACGGGGGCTAAACTTG‑AACCAGAACAATTCAAGAATa            >  1:1776116/1‑150 (MQ=11)
              cAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAg           >  1:3061570/1‑150 (MQ=11)
              cAAGCCGGCGTTATCCCGTGAGGCGCATTGTTGCGGGGGCTTTATCCCTGGTGGCATTGGGTGCTGGAAAGAGAAAACCCCCGCACGTTGACAGGTACAACCTGACAACACAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAg           >  1:166518/1‑150 (MQ=11)
               aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGGTAATCTTG‑AACCAGAACAATTCAAGAATAGc          >  1:3229157/1‑150 (MQ=11)
               aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          >  2:3976911/1‑150 (MQ=11)
                  ccGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGc       <  2:511945/150‑1 (MQ=11)
                    ggCGTTATCCGGTGAGGCGGAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaa     >  1:2661719/1‑150 (MQ=11)
                    ggCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaa     >  2:3848852/1‑150 (MQ=11)
                    ggCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaa     <  2:747957/150‑1 (MQ=11)
                     gCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCaaa    <  1:3302701/150‑1 (MQ=11)
                       gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  2:437154/1‑150 (MQ=11)
                       gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  2:1046420/1‑150 (MQ=11)
                       gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  <  1:215077/150‑1 (MQ=11)
                       gTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGCCGCAAAAt  >  1:4137925/1‑150 (MQ=11)
                                                                                                                                                     |                        
CGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17490‑17662

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