Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 136,062 T→G K135Q (AAG→CAG)  speE ← spermidine synthase (putrescine aminopropyltransferase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913136,0620TG100.0% 199.2 / NA 53K135Q (AAG→CAG) speEspermidine synthase (putrescine aminopropyltransferase)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/31);  total (22/31)

GATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTTAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTCGCCGCC  >  NC_000913/135915‑136209
                                                                                                                                                   |                                                                                                                                                   
gATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGaaa                                                                                                                                                  <  2:493983/151‑1 (MQ=255)
   tCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAgc                                                                                                                                                <  2:274505/150‑1 (MQ=255)
       gCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGa                                                                                                                                           <  2:175412/151‑1 (MQ=255)
           cGTTTGCAGCCTTCATAAAATGCCGAAGTGAATAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtcgt                                                                                                                                       >  1:310953/1‑151 (MQ=255)
           cGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtcgt                                                                                                                                       >  1:475607/1‑151 (MQ=255)
            gTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtcgtc                                                                                                                                      >  2:837927/1‑151 (MQ=255)
              ttGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGAtcgtcg                                                                                                                                     <  1:105677/150‑1 (MQ=255)
                       tCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTAcc                                                                                                                             <  1:496405/149‑1 (MQ=255)
                       tCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCgg                                                                                                                           >  2:130536/1‑151 (MQ=255)
                       tCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCgg                                                                                                                           >  1:683064/1‑151 (MQ=255)
                        cATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGc                                                                                                                          <  1:110231/151‑1 (MQ=255)
                            aaaTGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGtt                                                                                                                       >  2:257945/1‑150 (MQ=255)
                             aaTGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTAt                                                                                                                     <  1:799726/151‑1 (MQ=255)
                               tGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATgg                                                                                                                   >  1:720532/1‑151 (MQ=255)
                                gCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATgg                                                                                                                   <  1:294164/150‑1 (MQ=255)
                                             ggCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTg                                                                                                     <  1:199866/151‑1 (MQ=255)
                                                tttCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACg                                                                                                  >  1:608670/1‑151 (MQ=255)
                                                   cGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGc                                                                                                >  2:323144/1‑150 (MQ=255)
                                                   cGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGc                                                                                                <  2:354504/150‑1 (MQ=255)
                                                        ggACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACg                                                                                          <  1:901214/151‑1 (MQ=255)
                                                          aCCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGAt                                                                                        <  2:760609/151‑1 (MQ=255)
                                                          aCCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGAt                                                                                        <  2:407456/151‑1 (MQ=255)
                                                                   aTCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACAccc                                                                               <  1:890102/151‑1 (MQ=255)
                                                                    tCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCg                                                                              >  1:246286/1‑151 (MQ=255)
                                                                            gTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGAtt                                                                      <  2:83957/151‑1 (MQ=255)
                                                                              cGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTc                                                                    >  1:670385/1‑151 (MQ=255)
                                                                              cGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTc                                                                    <  2:819763/151‑1 (MQ=255)
                                                                                 agaTAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTccac                                                                 >  1:820986/1‑151 (MQ=255)
                                                                                 agaTAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTccac                                                                 <  2:578217/151‑1 (MQ=255)
                                                                                      aTGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCg                                                            >  1:562633/1‑151 (MQ=255)
                                                                                        gACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTg                                                          >  2:915869/1‑151 (MQ=255)
                                                                                        gACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTg                                                          >  1:301851/1‑151 (MQ=255)
                                                                                        gACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTg                                                          <  1:10469/151‑1 (MQ=255)
                                                                                           aTCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGAtt                                                       >  2:828835/1‑151 (MQ=255)
                                                                                               aaGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACt                                                   >  2:75122/1‑151 (MQ=255)
                                                                                                                   aaCGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCGGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGtt                               <  2:194156/151‑1 (MQ=255)
                                                                                                                    aCGAAATTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGtt                               >  2:692389/1‑150 (MQ=255)
                                                                                                                       aaaTTGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACt                            <  1:168561/150‑1 (MQ=255)
                                                                                                                           tGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCAc                        <  2:651796/150‑1 (MQ=255)
                                                                                                                           tGACGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCAc                        <  2:683083/150‑1 (MQ=255)
                                                                                                                             aCGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACgca                     <  2:436176/151‑1 (MQ=255)
                                                                                                                             aCGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACgc                      >  1:580505/1‑150 (MQ=255)
                                                                                                                              cGCCATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACgcag                    <  1:815357/151‑1 (MQ=255)
                                                                                                                                ccATCGTCGATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACgcagc                   <  1:539444/150‑1 (MQ=255)
                                                                                                                                    cgtcgATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGGCACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATgg               <  1:62599/150‑1 (MQ=255)
                                                                                                                                     gtcgATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGc              >  2:465436/1‑150 (MQ=255)
                                                                                                                                     gtcgATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGc              >  2:868268/1‑150 (MQ=255)
                                                                                                                                        gATCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCg          <  1:761038/151‑1 (MQ=255)
                                                                                                                                         aTCACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGt         <  2:232449/151‑1 (MQ=255)
                                                                                                                                           cACCAGCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTc        <  2:874698/150‑1 (MQ=255)
                                                                                                                                             ccAGCTGAAAGCGTGAATCGTCGTACCTACCGGCGTTATGGCTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGCGTTACCTCACGGAGCAAGCCAc            >  1:312672/1‑144 (MQ=255)
                                                                                                                                                gCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTcgccgcc  <  1:267526/151‑1 (MQ=255)
                                                                                                                                                gCTGAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTcgccgcc  <  2:327776/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GATTCAGGCAACGTTTGCAGCCTTCATAAAATGCCGAAGTGAAAAGGCTTTCGCCGGGACCGATAGGATCGGTGCAGTCGGAGATAATGACATCAAAGGTCTGGCTGGTTTGATTAACGAAATTGACGCCATCGTCGATCACCAGCTTAAAGCGCGGATCGTCGTAGCTACCGGCGTTATGGTTGGGTAGATACTGACGGCAGAACGATACGACACCCGCATCGATTTCCACCATCGTGATTGACTCAACGTTTTTATGTCGGGTTACTTCACGCAGCATGGCACCGTCGCCGCC  >  NC_000913/135915‑136209

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