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,7440GA97.9% 146.7 / ‑6.8 48*187* (TAG→TAAspeGspermidine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (23/24);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCT  >  NC_000913/1656621‑1656866
                                                                                                                           |                                                                                                                          
gATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcg                                                                                                             >  2:1122957/1‑139 (MQ=255)
     aCGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGt                                                                                                        <  2:637543/139‑1 (MQ=255)
             tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  1:1122957/139‑1 (MQ=255)
               tattaATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTTCTGAAGCCGACCGCCCAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  1:885660/137‑1 (MQ=255)
               tattaATGGTCAATATCGTAATGCCATGCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  2:885660/1‑137 (MQ=255)
                    aTGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  2:1301898/132‑1 (MQ=255)
                    aTGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  1:1301898/1‑132 (MQ=255)
                               cGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGa                                                                              <  2:101964/139‑1 (MQ=255)
                                         ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGttt                                                                    >  1:1091087/1‑139 (MQ=255)
                                                     tatTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcatc                                                        >  2:381702/1‑139 (MQ=255)
                                                         ccAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATTGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaaca                                                    <  2:104388/139‑1 (MQ=255)
                                                           aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTaataataat                                                                                                               >  2:992834/1‑78 (MQ=255)
                                                           aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTaataataat                                                                                                               <  1:992834/78‑1 (MQ=255)
                                                           aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  2:186540/93‑1 (MQ=255)
                                                           aGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  1:186540/1‑93 (MQ=255)
                                                            gCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATg                                                 >  1:56372/1‑139 (MQ=255)
                                                             cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAtt                                                                                                                        <  2:1351658/67‑1 (MQ=255)
                                                             cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACAtt                                                                                                                        >  1:1351658/1‑67 (MQ=255)
                                                             cATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGa                                                >  2:603958/1‑139 (MQ=255)
                                                               tCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                          <  1:296190/63‑1 (MQ=255)
                                                               tCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                          >  2:296190/1‑63 (MQ=255)
                                                                cAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTa                                             >  2:101770/1‑139 (MQ=255)
                                                                  gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGc                                           >  2:225586/1‑139 (MQ=255)
                                                                             agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  1:189742/75‑1 (MQ=255)
                                                                             agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  2:189742/1‑75 (MQ=255)
                                                                             agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  2:365654/75‑1 (MQ=255)
                                                                             agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  2:444587/75‑1 (MQ=255)
                                                                             agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  1:444587/1‑75 (MQ=255)
                                                                             agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  1:365654/1‑75 (MQ=255)
                                                                               cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                          >  1:1185685/1‑63 (MQ=255)
                                                                               cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                          <  2:1185685/63‑1 (MQ=255)
                                                                               cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                          >  2:527915/1‑95 (MQ=255)
                                                                               cacaAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                          <  1:527915/95‑1 (MQ=255)
                                                                                  aaaaCACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  2:887946/1‑70 (MQ=255)
                                                                                  aaaaCACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  1:887946/70‑1 (MQ=255)
                                                                                          gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                   <  2:56372/139‑1 (MQ=255)
                                                                                          gggTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                   <  1:101770/139‑1 (MQ=255)
                                                                                            gTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                 <  1:1004138/139‑1 (MQ=255)
                                                                                               aGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTc              <  2:1091087/139‑1 (MQ=255)
                                                                                                     tcCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                   >  1:760661/1‑128 (MQ=255)
                                                                                                     tcCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                   <  2:760661/128‑1 (MQ=255)
                                                                                                      ccTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTca       <  1:225586/139‑1 (MQ=255)
                                                                                                         gAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                <  1:904520/47‑1 (MQ=255)
                                                                                                         gAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                >  2:904520/1‑47 (MQ=255)
                                                                                                                 ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                <  1:831362/71‑1 (MQ=255)
                                                                                                                 ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACt                                                                >  2:831362/1‑71 (MQ=255)
                                                                                                                 ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGct  <  2:1043182/133‑1 (MQ=255)
                                                                                                                 ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGct  >  1:1043182/1‑133 (MQ=255)
                                                                                                                           |                                                                                                                          
GATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCT  >  NC_000913/1656621‑1656866

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