Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 925,899 (G)7→8 noncoding (73/88 nt) serW ← tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,8921.G93.3% 102.7 / 0.9 30noncoding (80/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (11/17);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGG  >  NC_000913/925756‑926022
                                                                                                                                         |                                                                                                                                  
cACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgg                                                                                                                                  >  1:529614/1‑139 (MQ=255)
           aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCg                                                                                                                        <  1:463194/139‑1 (MQ=255)
                 cTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAAcccc                                                                                                                  >  2:440165/1‑139 (MQ=255)
                   tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTGCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCg                                                                                                                <  1:440165/139‑1 (MQ=255)
                        aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACg                                                                                                           <  2:672366/139‑1 (MQ=255)
                           cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGa                                                                                                                       <  2:381148/124‑1 (MQ=255)
                           cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGa                                                                                                                       >  1:381148/1‑124 (MQ=255)
                                    cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                               <  1:44476/139‑1 (MQ=255)
                                      tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATacaca                                                                                             <  1:600384/139‑1 (MQ=255)
                                             cccGCCCACCATGTTCACTGCCGTACAGCCAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAACTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGtt                                                                                                         >  2:433165/1‑120 (MQ=21)
                                             cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGtt                                                                                                         <  1:433165/120‑1 (MQ=255)
                                             cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAg                                                                                     >  1:542863/1‑139 (MQ=255)
                                                                          cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCt                                                         <  2:657982/139‑1 (MQ=255)
                                                                            gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACg                                                                                                           <  1:491462/87‑1 (MQ=255)
                                                                            gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACg                                                                                                           >  2:491462/1‑87 (MQ=255)
                                                                             aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCAc                                                      >  1:93836/1‑139 (MQ=255)
                                                                                       aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGtt                                                                                                         <  2:22511/78‑1 (MQ=255)
                                                                                       aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGtt                                                                                                         >  1:22511/1‑78 (MQ=255)
                                                                                             aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTAcc                                     <  2:542863/139‑1 (MQ=255)
                                                                                                      cTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                             <  2:93836/139‑1 (MQ=255)
                                                                                                          cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGTACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATggg                         >  1:494783/1‑139 (MQ=255)
                                                                                                                aGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACTCTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCaa                   <  2:467875/139‑1 (MQ=255)
                                                                                                                  ctctTCTTTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                 <  2:494783/139‑1 (MQ=255)
                                                                                                                       cttTAAATATGGCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg            <  2:529614/139‑1 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCa                                    <  2:428199/105‑1 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCa                                    >  1:428199/1‑105 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg            >  2:35511/1‑129 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg            <  1:679106/129‑1 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg            >  2:679106/1‑129 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgcg            <  1:35511/129‑1 (MQ=255)
                                                                                                                                 ggCGGTGAGGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAgg  >  1:688673/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                  
CACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGG  >  NC_000913/925756‑926022

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: