Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,511,9230AC44.9% 111.1 / 13.9 78V141G (GTT→GGT) yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (31/12);  new base C (0/35);  total (31/47)
Fisher's exact test for biased strand distribution p-value = 3.16e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGC  >  NC_000913/3511784‑3512069
                                                                                                                                                   |                                                                                                                                                  
cAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATg                                                                                                                                         >  2:837994/1‑151 (MQ=255)
   ccccACGCGCTTCACGTTGGGGGACTCGCCTTTACGGAAAACCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAC‑CTCCAGCGGGAGGAATAAG‑CCTTTTGGGGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGGGAACCCTCCGCCATGTCa                                                                                                                                      <  2:761744/151‑1 (MQ=255)
             ttCACTTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCACCGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGGGGCATAAG‑CTTTTTTGCTT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGcc                                                                                                                              <  2:2092/149‑1 (MQ=255)
               cACGTTTGCGACCCCGCTTTTACGCCAAATAAACGGGTG‑CCCGCCAACGAGGGTGGGGGCCCCCAC‑CTCCAGCGGGGGACATAAG‑CTTTTTGGCGG‑G‑AGATAAC‑TTC‑CC‑GCCGCAAGTCGAAAACCCCCCGCCATGTCATACACGCCAtt                                                                                                                           <  1:14967/150‑1 (MQ=255)
                   tttGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGc                                                                                                                       >  2:364952/1‑150 (MQ=255)
                     tggggcctcccttttCCGGAAAACCAACTGGTG‑CCCGCAACCTTAGGTGGCCGCCCCCAC‑CCCAAGCGCGGGACATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTC‑CC‑GCCGTTAGTCGCAAAACCACCGCCATGTCATACACGCCATTAAGCa                                                                                                                      <  2:573182/137‑1 (MQ=255)
                        ggccTCGCCTTTACGGCAAATCAACTGGTG‑CCCGCAAACGTAGGTGTGCGCCCCCAA‑CCCCAGCGCGAGGCATAAG‑CTTTTTTGCGT‑G‑AGATAAT‑TTC‑CC‑GCCGCTAGTCGCGAAACCCCCGCCATGTCATACACGCCATTAAGCAGaca                                                                                                                  <  1:738454/147‑1 (MQ=255)
                         gACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGacac                                                                                                                 >  2:496092/1‑150 (MQ=255)
                           cccGCCTTT‑CGGGAAATCCAATGGGGCCTCCCCACCGTAGGTG‑GGGGCCCCAACCCCCAGCGGGGGGCATAGG‑CTTTTTTGCGG‑G‑AAATACC‑TTT‑CC‑GCCGCTAGCCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTc                                                                                                               <  2:409439/148‑1 (MQ=255)
                            tCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTcc                                                                                                              >  2:160999/1‑150 (MQ=255)
                             cGCTTTTCCGGAAAATCAACGGGTG‑CTCCCCAACGTAGGTGTGCGCCCCCAA‑CCCCAGCGCGAGGCATAAG‑CTTTTTTGCGT‑G‑AGATAAT‑TTC‑CC‑GCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCAc                                                                                                            <  1:817227/151‑1 (MQ=255)
                             cGCTTTTACGGAAAATCAACTGGTG‑CTCGCCACCGTAGGTGTGCGGCCCAAA‑CCCCAGCGCGGGGCAAAAG‑CTTTTTTGCGT‑G‑AGATACC‑TTT‑CC‑GCCGCTAGTCGCAAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTcc                                                                                                              <  1:874134/149‑1 (MQ=255)
                             cGCTTTTACGGAAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCAAA‑CCCCAGCGCGGGGCATAAG‑CTTTTTTGCGT‑G‑AGATACC‑TTT‑CC‑GCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTcc                                                                                                              <  1:122691/149‑1 (MQ=255)
                              gCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCAcg                                                                                                           >  2:59741/1‑151 (MQ=255)
                                  ttACGGAAAATCAACGGGTG‑CCCCCCACCGTAGGTGTGCGGCCCCAA‑CCCCAGCGCGAGGCATAAG‑CTTTTTGGCGT‑G‑AGATAAC‑TTC‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGa                                                                                                        <  1:902015/150‑1 (MQ=255)
                                         aattCACCTGGTC‑CTCGCCAACGTGGTTGTGCGCCCCCAC‑CCCCACCGCGAGGCATAAG‑CTTTTTTGGGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTACTCGCAAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCcgc                                                                                                <  1:997006/148‑1 (MQ=255)
                                         aaaTCACCTGGTG‑CTCGCCACCGTGGGTGTGCGCCCCCAC‑CCCCAGCGCGGGGCATAAG‑CTTTTTTGCGT‑G‑AGATACC‑TTC‑CC‑GCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCcg                                                                                                 <  1:1078801/150‑1 (MQ=255)
                                         aaaTCAACTGGTG‑CTCGCCAACGTAGGTGTGCGCCCCCAA‑CCCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATcc                                                                                                  <  1:142948/149‑1 (MQ=255)
                                              aaCTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCAc                                                                                           >  1:670114/1‑151 (MQ=255)
                                              aaCTGGGG‑CTCGCCACCGTAGTTGTGCGGCCCCAC‑CCCCAGCGCGGGGCATAAC‑CTTTTTTGGGT‑G‑AAATAAT‑TTT‑CC‑GCCCCTATCCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCa                                                                                            <  2:321834/150‑1 (MQ=255)
                                                   gggccTCGCAACCGTGGGTGTGCGCCCCAAC‑CCCCGGCGCGAGGAAAAAG‑CTTTTTTGGGG‑A‑AGAAAAC‑TTC‑CC‑GCCGCTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAt                                                                                       <  2:857106/147‑1 (MQ=255)
                                                     g‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgtg                                                                                    >  1:760170/1‑151 (MQ=255)
                                                       ccggCCAAAGTAGGTGTGCGCCCCCAC‑CCCAAGCGAGGGCAATAGG‑CATTTTCGGGA‑G‑AGATAAC‑TTG‑CC‑CCGGATAGTCCCGAACCCACCTCCATGTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtg                                                                                      <  2:406595/145‑1 (MQ=255)
                                                        tCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTgcg                                                                                  >  1:931432/1‑151 (MQ=255)
                                                        tCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTgcg                                                                                  >  1:788014/1‑151 (MQ=255)
                                                        cccccAACGTGGGTGTGCGCCCCAAC‑CCCCGGCGCGGGGAATAAG‑CTTTTTGGCTT‑G‑AAATATC‑TTC‑CC‑GCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgt                                                                                     <  2:1004407/145‑1 (MQ=255)
                                                         caccatcttgggTGTGCGCCCCCAA‑CCCCAGCCCGGGGC‑TAAGCCCTTTTTGGGT‑G‑AGATACC‑TTC‑CC‑CCCGCTTGTCGCAAAACCACCGCCATTTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGc                                                                                         <  2:289298/133‑1 (MQ=255)
                                                           ccAAGGTGGGTGGGGGCCCCAAA‑CCCCAGCGGGGGGAATAAC‑CCTTTTGGCGG‑G‑AAATAAC‑TTC‑CC‑GCCGCTAGTCGCAAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTgc                                                                                   <  1:458667/147‑1 (MQ=255)
                                                            cagcGTAGGGGGGC‑GCCCCAA‑CCCCCGGGGGGGGGAAAAA‑CCTTTTTGGGG‑G‑GAAAAACGTTC‑CC‑GCCGCTAGCGGCAAACCCTCCGCCATGTAATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgt                                                                                     <  2:859340/141‑1 (MQ=255)
                                                            cAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                               >  1:758358/1‑150 (MQ=255)
                                                             aaCGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTC‑CC‑GCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATa                                                                              <  1:248080/150‑1 (MQ=255)
                                                                gTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTa                                                                           >  2:415726/1‑150 (MQ=255)
                                                                   gttgtgGGCCCCAAA‑CTCCAGCGAGGGAAAAAAG‑CTTTTTTGCGT‑G‑AGATACC‑TTC‑CC‑GCCACTATTCGCAAACCCCCCGCAATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaa                                                                        <  1:1103979/148‑1 (MQ=255)
                                                                   ggTGTGGGCCCCCAC‑CTCCAGCGGGGGGCATAAC‑CTTTTTGGGGT‑G‑AGATACC‑TTC‑CC‑GCCGCTAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaaa                                                                       <  2:690020/151‑1 (MQ=255)
                                                                   ggTGGGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CTTTTTTGCGT‑G‑AGATAAC‑TTC‑CC‑GCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaa                                                                         <  1:559658/149‑1 (MQ=255)
                                                                    gtggGCGGCCCCAA‑CCCCAGCGCGGGGAATAAC‑CTTTTTTGCGT‑G‑AAAAAAC‑TTT‑CC‑GCCGCTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGt                                                                            <  2:276251/145‑1 (MQ=255)
                                                                      gtgCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAgg                                                                     >  1:280485/1‑150 (MQ=255)
                                                                         cGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATg                                                                 >  1:248419/1‑151 (MQ=255)
                                                                            ccccAC‑CCCCAGGGGGGGGCAAAAC‑CTTTTTTGGGT‑G‑AAATACT‑TTC‑CC‑CCCGCTAGCCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGa                                                                    <  2:251898/145‑1 (MQ=255)
                                                                            ccccAC‑CCCCAGCGCGGGGCAAAAC‑CTTTTTGGCGT‑G‑AAATACC‑TTC‑CC‑GCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAAt                                                                  <  1:67867/147‑1 (MQ=255)
                                                                               ccacccccAGCGCGAGGCAAAAC‑CTTTTTTGCGT‑A‑AGAAAAC‑TTC‑CC‑GCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATc                                                            <  1:335518/145‑1 (MQ=255)
                                                                                aa‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCa                                                           >  1:202906/1‑150 (MQ=255)
                                                                                 acccccAGCGGGGGGAAAAAC‑CTTTTTTGCGT‑G‑AAATACC‑TTC‑CC‑GCCGCTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCa                                                           <  1:752268/146‑1 (MQ=255)
                                                                                   ccccAGCGCGAGGAAAAAG‑CTTTTTTGCGG‑G‑AAATAAC‑TTC‑CC‑GCCGCTAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCaaa                                                         <  1:218344/148‑1 (MQ=255)
                                                                                    tCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACt                                                       >  2:697343/1‑151 (MQ=255)
                                                                                         cgcgAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCt                                                  >  2:445509/1‑151 (MQ=255)
                                                                                           gggggCAAAAG‑CCTTTTTGCGG‑GAAGATACT‑TTC‑CC‑CCCGCTAGTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTggg                                                    <  2:766982/145‑1 (MQ=255)
                                                                                           gggaggaataG‑CCTTTTTGGGTGG‑AAAAAAT‑TTC‑CC‑GCCGCTGGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGc                                                   <  2:196479/140‑1 (MQ=255)
                                                                                            gAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTc                                               >  2:701358/1‑151 (MQ=255)
                                                                                            gAGGAATAAG‑CTTTTTTGCGG‑G‑AGATACC‑TTC‑CC‑GCCGCTAGTCGCAAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTggg                                                    <  1:304805/146‑1 (MQ=255)
                                                                                                 aTAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATcc                                          >  1:935662/1‑151 (MQ=255)
                                                                                                   taacctttttgGCGT‑G‑AGATACT‑TTT‑CC‑CCCGCTGGGCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCt                                                  <  1:966981/136‑1 (MQ=255)
                                                                                                       ccTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTc                                      >  2:322450/1‑150 (MQ=255)
                                                                                                         tttttGGGT‑G‑AGATACT‑TTC‑CC‑GCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTc                                      <  2:771126/148‑1 (MQ=255)
                                                                                                         tttttGCGT‑G‑AGA‑AAC‑TTTCCC‑GCCGCTGGCCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTc                                      <  1:476693/148‑1 (MQ=255)
                                                                                                          ttttGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGcc                                   >  2:1079256/1‑150 (MQ=255)
                                                                                                           tttGGGA‑G‑AGATACT‑TTC‑CC‑GCCCTTAGCGCCGAACCCACCGCCATGTCGTACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCa                                  <  2:1074143/150‑1 (MQ=255)
                                                                                                           tttGCGT‑G‑AGAAACC‑TTC‑CC‑GCCGTTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCa                                  <  1:833561/150‑1 (MQ=255)
                                                                                                            ttttggtgtggaTAAC‑CTT‑CCGGCGGTTAGTCGCGAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATTCCATCAAACTGGGCTTTCAATCCGCTCg                                     <  2:1073419/139‑1 (MQ=255)
                                                                                                            ttGCGT‑G‑AAATACC‑TTC‑CC‑GCCGCTGGCCGCGAACCCCCCGCCATGTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAg                                 <  2:710254/150‑1 (MQ=255)
                                                                                                             ttgggtggaaaaacc‑TTC‑CC‑GCCGCTGGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCa                                  <  2:1104150/136‑1 (MQ=255)
                                                                                                             tGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGc                               >  1:229447/1‑151 (MQ=255)
                                                                                                                     agaTAAC‑TTC‑CC‑GCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGt                         <  1:654290/151‑1 (MQ=255)
                                                                                                                       aTAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTAc                       >  1:828991/1‑151 (MQ=255)
                                                                                                                         aaC‑TTC‑CC‑GCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTAc                       <  1:968080/149‑1 (MQ=255)
                                                                                                                          aC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACaa                     >  1:648540/1‑150 (MQ=255)
                                                                                                                              ttccc‑gccgCTAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACaaa                    <  1:26318/146‑1 (MQ=255)
                                                                                                                              tt‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCATACAAACgg                 >  1:50709/1‑151 (MQ=255)
                                                                                                                               t‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg                >  2:759645/1‑151 (MQ=255)
                                                                                                                                 cc‑gccgCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg                >  1:498645/1‑150 (MQ=255)
                                                                                                                                    gccgcTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtat             >  2:579701/1‑151 (MQ=255)
                                                                                                                                     ccgcTGGCCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtata            <  2:168551/151‑1 (MQ=255)
                                                                                                                                     ccgcTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtat             >  1:916864/1‑150 (MQ=255)
                                                                                                                                     ccgcTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtata            <  1:84011/151‑1 (MQ=255)
                                                                                                                                      caatAGCCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACCCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtata            >  2:1032699/4‑150 (MQ=255)
                                                                                                                                       gcTAGTCGCAAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATag           <  2:281438/150‑1 (MQ=255)
                                                                                                                                         tAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggc        >  1:40320/1‑151 (MQ=255)
                                                                                                                                             cgcgAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGCATagggcaagg    >  2:42208/1‑151 (MQ=255)
                                                                                                                                                 aaaCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggcaagggc  <  1:351959/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGGCCCCAA‑CTCCAGCGCGAGGCATAAG‑CCTTTTTGCGT‑G‑AGATAAC‑TTT‑CC‑GCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGC  >  NC_000913/3511784‑3512069

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