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.6% 125.1 / NA 41*187* (TAG→TAAspeGspermidine N(1)‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/24);  minor base T (0/1);  total (16/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

TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTA  >  NC_000913/1656608‑1656879
                                                                                                                                        |                                                                                                                                       
ttGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACa                                                                                                                                       <  1:203216/139‑1 (MQ=255)
            tGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAACACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatc                                                                                                                           <  2:721822/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  2:725149/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  2:689129/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  2:464365/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  2:321158/139‑1 (MQ=255)
                          tttATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  1:619656/139‑1 (MQ=255)
                                      cAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                       >  2:635414/1‑117 (MQ=255)
                                      cAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAatcgatc                                                                                                                       <  1:635414/117‑1 (MQ=255)
                                                   ccATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  1:13966/114‑1 (MQ=255)
                                                   ccATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             >  2:13966/1‑114 (MQ=255)
                                                   ccATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCggg                                                                                    >  2:149300/1‑139 (MQ=255)
                                                   ccATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCggg                                                                                    >  2:538423/1‑139 (MQ=255)
                                                     aTTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCgggct                                                                                  >  1:68478/1‑137 (MQ=255)
                                                      ttCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGttt                                                                                 <  1:538423/139‑1 (MQ=255)
                                                       tCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGtttt                                                                                <  1:397508/139‑1 (MQ=255)
                                                                tataTTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtca                                                                       >  1:112002/1‑139 (MQ=255)
                                                                      ccAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATcaaca                                                                 <  1:305982/139‑1 (MQ=255)
                                                                               gTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             <  1:452892/86‑1 (MQ=255)
                                                                               gTATCTGGCAGAGCACAAAAAACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGa                                                                                                             >  2:452892/1‑86 (MQ=255)
                                                                                    tGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACa                                                   <  1:149300/139‑1 (MQ=255)
                                                                                     ggCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcat                                                                      <  1:154824/119‑1 (MQ=255)
                                                                                     ggCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTtcat                                                                      >  2:154824/1‑119 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTaataataat                                                                                                                            >  2:624771/1‑60 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTaataataat                                                                                                                            <  1:624771/60‑1 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCAc                                             >  2:389722/1‑139 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCACACAATAACATTAATAATAatcgatc                                                                                                                       <  2:264093/65‑1 (MQ=255)
                                                                                          agCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCACACAATAACATTAATAATAatcgatc                                                                                                                       >  1:264093/1‑65 (MQ=255)
                                                                                                     ccGGGTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGt                                  <  2:68478/139‑1 (MQ=255)
                                                                                                         gTCAGACTCTCCTGAAGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACt                              >  2:269475/1‑139 (MQ=255)
                                                                                                                       aaGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGc                <  1:269475/139‑1 (MQ=255)
                                                                                                                       aaGCCGACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGc                <  1:389722/139‑1 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                       <  2:511488/63‑1 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGAcggc                                                                                       >  1:511488/1‑63 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                                >  2:295292/1‑118 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGa                                <  1:295292/118‑1 (MQ=255)
                                                                                                                            gACCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGt           >  1:329393/1‑139 (MQ=255)
                                                                                                                             aCCGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGtt          <  2:112002/139‑1 (MQ=255)
                                                                                                                              ccGCACAATAACATTAATAATAATCGATCGTATTTTTGATGGTTTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGttt         >  1:167531/1‑139 (MQ=255)
                                                                                                                                gCACAATATCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTATCGTGACTGTCATATTCACAGCTCTGTTTTa       <  2:167531/139‑1 (MQ=255)
                                                                                                                                     aTAACATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTa  <  1:166592/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                       
TTGAAGGTGAATTGATGCACGAGTTCTTTATTAATGGTCAATATCGTAATGCCATTCGCATGTGTATATTCCAGCATCAGTATCTGGCAGAGCACAAAACACCGGGTCAGACTCTCCTGAAGCCGACCGCACAATAGCATTAATAATAATCGATCGTATTTTTGATGGTGTAAACCCGTTCGACGGCGGGTTTTACTCCTTCATCAACAATGATTAGCTGACAGTCCACCGGATTAGCGTGACTGTCATATTCACAGCTCTGTTTTACATTA  >  NC_000913/1656608‑1656879

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: