Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913925,9000AG41.8% 119.9 / 19.8 79noncoding (72/88 nt)serWtRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (7/39);  new base G (33/0);  total (40/39)
Fisher's exact test for biased strand distribution p-value = 1.17e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925757‑926050
                                                                                                                                                |                                                                                                                                                       
aCTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCgga                                                                                                                                                   >  1:722058/1‑148 (MQ=255)
  tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGT‑GGGGTGGGttgcgctac                                                                                                                                                >  2:1165892/1‑144 (MQ=255)
   gCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAccc                                                                                                                                                <  2:893224/150‑1 (MQ=255)
    ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAccccc                                                                                                                                              <  2:453782/151‑1 (MQ=255)
        aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                           <  2:43526/150‑1 (MQ=255)
        aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                           <  2:1115946/150‑1 (MQ=255)
           ggCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑t                                                                                                                                      <  2:315095/151‑1 (MQ=255)
            gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCCTACG‑tt                                                                                                                                     >  1:1040922/1‑151 (MQ=255)
              aGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGA‑GGGGGGGGGTCACAccccccgtccgtt                                                                                                                                      >  1:449163/1‑141 (MQ=37)
                 ttAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAGTTGAACCCCCCATACG‑TTGTCg                                                                                                                                 >  2:375178/1‑150 (MQ=255)
                  tAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCGACACggttac                                                                                                                                   >  2:682842/1‑143 (MQ=255)
                       aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGATCGAACCCCCCATACG‑TTGCCGGATaca                                                                                                                           >  1:829139/1‑150 (MQ=255)
                        ttCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACG‑TTTCCGGATacac                                                                                                                          >  1:132676/1‑150 (MQ=255)
                         tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACG‑TTGGCCGATacac                                                                                                                          >  1:1274444/1‑149 (MQ=255)
                            caGGTAACGTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACttt                                                                                                                     <  1:93505/151‑1 (MQ=255)
                            caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCAACG‑TTTCCCTATAAACACttt                                                                                                                     >  1:241044/1‑151 (MQ=25)
                              ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCCTACGTTTGCCG‑ATACACACctt                                                                                                                     >  1:665190/1‑149 (MQ=25)
                                 aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAg                                                                                                                 <  1:790808/150‑1 (MQ=255)
                                       cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTGTGAACCCCCCATACG‑GTGCCGGATAAACACCTTcccagcggggt                                                                                                          >  1:587362/1‑142 (MQ=21)
                                         ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCt                                                                                                          <  1:498694/149‑1 (MQ=255)
                                         ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTc                                                                                                         <  2:1025024/150‑1 (MQ=255)
                                         ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTc                                                                                                         <  1:329557/150‑1 (MQ=255)
                                         ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTc                                                                                                         <  1:1199068/150‑1 (MQ=255)
                                             ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTCGAACCCCCCATACG‑TTGCCCCATACCCACCCTCCAcgcggg                                                                                                            >  1:1676/1‑141 (MQ=21)
                                              cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGAACCCCCATACG‑TTGCCGTATACACACTTTTCAGGCGGGCTCCtt                                                                                                      >  1:944785/1‑148 (MQ=18)
                                              cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGACCCCCCCTACG‑GTTCCCTATACACACCTTCCCGGCGTGCTccct                                                                                                      >  1:1163243/1‑146 (MQ=18)
                                              cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCACACG‑TTGCCGTATACACACTTTCCAGGCGGGCTccct                                                                                                      >  1:552754/1‑146 (MQ=255)
                                                  caccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                               <  1:376690/151‑1 (MQ=255)
                                                   accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCCGACCCCCGCAACG‑GTTCCGGATACACACCTTCCCGGCGGGGTCCTTcttc                                                                                                  >  1:234381/1‑144 (MQ=18)
                                                   accaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACt                                                                                              <  2:898951/151‑1 (MQ=255)
                                                    ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTc                                                                                             <  1:11960/151‑1 (MQ=255)
                                                    ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑AGGGGGGGGTTGGACCCCCGACACG‑GTTTCGCATAAAACCTTTCCCAGCGGGGTCCTTCAGccccac                                                                                             >  2:484057/1‑147 (MQ=16)
                                                     caTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCg                                                                                            <  2:1128579/151‑1 (MQ=255)
                                                      aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCgg                                                                                           <  2:716760/151‑1 (MQ=255)
                                                      aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGTTTGTACCCCCCATAAG‑TTGCCGCATAAACACTTTCCCGGCGGGGGCCTTCAGCCACTc                                                                                             >  2:1313962/1‑149 (MQ=18)
                                                              tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACG‑TTTTCGCATACACACTTTCCAGGCGGGGCCCCTCCGCCACCCCGACAcccc                                                                                    >  2:1171624/1‑148 (MQ=17)
                                                              tGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGAACCCCCCCTACG‑TTGGCGTATACCGAATTTCCCGGCGGGCGCCCTCAGCCCCTCAGCCAcccc                                                                                    >  2:916389/1‑148 (MQ=17)
                                                               gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCAc                                                                                  <  1:1275567/151‑1 (MQ=255)
                                                                 cGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGAAACCCCGAAAAG‑GTGCCCGATACACACTTTTTAGGCGGGCCCCCTCAGCCACCCGGgcaca                                                                                      >  1:404791/1‑144 (MQ=17)
                                                                   tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTAA‑GGGGGGGGTTTCAACCCCCCAGACG‑TTGGCGTATACACACTTTTCCGGGGTGCTCCTTCCGCCCCTCCGTCACCCCAccca                                                                               >  2:232644/1‑148 (MQ=17)
                                                                     cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTAA‑GGGGGGGAGTTTTACCCCCGATACG‑ATGCCGTTTAAAAAATTTTCAGGCGGGGTCCTTCAGCCCCTCCGCAACCTccccca                                                                               >  2:57894/1‑146 (MQ=14)
                                                                          agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCACAAG‑TTGTCGTATACAAACTTTTCCGGCGGGCTCCCTTACCCCCTCCGCAACCTCCCCCAAAtttttt                                                                       >  2:1138986/1‑151 (MQ=9)
                                                                          agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCGCTACG‑GTGCCGTATACAAACTTTTCAGGCGGGGTCCCTCAGCCACCCGGACACCCCACCCAAAtttttt                                                                       >  1:145221/1‑151 (MQ=14)
                                                                             tAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTGGACCCCCCATACG‑GTGTCGCATACACACCTTCCAGGGGTGCGCCCTCCCCCACTCGGACACCCCACCaaaaatttttttt                                                                    >  2:43339/1‑141 (MQ=14)
                                                                                 aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTA‑GGGGGGGGTTTGAACCCCCCACAAG‑GTGCCGTATACACACTTTTCAAGCGGGCTCCTTCAGCCACCCCGGCACCCCCCCCAAATTTTTTTCTAcc                                                                 >  1:174941/1‑150 (MQ=14)
                                                                                 aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGATTGAAACCCCCCTACG‑GTTTCGTAATCACACTTTTCCGGGGCGCTCCTTCACCCACTCGGGCAGCCCCCCaacatggtttcttaac                                                                 >  2:909020/1‑136 (MQ=2)
                                                                                   gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGTCGGTGA‑GGGGGGGATTCTAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCTCCAAATTGTTTTGCTACCaa                                                               <  2:864124/150‑1 (MQ=255)
                                                                                      aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGGTTGGAACCCCC‑ACAAG‑TTGTCCTATACACACATTTCAGGCGGGGTCCTTCATCCACCCGGGCACCCCCCCAAATTTTTTTTCTTCCCAAc                                                             >  2:2299/1‑149 (MQ=2)
                                                                                       aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTc                                                          <  2:1100236/151‑1 (MQ=255)
                                                                                         aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAGCCCCCGATACG‑TTACCTTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACTCCTCACCAAATTGTTTTTCTACCAAGCCTc                                                          <  2:209780/149‑1 (MQ=14)
                                                                                         aaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCAt                                                        <  1:407871/151‑1 (MQ=255)
                                                                                             aGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCGCCAAATTGTTTTGCTACCAAACCTCATGGGt                                                    <  1:206679/151‑1 (MQ=255)
                                                                                                 cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGc                                                 <  1:113854/150‑1 (MQ=255)
                                                                                                  gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCGTACG‑TTTGCGTATAAACACTTTTCCGGCGTGCTCCTTCCGCCACCCCGACACCCCACCCAATTGTTTTTCTTCCAAAACCCCTgggggg                                                  >  2:941767/1‑148 (MQ=2)
                                                                                                  gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCATACG‑TTGCCGTATACACACTTTTCCGGCGTGCTCCTTCAGCCCCCCGGACACCCCCCCCAAATGTTTTTCTACCAAACCCCATggggggg                                                 >  1:570259/1‑148 (MQ=14)
                                                                                                    aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCCAATG‑TTTCCGGATACACACCTTCCCGGCGTGCCTCTTCAGCCACTCGGGCCCCCCACCCAAATTTTTTTTTACCCAACCCCATCGGGGGCaa                                               >  2:885443/1‑149 (MQ=2)
                                                                                                      tttCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCCAACG‑TTGCCGTATACACACTTTCCAGGCGGGCTCCTTCAGCCCCCCCGGCACCCCACCCAAATGTTTTTCTACCAAAACTCCTGGGGGGGaaacgg                                           >  2:1053943/1‑147 (MQ=9)
                                                                                                        tCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCCCACG‑TTGTCGTATACACACTTTTCCGGCGGGGTCCTTCCGCCACCCGGACACCCCACCAAAATGTTTTGCTTCCAAAACCCCTGGGGGGGAACgggg                                          >  2:596618/1‑150 (MQ=9)
                                                                                                         cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACgggg                                          <  2:703360/149‑1 (MQ=255)
                                                                                                            aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGcta                                     <  1:1137923/151‑1 (MQ=255)
                                                                                                            aCGAGCTCTTCTTTAAATATGGCGGAGG‑GGGGGGGAGTCGGACCCCCCATAAG‑TTTTCGTATACACACTTTCCAGGGGGGCTCCTTCAGCCCCCCCGCCACCACACCCAAATGTTTTTTTACCAAACCCCCTGGGGGGcgaaaggggggtc                                     >  2:966216/1‑139 (MQ=2)
                                                                                                                   cttcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACtata                                <  1:1269965/149‑1 (MQ=255)
                                                                                                                     tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTGCCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                               <  1:842336/148‑1 (MQ=255)
                                                                                                                      cttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCCATACG‑GTGGCGTTTACACACTTTTCAGGCGTGCTCCTTCAGCCACCCGGACACCCCACCCAAATTTTTTTCTTACAAAACTCCTGGGGGGCAAAAGGGGGCCACtattagggg                           >  1:405512/1‑145 (MQ=2)
                                                                                                                         tAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  1:231573/150‑1 (MQ=255)
                                                                                                                         tAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                         <  1:566596/150‑1 (MQ=255)
                                                                                                                          aaaTATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                       <  2:470877/151‑1 (MQ=255)
                                                                                                                           aaTATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa                      <  1:1120214/151‑1 (MQ=255)
                                                                                                                              atGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                   <  1:260704/151‑1 (MQ=255)
                                                                                                                                ggCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa                 <  2:722392/151‑1 (MQ=255)
                                                                                                                                ggCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                  <  1:682546/150‑1 (MQ=255)
                                                                                                                                      gA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           <  1:980616/151‑1 (MQ=255)
                                                                                                                                      gA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc           >  2:1333212/1‑151 (MQ=255)
                                                                                                                                       a‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          <  1:1165225/151‑1 (MQ=255)
                                                                                                                                         gggggggATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa         <  1:957921/151‑1 (MQ=255)
                                                                                                                                           gggggATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa          >  1:1204429/1‑148 (MQ=255)
                                                                                                                                             gggATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga       >  2:336441/1‑149 (MQ=255)
                                                                                                                                               gATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa   <  1:1188946/151‑1 (MQ=255)
                                                                                                                                               gATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa   <  2:263775/151‑1 (MQ=255)
                                                                                                                                                aTTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:357600/151‑1 (MQ=255)
                                                                                                                                                aTTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAt  <  1:134146/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                       
ACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACG‑TTGCCGTATACACACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAAT  >  NC_000913/925757‑926050

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: ‑