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.96.7% 50.0 / ‑1.5 30intergenic (+1418/+312)chaA/kdsAsodium‑potassium/proton antiporter ChaA/3‑deoxy‑8‑phosphooctulonate synthase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (12/17);  total (13/17)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17485‑17662
                                                                                                                                                   |                               
tGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaac                               >  2:469521/1‑149 (MQ=37)
tGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaac                               <  1:2946724/150‑2 (MQ=37)
 gAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGaacc                              >  1:2002375/1‑148 (MQ=37)
      gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       >  2:3961047/1‑148 (MQ=255)
      gtagtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aac                       <  2:2267059/150‑3 (MQ=255)
         gtaTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGGTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGGACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    >  1:1170490/1‑145 (MQ=11)
         gtaTCCACACCAACTCGGCGTTATGCGGTGAGGCGCAATGTTGCGGGGGCTTTATCACTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaat                    <  1:58847/150‑6 (MQ=11)
           aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTGatccagaacaattc                    >  2:4054921/1‑137 (MQ=11)
           aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  >  1:4051866/1‑143 (MQ=11)
           aTCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAG‑aacaattc                  <  1:3650706/150‑8 (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:3473510/150‑1 (MQ=11)
                cacCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGaa              <  1:955401/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          <  2:2753101/150‑1 (MQ=11)
                    aaGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTACAACCTGACAACAGAACGGGGGCTAATCTTG‑AACCAGAACAATTCAAGAATAGc          <  1:1860323/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)
                                                                                                                                                   |                               
TGAGTCGTAGTATCCACACCAAGCCGGCGTTATCCGGTGAGGCGCAATGTTGCGGGGGCTTTATCCCTGGTGGCATTGGTTGCTGGAAAGAGAAAACCCCCGCACGTTGAGAGGTAAAACCTGACAACAGAACGGGGGCTAATCTTGAAACCAG‑ACACTTCAAGAATAGCCGCAAAAT  >  NODE_4/17485‑17662

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