Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG49.1% 67.7 / 14.6 53noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (4/23);  new base G (26/0);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 3.23e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925760‑926047
                                                                                                                                            |                                                                                                                                                     
gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAcccc                                                                                                                                           <  1:162884/151‑1 (MQ=255)
gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAAcc                                                                                                                                             <  1:50886/149‑1 (MQ=255)
    tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGa                                                                                                                                        <  1:220067/150‑1 (MQ=255)
        ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGt                                                                                                                                   <  2:809330/151‑1 (MQ=255)
          cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTg                                                                                                                                 <  1:882359/151‑1 (MQ=255)
              ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTGcc                                                                                                                               >  1:1017757/1‑149 (MQ=255)
               tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATATTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGt                                                                                                                             <  2:488159/150‑1 (MQ=255)
                    aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAACCCCCCATACGTTGCCGTATaca                                                                                                                        >  2:424262/1‑150 (MQ=255)
                       cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACt                                                                                                                    <  1:56101/151‑1 (MQ=255)
                                     tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGATCGCACCCCCGATACGTTTTCCTATACACACTTTTCAGGCgag                                                                                                         >  2:320386/1‑146 (MQ=21)
                                          ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCGTACGGTTCCGTATACACACTTTTCAGGCGGGCTCCtt                                                                                                   >  2:899036/1‑149 (MQ=18)
                                            gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTCGAAACCCCCATACGTTGTCGTATACAAACTTTCCCGGGGTGccccct                                                                                                   >  1:388527/1‑145 (MQ=18)
                                            gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCCGAACCCCCCTACGTTTCCGTATACACACTTTTCCGGCGGGCCCCCTCAGc                                                                                               >  1:566036/1‑151 (MQ=18)
                                              ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCa                                                                                             <  1:557711/151‑1 (MQ=255)
                                                accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACGTTTCCGTATACACACTTTCCCGGCGGGCTCCTTCAGcc                                                                                              >  2:430830/1‑148 (MQ=255)
                                                accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACt                                                                                           <  2:535457/151‑1 (MQ=255)
                                                 ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGGACCCCCGCTACGTTGCCGTATACACACTTTTCCGGCGTGCGCCTTTAGCCAccc                                                                                          >  1:842254/1‑149 (MQ=18)
                                                        cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGCACCCCCCCTACGTTGCCGTATACACACCTTTCAGGCGGGCTCCCTCCGCCACTCGGACAcc                                                                                   >  1:320360/1‑151 (MQ=17)
                                                           tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGAACCCCCAAACGTTTCCGTAAAAACACTTTTCAGGGGTGCTCCTTCCGCCCCCCGGACACCTc                                                                                 >  2:399004/1‑150 (MQ=17)
                                                            gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATACACACCTTTCAGGGGGGCTCCCTCCGCCAccccggacccc                                                                                  >  1:571995/1‑138 (MQ=17)
                                                             ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGGACCCCCCATACGTTGCCGTATACACACATTCCCGGCGGGCTCCTTCCGCCACCCCGACACCTCAcc                                                                              >  2:451990/1‑151 (MQ=17)
                                                              cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCGATAAGTTGCCGTATACACACCTTCCAGGCGTGCTCCTTCAGCCACCCGGACACCTCAcc                                                                              >  2:264773/1‑150 (MQ=17)
                                                                 acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATAAGTTGCCGTATACACACTTTTCCGGCGTGCTCCTTCAGCCACCCGGACACCCCACCaaa                                                                           >  2:139809/1‑150 (MQ=14)
                                                                  cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGGGGGGGGGGTTTGAAACCCCGCTAAGCTGCCCGATACACACCTTTCCGGCGCCCTCCTCTCGCCCCCCCGGCACCTCAcc                                                                              >  2:361/1‑146 (MQ=9)
                                                                    gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGt                                                                       <  2:113689/151‑1 (MQ=255)
                                                                     acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCAAACGTTTCCGTATACACACTTTTCCCGCGGGGTCCTTCAGCCCCCCGGGCACCCCCCCAAAATGtt                                                                      >  2:679732/1‑151 (MQ=14)
                                                                      cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTTCCGTATAAACACTTTCCCGGGGGGGTCCTTCAGCCACCCCGACACCCCCCCAAAATtttt                                                                     >  1:206846/1‑151 (MQ=14)
                                                                      cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTGCCGTATACACACCTTCCAGGCGGGGTCCCTCCGCCCCCCCGGCACCCCCCCCAAATGttt                                                                     >  2:623608/1‑151 (MQ=14)
                                                                      cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCCTACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCCGACACCCCCCCCAAATGttt                                                                     >  1:43493/1‑151 (MQ=14)
                                                                           aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCCTACGTTTCCGGATAAACACTTTCCCGGCGGGGCCCTTCAGCCACTCCGACAACCCACCAAAATGTTTTTCt                                                                 >  2:924742/1‑150 (MQ=14)
                                                                              aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCCAACCCCCCCTACCTTGCCGCCTACACACCTTCCCCGCGCGCTCCCTCCCCCCCCCCGCCACCCCCCCCAAATGTTTTTCTccc                                                              >  1:516239/1‑150 (MQ=2)
                                                                                  gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCATAAGTTGGCGTATACACAATTTTCAGGCGGGGTCCCTCAGCCACTCGGGCACCCCACCCAAATGTTTTGCTACCCAAcc                                                         >  1:342564/1‑151 (MQ=14)
                                                                                     aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCa                                                      <  2:731467/151‑1 (MQ=255)
                                                                                         aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCCTACGTTGCCGTATAAACACTTTCCAGGCGTGGTCCTTCAGCCACCCGGACAACCCACCAAAATGTTTTGCTACCCAACCCCATggg                                                  >  2:54250/1‑151 (MQ=14)
                                                                                           gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTg                                                <  2:794098/151‑1 (MQ=255)
                                                                                            cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTTTCGACCCCCCATACGTTGCCGGATATACACCTTTCAGGCGTGCCCCCTCAGCCACTCGGACCCCCCCCCAAATTGTTTTTCTACCCCACCACCTgggggg                                               >  1:220383/1‑151 (MQ=2)
                                                                                                cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAAc                                           <  2:953990/151‑1 (MQ=255)
                                                                                                 aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                                          <  2:923253/151‑1 (MQ=255)
                                                                                                     tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGCTGCCGGATACACACCTTCCAGGCGGGCTCCTTCAGCCACTCCGACACCCCACCAAATTTTTTTTCTACCAAAACTCATGGGGGGGAACgggg                                       >  2:150518/1‑150 (MQ=14)
                                                                                                      cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGgc                                      <  2:401248/150‑1 (MQ=255)
                                                                                                       gTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGGAACCCCCGTACCTTGCCGTATACACACTTTCCAGGGGTGCTCCTTCAGCCACCCGGGCACCCCACCAAAATGTTTTTCTACCAAAACTCCTGGGGGGCAAAGGGgcg                                     >  1:683448/1‑150 (MQ=14)
                                                                                                              ctctTCTTTAAATATGGCGGGGAGGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACACACCTTTCAAGCGCGCTCCCTCCGCCCCTCCGGCACCCCACCCAATTTTTTTTTTTCCCAACCCCATGGGGGGGAAAGGGGGGccacta                               >  1:546341/1‑149 (MQ=2)
                                                                                                                     ttAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAg                        <  1:551658/149‑1 (MQ=255)
                                                                                                                      tAAATAAGGCGGTGAGGGGGGGATTCGGACCCCCGATTCGTTGGCGTATACCCCCTTTCCAGGCGTTCCCCCTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                      <  2:718536/150‑1 (MQ=21)
                                                                                                                           atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCCAACCTCCTGGGTGGCAACGGGGCGCTACTATAGGGGGTTGGAg                  >  1:54867/1‑149 (MQ=255)
                                                                                                                           atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                <  2:854477/151‑1 (MQ=255)
                                                                                                                           atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGGTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGGGTTGGAg                  >  1:1027272/1‑149 (MQ=255)
                                                                                                                             ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa               <  1:11308/150‑1 (MQ=255)
                                                                                                                              gCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa             <  2:393673/151‑1 (MQ=255)
                                                                                                                                  tGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         >  1:257167/1‑151 (MQ=255)
                                                                                                                                   gAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         <  2:324673/150‑1 (MQ=255)
                                                                                                                                     gggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa      <  1:590624/151‑1 (MQ=255)
                                                                                                                                         gggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag  <  1:559228/151‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
GCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925760‑926047

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