Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG42.9% 114.5 / 24.5 77noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (3/41);  new base G (33/0);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.14e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAG  >  NC_000913/925755‑926044
                                                                                                                                                 |                                                                                                                                                   
tCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGtttg                                                                                                                                                 >  1:1740/1‑148 (MQ=255)
       cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCg                                                                                                                                         <  2:770583/151‑1 (MQ=255)
        gTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAacccc                                                                                                                                          >  2:883865/1‑149 (MQ=255)
         tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTGCGataacccc                                                                                                                                         >  2:932018/1‑142 (MQ=38)
          aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATa                                                                                                                                      <  1:1111471/151‑1 (MQ=255)
           cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATAc                                                                                                                                     <  1:330242/151‑1 (MQ=255)
           cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATAc                                                                                                                                     <  2:574248/151‑1 (MQ=255)
            aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGGCGAACCCCCGATa                                                                                                                                      >  1:920232/1‑149 (MQ=255)
            aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACg                                                                                                                                    <  1:125462/151‑1 (MQ=255)
             ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACg                                                                                                                                    <  1:895860/150‑1 (MQ=255)
              gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGtt                                                                                                                                  <  1:14287/151‑1 (MQ=255)
                aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGTGGGGGGGGTTTGGAACCCCCCTACGtt                                                                                                                                  >  2:1109082/1‑149 (MQ=37)
                 gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGTTTGAAACCCCCATACGTTGcc                                                                                                                               >  1:1130425/1‑151 (MQ=255)
                  cTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTATAAcccccccgaccgtggc                                                                                                                               >  1:718194/1‑139 (MQ=37)
                   ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCGATACGTTGcc                                                                                                                               >  1:191228/1‑149 (MQ=255)
                    tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTCGGACCCCCGCTACGTTGCCGt                                                                                                                             >  2:862236/1‑150 (MQ=255)
                      gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGtata                                                                                                                          <  2:569743/151‑1 (MQ=255)
                       aaaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCATACGTTGCCGtat                                                                                                                           >  1:564120/1‑149 (MQ=255)
                        aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATaca                                                                                                                        <  2:1160522/151‑1 (MQ=255)
                         aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCCTAAGTTGCCGcat                                                                                                                           >  2:743425/1‑147 (MQ=25)
                            cacaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATacacac                                                                                                                    <  1:749519/150‑1 (MQ=255)
                             acaGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGGACCCCCCACACGGTGCCGTATaaacac                                                                                                                    >  2:766928/1‑149 (MQ=25)
                              caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACtt                                                                                                                  <  2:565264/150‑1 (MQ=255)
                                  tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAACCCCCCATACGTTGCCGTATAAACACCTTcc                                                                                                               >  1:1029476/1‑149 (MQ=255)
                                    aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCGATACGTTGCCGTATACACACTTTTCacg                                                                                                            >  2:791764/1‑148 (MQ=255)
                                    aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGc                                                                                                           <  2:684310/151‑1 (MQ=255)
                                    aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGc                                                                                                           <  1:261999/151‑1 (MQ=255)
                                     cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGAAACCCCCCTACGGTTCCGTATACACAATTTCCCgg                                                                                                            >  1:663042/1‑149 (MQ=21)
                                     cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATAAGTTGCCGTATACACACTTTCCAGGCg                                                                                                          <  1:628383/151‑1 (MQ=255)
                                        aCTCCACCCGCCCACCATATTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTCGGAACCCCCCTACGTTGCCCTATACACACCTTcc                                                                                                               >  2:686206/1‑143 (MQ=21)
                                         cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGc                                                                                                       <  1:1202951/150‑1 (MQ=255)
                                          tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAAGGGGGGGTTCGCAACCCCGGTAAGTTGTCGTGTACAAACATTTCATGCGTGCt                                                                                                      >  2:350530/1‑150 (MQ=21)
                                               ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCa                                                                                               <  1:921461/151‑1 (MQ=255)
                                                   ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTCGAACCCCCCATACGTTGCCGTATACACACCTTCCCGGCGTGCTCCTTCAGcc                                                                                            >  2:37372/1‑150 (MQ=18)
                                                    caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGGTCGGAACCCCCATAAGTTGCCGCATACACAACTTTCCAGCGGGGTCCTTCAGcc                                                                                            >  1:303277/1‑149 (MQ=18)
                                                     accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGACCCCCGATACGGTTCCGTATAAACACTTTTCAGGCGGGCTCCTTTAGcccc                                                                                          >  1:734746/1‑148 (MQ=18)
                                                         tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCGAAACGTTGCCGTATACACACCTTCCCGGCGGGGTCCTTCAGCCACCCgg                                                                                      >  1:933181/1‑150 (MQ=18)
                                                            tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                   <  2:170386/150‑1 (MQ=255)
                                                             cACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGacac                                                                                  <  2:161878/150‑1 (MQ=255)
                                                               cTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGTTACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCt                                                                                <  2:955360/150‑1 (MQ=255)
                                                                     tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCGATACGTTGCCGTATACACACTTTTCCGGGGTGCTCCTTTCGCCACCCGGGCACCCCAccca                                                                          >  1:891684/1‑148 (MQ=17)
                                                                         gacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCCCCCCCAATTGt                                                                    >  1:739096/1‑151 (MQ=255)
                                                                          acagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGtt                                                                   <  2:156320/151‑1 (MQ=255)
                                                                           cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCCATACGTTGCCGTATACACACTTTTCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGttt                                                                  >  1:674255/1‑151 (MQ=255)
                                                                             gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTGCCGTATACACACTTTCCAGGCGGGCTCCTTCAGCCACCCGGGCACCCCCCCAAAATGtttt                                                                 >  1:1103320/1‑150 (MQ=14)
                                                                                aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAACCCCCGCTACGTTGCCGGATACACACTTTTCAGGCGGGGTCCTTTAGCCACTCGGGCACCCCACCAAAATTTTTTGCt                                                              >  1:734920/1‑150 (MQ=14)
                                                                                     gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAAACCCCCCAACGTTGCCCTATAAAAACTTTTCCGGGGGGCTGCCTTTGCCACACGGACAACCCCCCAAATAGTTTTGCTAACaa                                                         >  2:389449/1‑150 (MQ=9)
                                                                                        aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAAcc                                                      <  2:821020/150‑1 (MQ=255)
                                                                                          aaaaaaGCTCGCACTTGCGTACGGGCTCTTCTTTAAAAATGGCCGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCCCCAAATTGTTTTGCTACCAAACCt                                                     <  2:1011178/149‑1 (MQ=255)
                                                                                           aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCAt                                                  <  2:548786/151‑1 (MQ=255)
                                                                                              aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCTGACCCCCCCCACCGTTCCGGTTACACCCTTTCCAGGCGTGCCCCTTCCCCCCCCCGGCCACCCCCCCCAAATGTTTTGCTTCCaaaaccccttgg                                               >  2:291528/1‑142 (MQ=9)
                                                                                              aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCGAAACCCCCCTACCTTTCCGTATAAACACTTTTCAGGCGGGCCCCTTCAGCCACCCCGACCCCCCCCCAAAATGTTTTTCTTCCAAACCTCCTggg                                               >  2:632970/1‑151 (MQ=14)
                                                                                              aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATggg                                               <  2:182059/151‑1 (MQ=255)
                                                                                              aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTGAACCCCCCCTACGTTTGCGTTTAAACACCTTCCACGGGGGGTCCCTCCGCCACTCGGGAACCCCACCCAAATGTTTTTCTTCCAAACCCCCTggg                                               >  2:253366/1‑151 (MQ=9)
                                                                                               aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCCCTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGt                                              <  2:1115959/151‑1 (MQ=255)
                                                                                               aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCCATACGGTTCCGTATACACACTTTCCAGGCGTGCTCCCTCAGCCACTCGGACACCTCACCAAAATTTTTTTCTACCAAACCTCCTggg                                               >  1:101168/1‑150 (MQ=14)
                                                                                                   cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTgg                                            <  2:494550/149‑1 (MQ=255)
                                                                                                    gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTCCAACCCCCCCTACGGTGCCCTTTACACACCTTTCAGGGGGGCCCCTTCAGCCACCCCGGCACCTCACCAAAATGTTTTTTTCCCaaaaccccttggggg                                            >  1:893898/1‑136 (MQ=2)
                                                                                                     cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGAACCCCCCATACGTTGCCGTATACACACTTTTCAGGCGGGCTCCCTCAGCCACCCGGGCACCCCACCAAAATTTTTTTCTACCCAACCCCATGGGGGGCaa                                         >  1:356188/1‑150 (MQ=14)
                                                                                                         ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgg                                      <  2:960964/149‑1 (MQ=255)
                                                                                                          tCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCCATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCCCCAAATTGTTTTTCTACCAAACCTCATGGGTGGCAACgggg                                    >  1:721812/1‑150 (MQ=255)
                                                                                                               cGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctac                              <  1:542722/151‑1 (MQ=255)
                                                                                                                 aGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGctacta                            <  1:854435/151‑1 (MQ=255)
                                                                                                                 aGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCGACCAAACCTCATGGGTGGCAACGGGGCGctact                             <  2:945562/150‑1 (MQ=255)
                                                                                                                  gCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTCGGACCCCCGCTTAGTTTTCCTATATACAATTTTCCGGCGGGCGCCTTCCACCCCCCCGGACCCCCACCAAAATTGTTTTTTTCCCAACCCCAATgggggtcacggg                                      >  2:899068/1‑132 (MQ=2)
                                                                                                                       tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAgg                        <  1:374901/149‑1 (MQ=255)
                                                                                                                       tcttTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                         <  2:652266/148‑1 (MQ=255)
                                                                                                                           tAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                   <  2:260497/150‑1 (MQ=255)
                                                                                                                                atGGCGGTGAGGGGGGGATTCGCACCCCCGATACGTTGCCGTATAAACACTTTCCAGGCGTGCTCCTTCAGCCCCTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGGGGCAACGGGGCGCTACTATAGGGAGTTGGAGt              >  2:918112/1‑150 (MQ=255)
                                                                                                                                atGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa             <  1:883888/151‑1 (MQ=255)
                                                                                                                                  ggCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa           <  1:758010/151‑1 (MQ=255)
                                                                                                                                   gCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaaa          <  1:1053666/151‑1 (MQ=255)
                                                                                                                                     ggTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg        <  1:213864/151‑1 (MQ=255)
                                                                                                                                     ggTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg        <  2:1130395/151‑1 (MQ=255)
                                                                                                                                         aGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa    <  2:573323/151‑1 (MQ=255)
                                                                                                                                         aGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa    <  2:279129/151‑1 (MQ=255)
                                                                                                                                           ggggggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaag  <  2:356180/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAG  >  NC_000913/925755‑926044

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