Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,901 G→T P96Q (CCG→CAG)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,9010GT96.2% 77.0 / NA 26P96Q (CCG→CAG) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (16/9);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTCGCAGACGGGTTTT  >  NC_000913/1358677‑1359145
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
tgatgGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCCGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCAACGAAGGTTATCAACCACAACCACGGCGGCACTTCCGGGAAGAGGGCGGGGAGATAGAATTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCCGCCGACCATAAATCCGACGGGCTGGTTAATCCACTTTTGCGGGTAggg                                                                                                                                                                                                                              >  2:258083/1‑249 (MQ=255)
  atgGAGATATGCTCCAGTACAAACAGGGTATTGAAGTTAGCGCCCAGGTTGACGCTCTTTAGAGTCGCGGCGGTTAAAATGGCGGCGAAGGTTACGTCCCACCGCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGTCAACAAGACTTTAATCATCGGAAAAAAGAGTCAAGCCAGCAGCGATGACCAGGGGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGAGCTTAGGtata                                                                                                                                                                                                                           <  1:42724/251‑1 (MQ=255)
             cACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATAGCGAAGAAGGTTAACACCCACACCCACGGCGGCCCTTCCGGGAAGAGGGCGGGGAGATAGATTTTCGCCAAAAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCCGCCGACCATAAATCCCACGTGCTGGTTAATCGGCTTTTTCGCGTAGGGTTAGGCCGCAACCg                                                                                                                                                                                                                >  2:144081/1‑251 (MQ=255)
                          aGGGTACTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGGATAAGCCG‑AACCGGCCTCCGGAAACTg                                                                                                                                                                                                  >  1:128041/1‑251 (MQ=255)
                                ttGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTAACAACCACCCCCCCGGCGGCCCTTACGGGACGAGGGCGGAGAGATAGCTTTTCGCCAAAAAGACGTTAATCATCGGCCAAAAGAGACACTACAGCAGCGATGACCAACCGACCATAAATCCGACGTGCTGGTTAATACCCTTTTTCGCGTAGGTATAGGCCG‑AACCGGCCCCCGGAAACTGGCGAAc                                                                                                                                                                                            >  2:275586/1‑251 (MQ=255)
                                       tAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATAGCGACGAAGGTTACCACCCACCCCCACGGAGGCACTTCCGGGAAGAGGGAGGGGAGATAGATTTTTGCCAACAAGACGTTAATCATCGGCCACAAGAGACAACCCAGCAGCGGTGACCAGCCGCCCATAAATCCGACGTGCTGGTTAATCCACTTTTGCGCGTAGGTATTAGCCG‑CACCGGCCTCCGGAAACCGACGAACCAGtatg                                                                                                                                                                                     >  2:82368/1‑248 (MQ=255)
                                           gACCAGGTTGACGCCCTTCAGATTCGCGGCGGATAAAATGGCGACGAAGGTTAACACCCACACCCCCGGCGGGACTTCCGGGAAGAGGGCGGAGAGATCGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGGTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCAAACCAGGTTggcgc                                                                                                                                                                                 >  1:150574/1‑247 (MQ=255)
                                                            tcggATTCGCGGGGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTACTCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAGCTGGCGAACCAGTTTGCCGTAGCTGATAGAGGTaaa                                                                                                                                                                 <  2:93867/247‑1 (MQ=255)
                                                                     cggcggTTAAAATGGCGACGAAGGGTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCCGTTTGCCGTAGCTGATAGCGGTAAACAGCACAccc                                                                                                                                                       >  2:192322/1‑251 (MQ=255)
                                                                            aaaaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAACCACCGGAAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcg                                                                                                                                                <  1:39790/250‑1 (MQ=255)
                                                                                        gagggtgaccaccCACACCCGCGGAGGCACTTCGGGGAAGAGGGGGGAGAGATAGATTTTCGCAAACAAGACGTTAAGCAGCGGAAAAAAGGGACAATCCAGGAGGGAGGTCCAGCGGAGGATAAATCCGATGTGGGGGTTAAGCGACTTGTGCGCGTGGGTATAGGCCG‑AACCGGCGGCTGGAAACTGGCGACCCAGTGTGCTGTGTGTGATATGGTTAAACAGCACACCCGCCAGCGCCAGCAAATAgg                                                                                                                                    <  2:41706/244‑1 (MQ=255)
                                                                                                            ccggcggCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACTTTAATCACCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGCCTTTTGCGCGTAGTTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCg                                                                                                                <  1:146644/250‑1 (MQ=255)
                                                                                                                        ccGGGAAGAGGGCGGAGATATAGATTTTCGCCAACAAGACTTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGCCCAGCGGACCATAAACCCGACGTGCTGGTTAATCGACTTTTGCGGGTAGTAAAAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCCCCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGgag                                                                                                    <  1:416531/251‑1 (MQ=255)
                                                                                                                          gggAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCGGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCCCCCCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGAc                                                                                                  <  1:448906/251‑1 (MQ=255)
                                                                                                                                     ggAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAg                                                                                       <  1:192322/251‑1 (MQ=255)
                                                                                                                                                 tttCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACgg                                                                           >  1:451159/1‑251 (MQ=255)
                                                                                                                                                   ccgcaaacaaGACTTTAATCACCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTGGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCATCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGt                                                                          <  1:305764/245‑1 (MQ=255)
                                                                                                                                                       caacaaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGCCGCTAATGCCGGAGACAATACCAAAAGAATCAAATACGGTCATCg                                                                     >  2:159940/1‑251 (MQ=255)
                                                                                                                                                        aacaaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCgg                                                                    >  1:420752/1‑251 (MQ=255)
                                                                                                                                                            aGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTCGGTATAGGCCG‑CACCGACCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTg                                                                >  1:290441/1‑251 (MQ=255)
                                                                                                                                                                               aaaGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCcatc                                              >  2:29037/1‑250 (MQ=255)
                                                                                                                                                                                       aaTCCAGCAGCGAGGACCAGCCGACCATAAATCCGACGGGCTGGTTAATCGTCTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCaccacc                                     <  1:29037/251‑1 (MQ=255)
                                                                                                                                                                                        aTCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑CACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTACTCCCGGAGACAATACCACAAGTATCAAAACCGGTCATCGGCGTGAGATAGGCAAGACCCCTCATCACCCCCt                                    >  2:344692/1‑251 (MQ=255)
                                                                                                                                                                                                  gATGACCAGCCGACCATAAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAAttt                          >  1:164788/1‑251 (MQ=255)
                                                                                                                                                                                                       acAGCCGACCATCAATCCGACGTGCTGGTTAATCGACTTTTGCGCGTAGGTATCGGCCG‑AACCGGCCTCCGGAAACTGGCGTACCAGTTTGCCGTAGCTGATAGCGGTAACCAGCACACCAGCCAGCGCCAGCAACTAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACACTACCACAAGTATCAAATACGGTCATAGACGTGAGATAGGCCAGACCCATCATCACCACCTGCCCCACTTTCAggg                     >  2:374004/2‑249 (MQ=255)
                                                                                                                                                                                                                          aCGTGCTGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTTGCAGACGGGtttt  >  2:343442/1‑251 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
TGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCG‑AACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCAGACCCATCATCACCACCTGCCACAATTTCAGTGATTTTCGCAGACGGGTTTT  >  NC_000913/1358677‑1359145

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