Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,859 T→C Y110C (TAT→TGT)  puuP ← putrescine:H(+) symporter PuuP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,8590TC96.2% 76.5 / NA 26Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (13/12);  minor base A (1/0);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

GCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCA  >  NC_000913/1358621‑1359093
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          
gCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCGCACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAGCAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCg                                                                                                                                                                                                                                <  1:370321/251‑1 (MQ=255)
                  aggCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACTCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAATACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATaa                                                                                                                                                                                                              <  1:220464/249‑1 (MQ=255)
                      aGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGAAGAAGGTTCCCACCCACACCCACGGCGGCACCTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACCAGAAGTTAATCATCGGCAAAAAAAGACAATCCAGAAACGATGAACAGCCGACCATAAATcc                                                                                                                                                                                                          >  1:173719/1‑251 (MQ=255)
                        tCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCCGGGAGAGGGCGGAGAGATAGCTTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCg                                                                                                                                                                                                         >  1:364730/1‑250 (MQ=255)
                           cTGAAATACCCGGAAGATAAACACCACCATGCTGGAGCTTTGCACCAGTACAAACAGGGAATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGGTTCCACCCACACCCACGGGGGCACTTCCGGGAAAGGGGGGGAGAAATAGATTTTCGCCCACAAGACGCTAATCATAGGCAAAAAGAGAAAATCCAACAGCGATGACCAGCCGACCATAAATCCCACgc                                                                                                                                                                                                     >  2:244330/1‑250 (MQ=255)
                                  aCCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGt                                                                                                                                                                                               <  2:175411/250‑1 (MQ=255)
                                   ccAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTa                                                                                                                                                                                             <  1:142504/251‑1 (MQ=255)
                                         aGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAAAGGCGACGAAGGGTTACCCCCACACCCACGGCGGCACTTACGGGAAGAGGGCGGAGAGATAAATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAACGAAGACAAGCCGACCCTAAATCCGCCGTGCGGGTTAATAGAc                                                                                                                                                                                       >  2:123666/1‑251 (MQ=255)
                                                caccacCATGAGGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGTCGGCACTTCTGGGAAGAGGGCGGATAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTgcg                                                                                                                                                                                <  1:345439/251‑1 (MQ=255)
                                                    accaTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAAGGGGGGGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCCTAAAACCGACCTGCGGGTTAATCGACTTTTGCGCGt                                                                                                                                                                             >  2:302905/1‑250 (MQ=255)
                                                      caTGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAAACGTTAAACATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGGTAATCGACTTTTGCGCGTAgg                                                                                                                                                                          >  1:3547/1‑251 (MQ=255)
                                                                              aaaCAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGaa                                                                                                                                                  <  1:191939/251‑1 (MQ=255)
                                                                                                    aCCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAAACCGACGGGCGGGTTAATCGACTTTTGCGCGTAGGGATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGCTTGCCGt                                                                                                                             >  1:322258/1‑250 (MQ=255)
                                                                                                                  ctTCAGATTCGCGGCGGTTTAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATc                                                                                                                                                                                          >  2:49783/1‑175 (MQ=255)
                                                                                                                  ctTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATc                                                                                                                                                                                          <  1:49783/175‑1 (MQ=255)
                                                                                                                       gATTCGCGGCTGTTAAAATGGCGACGAAGGTTACCACCCACACCCACTGCGGCGCTTCCGGGAAGAGGGCGGATAGATATATTTTCGCCAACAAGACGTTAATCATCGGCAAAAATAGACAATCCAGCAGCGATGGCCAGCCGCCCAGAAAGCCGACGTGCGGGTTAATCGACTTTTGCGCTGAGGTATTGGCCGAACCGGTGTCCGGAAACTGTGGAATCGGTTTGCCTGAGGTGATACCGGTAAACAGc                                                                                                         <  2:28405/251‑1 (MQ=255)
                                                                                                                       gATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGc                                                                                                         <  2:327004/251‑1 (MQ=255)
                                                                                                                           cgcgGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGGTTGCCGTAGCTGATAGCGGTAAACAGcaca                                                                                                      <  2:356817/250‑1 (MQ=255)
                                                                                                                                      aaaTGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCcgccagcgcc                                                                                          >  1:266926/1‑251 (MQ=255)
                                                                                                                                             gacgaAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAAt                                                                                   >  2:185797/1‑251 (MQ=255)
                                                                                                                                                 aaGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAgg                                                                                >  1:140166/1‑250 (MQ=255)
                                                                                                                                                        accccccacccccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACGGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGAACACCCGCCAGCGCCAGCAAATAGGACGCCg                                                                          <  1:167939/246‑1 (MQ=255)
                                                                                                                                                                   ccacGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGt                                                              >  1:217508/1‑251 (MQ=255)
                                                                                                                                                                                    ggAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGACCCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACa                                             <  1:294389/251‑1 (MQ=255)
                                                                                                                                                                                                agagaTAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGt                                  >  1:245063/1‑250 (MQ=255)
                                                                                                                                                                                                                               catcGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATCGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATCCGGTCATCGGCGTGAGATAGGccc  >  2:232375/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                          
GCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGATAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCAGTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCCA  >  NC_000913/1358621‑1359093

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: