Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,656,744 G→A *187* (TAG→TAA speG → spermidine N(1)‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,656,7440GA100.0% 173.6 / NA 55*187* (TAG→TAAspeGspermidine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (27/28);  total (27/28)

TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTA  >  NC_000913/1656608‑1656874
                                                                                                                                        |                                                                                                                                  
ttGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  <  2:54906/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTTATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:633063/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:140614/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  2:76731/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  2:156735/139‑1 (MQ=255)
                                 aTGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATtaataataa                                                                                                                        >  1:824806/1‑116 (MQ=255)
                                 aTGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATtaataataa                                                                                                                        <  2:824806/116‑1 (MQ=255)
                                    gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAAc                                                                                              >  2:801383/1‑139 (MQ=255)
                                    gTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGCCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAAc                                                                                              >  1:992606/1‑139 (MQ=255)
                                              tAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  >  2:901472/1‑93 (MQ=255)
                                              tAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  <  1:901472/93‑1 (MQ=255)
                                                aTGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                  <  2:229733/139‑1 (MQ=255)
                                                 tGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggcg                                                                                 <  1:179537/139‑1 (MQ=255)
                                                      ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  >  1:378893/1‑85 (MQ=255)
                                                          cATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                        <  1:801383/139‑1 (MQ=255)
                                                                    ttCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaa                                                              <  1:337800/139‑1 (MQ=255)
                                                                     tCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  2:323151/96‑1 (MQ=255)
                                                                     tCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  1:323151/1‑96 (MQ=255)
                                                                       cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  2:743625/94‑1 (MQ=255)
                                                                       cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  1:743625/1‑94 (MQ=255)
                                                                       cAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaacaa                                                           >  1:886363/1‑139 (MQ=255)
                                                                        aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAATCCCGTTCGACGGCGGGTTTTACTCCTtc                                                                   <  1:664557/130‑1 (MQ=255)
                                                                        aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtc                                                                   >  2:664557/1‑130 (MQ=255)
                                                                        aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAAt                                                          <  2:992606/139‑1 (MQ=255)
                                                                         gCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATg                                                         >  2:142087/1‑139 (MQ=255)
                                                                          cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                  <  2:91244/65‑1 (MQ=255)
                                                                          cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                        >  1:34973/1‑139 (MQ=255)
                                                                          cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                        >  1:601733/1‑139 (MQ=255)
                                                                          cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTGTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                        >  1:711678/1‑139 (MQ=255)
                                                                            tCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGAtt                                                      <  2:601733/139‑1 (MQ=255)
                                                                                    tGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACa                                              >  2:660138/1‑139 (MQ=255)
                                                                                      gCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGt                                            >  2:794063/1‑139 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggcg                                                                                 <  2:755356/98‑1 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggcg                                                                                 >  1:755356/1‑98 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                        >  1:242206/1‑139 (MQ=255)
                                                                                              caAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGa                                    <  2:711678/139‑1 (MQ=255)
                                                                                                   cacCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGc                               <  1:1074060/139‑1 (MQ=255)
                                                                                                       gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                           <  2:34973/139‑1 (MQ=255)
                                                                                                         gTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                         >  2:900782/1‑139 (MQ=255)
                                                                                                                       aaGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCGTCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAg            >  1:690713/1‑138 (MQ=255)
                                                                                                                       aaGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAg            <  2:690713/138‑1 (MQ=255)
                                                                                                                        aGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGct          <  1:142087/139‑1 (MQ=255)
                                                                                                                           cGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTg       <  1:340462/139‑1 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                         >  2:225049/1‑120 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                         <  1:225049/120‑1 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      >  2:27651/1‑139 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      >  2:824018/1‑139 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      >  2:493833/1‑139 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt      >  1:1028323/1‑139 (MQ=255)
                                                                                                                             aCCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        >  2:332491/1‑40 (MQ=255)
                                                                                                                             aCCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                        <  1:332491/40‑1 (MQ=255)
                                                                                                                              ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTAACAGCTCTGttt    >  2:642005/1‑139 (MQ=255)
                                                                                                                                gCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTa  <  1:27651/139‑1 (MQ=255)
                                                                                                                                  acaATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                         <  2:219098/114‑1 (MQ=255)
                                                                                                                                  acaATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                         >  1:219098/1‑114 (MQ=255)
                                                                                                                                        |                                                                                                                                  
TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTA  >  NC_000913/1656608‑1656874

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