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.4% 79.7 / ‑5.8 28Y110C (TAT→TGT) puuPputrescine:H(+) symporter PuuP
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/16);  minor base A (1/0);  total (12/16)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCC  >  NC_000913/1358613‑1359092
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
cAGACGGTGCCAACGCCGTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCTACGAAGGTTACCTCCCACACCCTCGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATc                                                                                                                                                                                                                                        <  1:57670/251‑1 (MQ=255)
  gACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCca                                                                                                                                                                                                                                      <  1:320068/251‑1 (MQ=255)
    cGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCcagc                                                                                                                                                                                                                                    <  1:229498/251‑1 (MQ=255)
          cAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGAt                                                                                                                                                                                                                              <  2:255314/251‑1 (MQ=255)
             ccccTTCTCCTCTATGCCGTCACTGAACCACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAATAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGAc                                                                                                                                                                                                                           <  1:28337/249‑1 (MQ=255)
              gCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGGAGAGGGGGGAG‑GACTAGATTTTCGCCCACAAGCCGCTAATCATCGGCAAAAAAAGACAATCCCGCCGAGATGAc                                                                                                                                                                                                                           >  2:66642/1‑250 (MQ=255)
                     ccagtATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGTTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGTCAAAAAGAGCCAATCCAGCAGCGATGACCAGCCg                                                                                                                                                                                                                     <  2:149653/245‑1 (MQ=255)
                          aTGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAAGTTACCCCCCCCACCCAAGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTTGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAAACCAGCAGCGATGACCAGCCGACCATa                                                                                                                                                                                                               >  1:48188/1‑250 (MQ=255)
                                       aCTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCgg                                                                                                                                                                                                 <  1:179579/251‑1 (MQ=255)
                                                       acacCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCCCGGAGGCACTTCCGGGGAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACACGCCGACCATAAATCCGACGTGCGGGTTaacagccctttgg                                                                                                                                                                                 >  2:218615/1‑240 (MQ=255)
                                                       acacCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGATAAAAAGGCGACGAAAGTTACCACCCCCACCCACGGCGGCACTTCCGGGGAGAGGGCGGGGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGGGGGTTAATCGACTTTTgg                                                                                                                                                                                 >  1:238679/1‑250 (MQ=255)
                                                                       gaTTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCg                                                                                                                                                                 <  2:70389/251‑1 (MQ=255)
                                                                                          aGGGTCTTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTTAAATAGCGACGGAGGGTACCACCCACACCCCCCGCGGCACCTCCGGGGAGAGGGCGGAGAGA‑AAGATTTTCGCCAAAAAGACGTTAATCATTGGCAAAAAGAGAAAATCCAGCAGCGCAGACCAGCCGGCCATAAATCCGAAGAGCGGGATAATCGACTTTTGCGCGTAGGGATAGTCCGaaacacgcgccggaaaccg                                                                                                                                              >  2:170544/1‑234 (MQ=255)
                                                                                                    aGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCCCGGCGGCACTTCCGGGAAGAGGGCGGGGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAAACA‑Gt                                                                                                                                   >  1:286552/1‑251 (MQ=255)
                                                                                                                  ttGACGCTCTTCAGATTCGCGGCGGTTAAAATGGGGACGAAAGTTCCCACCCACACCCCCGGCGGCCCTTTCGGGGAGAGGGCGGGGAGA‑TAGATTTTCGCCAACAAAACGTTAATCATCGGCACAAAGAGACAATCCAGCAGCGATGACCAGCCACCCATAAAATCGACGAGAGGGTTAATCGACTTTTGCGCGTTGGTATAGGCCGAACCGGCCTCCAGAAACTCGCGCACCA‑GATTGTCGTCGCTGAt                                                                                                                     >  1:150926/1‑251 (MQ=255)
                                                                                                                         ccgtCAGATTCTCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACGCCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTGTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAAGCGACGTTTGCGCTTAGGTATAGGCCGAACCTGCCTTCGGAAACTTGCGAACCA‑GTTTGCCGTATCTGATAGCGGTa                                                                                                              <  2:185737/248‑1 (MQ=255)
                                                                                                                             cAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACa                                                                                                          <  2:150926/251‑1 (MQ=255)
                                                                                                                               gATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGc                                                                                                        <  1:49314/251‑1 (MQ=255)
                                                                                                                                 ttCGCGGCGGGTAAAATGGCGACGAAGGTTGCCACCCACACCCACGGCGGCACTTCCGGGAAGATGGCGGAGAGA‑TAGATTTTCGCCAACAAGGCGTTAATCATCGGCAAAAATAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGc                                                                                                                                                                                                   <  1:148031/159‑1 (MQ=255)
                                                                                                                                    gcggcgGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACAcc                                                                                                   <  1:247154/251‑1 (MQ=255)
                                                                                                                                     cggcggTTAAAATTGCGACGAAGGTTACCACCCACACCCACGGCGGCCCTTCCGGGGAGAGGGCGGAGAGA‑TAGAATTTCGCCAAAAAGACGGTAATCATCGGCAAAAAGAGACAATCCAGCAGCGCTGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGGCGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGTACCAAGTTT‑CCGTAgacggcagagggaaacaccacaccc                                                                                                  >  2:165021/1‑227 (MQ=255)
                                                                                                                                                         aaGGTTAGCACACACACCCACGGCGGCAGTTCAGGACAGAGGGCGGAGAGA‑TAGCTTTTCGCCAACAAGACGTTAATCATCGGAAAAAAGAGACAATCCAGCAGCGTTGAACAGACGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGACCCGGCCTCCGGCAACCGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCACCGCCAGCAAATAgg                                                                               >  2:207117/1‑250 (MQ=255)
                                                                                                                                                            gTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTCCGCCAACAAGGCGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGGCGTAGCGGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGc                                                                           <  1:170544/251‑1 (MQ=255)
                                                                                                                                                              tACCACCCACACCCACGGCGGCCATTCAGGGAAGGGGGTGGAGAAA‑AAGATTTTCGCCACAAAGACGTAAACCATCGGAAAAAAGAGACAACCCAGCAGCGATGCCCAGCCGACCATAAAGCCGCCGTGCGGGTTAATCGTCTTGGGCGCGTAGGTATATGCCGAACCGGCCTCAGGAAACTGGCGAACCG‑GTTTGCCGAAGCTGATGGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGcc                                                                          <  2:281249/250‑1 (MQ=255)
                                                                                                                                                                                             aaGAGGGAGGAGAGA‑TAGATTTTCGCCAAAAAGACGTAAATCATCGCCAAAAAGAGACAACCCAGCAGCGATGACCAGCCGACCATAAATCCTACGTGCGGGTTAATCGCCTTTGGCGCGTAGGAATAGGCCGAACCGGCCTCCGGAAACTGGAGACCCA‑GGTTGCCGTAGCTGATCGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAAt                                          <  1:63426/251‑1 (MQ=255)
                                                                                                                                                                                                     ggAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGCGACAATACCaaa                                    >  2:148024/1‑249 (MQ=255)
                                                                                                                                                                                                                        caacaaGACGTTAATCATCGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCg                >  1:244116/1‑251 (MQ=255)
                                                                                                                                                                                                                                        atcGGCAAAAAGAGACAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTAATCGGCGTGAGATAGGcc  >  1:16702/1‑249 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                          
CAGACGGTGCCAACGCCTTCTCCTTTATGCAGTCCCTGAACTACCAGGAAGATAAACACCACCATGATGGAGATTTGCACCAGTACAAACAGGGTATTGAAGTTAGCGACCAGGTTGACGCTCTTCAGATTCGCGGCGGTTAAAATGGCGACGAAGGTTACCACCCACACCCACGGCGGCACTTCCGGGAAGAGGGCGGAGAGA‑TAGATTTTCGCCAACAAGACGTTAATCATCGGCAAAAAGAGATAATCCAGCAGCGATGACCAGCCGACCATAAATCCGACGTGCGGGTTAATCGACTTTTGCGCGTAGGTATAGGCCGAACCGGCCTCCGGAAACTGGCGAACCA‑GTTTGCCGTAGCTGATAGCGGTAAACAGCACACCCGCCAGCGCCAGCAAATAGGACGCCGGAACGTGACCGTCGCTAATGCCGGAGACAATACCAAAAGTATCAAATACGGTCATCGGCGTGAGATAGGCC  >  NC_000913/1358613‑1359092

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