Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,511,9230AC46.6% 108.0 / 15.2 73V141G (GTT→GGT) yhfZputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (32/7);  new base C (0/34);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 3.07e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCATTG  >  NC_000913/3511773‑3512073
                                                                                                                                                          |                                                                                                                                                      
gAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGaaa                                                                                                                                                        <  1:828497/151‑1 (MQ=255)
 aaCGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAac                                                                                                                                                       >  1:986232/1‑151 (MQ=255)
      ttATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAaccacc                                                                                                                                                   <  1:224381/150‑1 (MQ=255)
           cAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAACCCTTTTT‑GCGTGAGATAACTTC‑CCGCCGCTAGTCGCGAAACCACCGCCATg                                                                                                                                             <  1:130029/151‑1 (MQ=255)
              ccccACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATACCTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCa                                                                                                                                          <  1:716233/151‑1 (MQ=255)
              ccccACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTc                                                                                                                                           >  2:1284634/1‑150 (MQ=255)
               cccACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCAt                                                                                                                                         <  2:124043/151‑1 (MQ=255)
                  aCGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATac                                                                                                                                       >  2:626665/1‑150 (MQ=255)
                    gcgcTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACg                                                                                                                                    >  1:1362399/1‑151 (MQ=255)
                       cTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGcc                                                                                                                                  >  1:325824/1‑150 (MQ=255)
                          ccctttggcgcccccGCTTTTACGAAAAATCACTGGGTG‑CTCCCAAACGTGGGTGTGCGC‑CCCCACCCCCGGCGCGAGGCAAAACCTTTTTG‑GCGTGAAAAACCTTC‑CCCCCGCTAGCCGCAAACCCCCCGCCATGTCATACACGCCAtt                                                                                                                               <  2:736283/136‑1 (MQ=255)
                          cACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCAtt                                                                                                                               >  1:1006080/1‑150 (MQ=255)
                             gTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGc                                                                                                                           >  2:978633/1‑151 (MQ=255)
                               tgggggcccccccTTTCCGGCAAATCCA‑TGGTGCCTCGCAAGCGTGGGTGTGCGC‑CCCCACCCCCGGCGCGGGGCATAACCTTTTTT‑GCGTGAGATACCTTT‑CCCCCGCTATTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCa                                                                                                                          <  1:294653/139‑1 (MQ=255)
                                   ggtcTCGCTTTTCCGCAAACTCACCGGGTC‑CTCGCCACCGTGGGTGTGCGC‑CCCCACCCCCGGGGGGGGAAATACCCTTTTTT‑GCGGGAGATACTTTC‑CCGCCGCTGGGCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGac                                                                                                                       <  1:23757/146‑1 (MQ=255)
                                   ggACCCCCCTTTCCGCCAAACAAACTGGTC‑CCCCCCAAGGTGGGTGGGCGC‑CCCAACCTCCAGCGGGGGGAATAACCTTTTTG‑GCGTGAAATACCTTC‑CCGCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAg                                                                                                                         <  1:394004/147‑1 (MQ=255)
                                    gtcccGCTTTTGCGACAAACCAACTGGTC‑CCCGCCACCGTGGGTGTGCGC‑CCCCAACCCCAGCGCGAGGCATAACCCTTTTG‑GCTTGAGAAAGTTTC‑CCGCCGTTAGTCGCAAAACCCCCGCCATGTCATACACGCCATTAAGCa                                                                                                                          <  1:1252874/141‑1 (MQ=255)
                                    gACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACt                                                                                                                    >  2:553791/1‑151 (MQ=255)
                                       tCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTcc                                                                                                                  >  1:566608/1‑150 (MQ=255)
                                       tCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTcc                                                                                                                  >  2:616637/1‑150 (MQ=255)
                                       tCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTcc                                                                                                                  >  2:67611/1‑150 (MQ=255)
                                            tttACGGCAAATCAACTGGTT‑CTAGCCAACGTAGGTGTGCGG‑CCACAGCTACAGCGCGAGGCATACGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGACGCGAGACCACCGCCATGTTATATAAGCCATTTCGCAGACACTCCa                                                                                                                 >  1:245527/1‑146 (MQ=255)
                                            tttACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCAcgcg                                                                                                             >  2:818174/1‑150 (MQ=255)
                                             ttACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGa                                                                                                            >  2:733064/1‑150 (MQ=255)
                                                cGGAAAACCACCTGGTG‑CTCGCCACCGTGGGTGTGCGC‑CCCCACCCCCAGCGCGGGCAATAAGCTTTTTT‑GCGTAAAATACCTTC‑CCCCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAat                                                                                                          <  1:132536/149‑1 (MQ=255)
                                                 ggaaattcactggggg‑CTCGCAACCGTAGGTGTGGGC‑CCCCAACTCCGGCGCGGGGCATACCCCTTTTT‑GCGTGAAATAACTTC‑CCGCCGCTAGTCGCGAAACCACCGCAATGTCATACACGCCATTAAGCAGACACTCCACGCGAatat                                                                                                        <  2:1065051/135‑1 (MQ=255)
                                                 ggCAAATAAACTGGTG‑CTCGCCAACGTAGGTGTGCGC‑CCCCAACCCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATACCTTC‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAatat                                                                                                        <  2:1124543/150‑1 (MQ=255)
                                                    aaaaCAACGGGTG‑CCCCCCAACGTAGGTGGGCGC‑CCCCACCCCCAGCCCGAGGCATAACCTTTTTG‑GCTTGAGAAAACTTC‑CCGCCGCTGGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCcgc                                                                                                    <  1:28866/151‑1 (MQ=255)
                                                    aaaCCAAGTGTTG‑CTCCCCAACGTGGGTTGCCGG‑CCCAAACCCCGGCGCGCGGAATAACCCTTTTG‑GCGTGAAATAATTTC‑CCGCCGCTCGCCGCGAACCCACCGCAATGTCATATACGCCATTAAGCAGACACTCCACGCGAATATcc                                                                                                      <  1:1078714/149‑1 (MQ=255)
                                                    aaaCCAACTGGTG‑CTCGCCAACGTAGGTGTGCGC‑CCCCACCCCCAGCGCGGGGCATAAGCTTTTTT‑GCGTGAGAAAACTTC‑CCCCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCcgc                                                                                                    <  2:238121/151‑1 (MQ=255)
                                                              ggtccccGCCAACGTAGTGGTGCGC‑CCCCAACCCCAGCGCGAGGCATAACCTTTTTT‑GCGGGAAATACTTTC‑CCGCCGCTAGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACg                                                                                              <  1:52256/142‑1 (MQ=255)
                                                              gTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtg                                                                                          >  2:693590/1‑151 (MQ=255)
                                                                     gccacggtggggtggcgc‑CCCAAACTCCAGCGCGGGGCAAAAGCCTTTTG‑GCGGGAGATACTTTC‑CCCCCGCTGGTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgtg                                                                                        <  1:2845/133‑1 (MQ=255)
                                                                      ccAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   >  1:252180/1‑151 (MQ=255)
                                                                      ccAACGT‑GGGGTGGGGCCCCCACCCCCGGCGCGGGGAAAAAGCTTTTTG‑GGGGGAGATAACTTC‑CCGCCGCTAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCAtgtg                                                                                        <  1:1355908/146‑1 (MQ=255)
                                                                        aaCGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   >  2:426704/1‑149 (MQ=255)
                                                                        aaCGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   >  1:615132/1‑149 (MQ=255)
                                                                        aaCGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   >  2:511095/1‑149 (MQ=255)
                                                                         aCGTGGGTGTGCGC‑CCCCACCCCCAGCGCGGGCCATAAGCTTTTTG‑GCGGGAAATAACTTC‑CCGCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCAt                                                                                   <  1:410725/148‑1 (MQ=255)
                                                                           gaaGGTGCGCGC‑CCCAAACCCCGGCGCGGGGCATATCCCTTTTG‑GCGAGAGAAAACTTC‑CCGCCGCTGGGCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTa                                                                               <  1:314760/148‑1 (MQ=255)
                                                                             aGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaa                                                                             >  2:256901/1‑150 (MQ=255)
                                                                              ggTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaaa                                                                            >  2:248729/1‑150 (MQ=255)
                                                                              ggTGTGCGC‑CCCAAACCCCAGCGCGGGGCATAAGCTTTTTT‑GCGTGAGATACCTTC‑CCCCCGCTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaa                                                                             <  2:734904/149‑1 (MQ=255)
                                                                              ggTGGGCGC‑CCCCACCCCCAGGGCGAGGAATAGCCTTTTTT‑GCGTGAAATAACTTC‑CCGCCGCTATTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTaaa                                                                             <  1:841875/149‑1 (MQ=255)
                                                                                      gcccccACCCCCAGCGCGAGGAAAAAGCTTTTTT‑GGGTGAGATACTTTT‑CCGCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGc                                                                    <  2:495614/149‑1 (MQ=255)
                                                                                        ccccACCCCCAGCGCGGGGAAAAAGCTTTTTG‑GCGTGAAAAAACTTC‑CCGCCGCTAGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGaa                                                                       <  1:133992/146‑1 (MQ=255)
                                                                                            acccccAGCGCGAGGCATAACCTTTTTT‑GCGTGAGATACCTTC‑CCGCCGCTAGCCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCa                                                               <  1:668775/146‑1 (MQ=255)
                                                                                            aaCCCCAGCGCGAGGCATAAGCTTTTTT‑GCGTGAAATACTTTC‑CCGCCGCTAGTCGCGAACCCCCCGCCATGTCATCCACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATc                                                                <  2:940683/149‑1 (MQ=255)
                                                                                              cTCCAGCGCGGGCCATAACCTTTTTG‑GCGTAAAATAACTTT‑CCGCCTCAAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCaaa                                                             <  1:840547/150‑1 (MQ=255)
                                                                                               cccAGCGCGGGGCATAACCTTTTTG‑GCGTGAGATACCTTC‑CCCCCGCTGGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACt                                                           <  2:14466/150‑1 (MQ=255)
                                                                                               cccAGCGCGAGCCA‑AAGCCTTTTTGGGGGGAAATACCTTC‑CCGCCGCTAGCCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACt                                                           <  2:860560/150‑1 (MQ=255)
                                                                                                ccggcgcGGGCCATAACCTTTTTT‑GGGTGAGATAATTTC‑CCCCCGCTGGCCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTg                                                          <  1:1281843/148‑1 (MQ=255)
                                                                                                    cgcgAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCt                                                      >  1:351226/1‑151 (MQ=255)
                                                                                                     gcgAGGCATAACCTTTTTT‑GCGTGAGATAACTTC‑CCGCCGCTAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCtt                                                     <  2:1255612/151‑1 (MQ=255)
                                                                                                         ggaaTAACCCTTTTT‑GGGTAAAATACCTTC‑CCGCCGCTAGCGGCAAAACCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCttt                                                    <  1:666604/145‑1 (MQ=255)
                                                                                                          ggaataGCCTTTTTGGCGT‑AAAAAACTTTCCCCCCGTTAGCCGGAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCttt                                                    <  1:1072934/143‑1 (MQ=255)
                                                                                                          gCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAAt                                                >  2:1221075/1‑151 (MQ=255)
                                                                                                           cATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCGATc                                               >  1:907393/1‑151 (MQ=255)
                                                                                                              aaGCTTTTTT‑GCGTGAGATAATTTC‑CCGCCGCTGGGCGCAAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATcc                                              <  2:736687/149‑1 (MQ=255)
                                                                                                                  cTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCg                                         >  2:909592/1‑150 (MQ=255)
                                                                                                                  cTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGc                                        >  2:894705/1‑151 (MQ=255)
                                                                                                                          cGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACCCTCCACGCGAATATCCGCGCCCCGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCt                                 >  1:499900/1‑151 (MQ=255)
                                                                                                                               aaataCTTTCCCGCCGCAAGTCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCtt                                <  1:305206/144‑1 (MQ=255)
                                                                                                                                aTAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTAc                           >  1:41278/1‑151 (MQ=255)
                                                                                                                                 tAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTAc                           >  2:288466/1‑150 (MQ=255)
                                                                                                                                  accTTC‑CCGCCGCTAGCCCCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTAc                           <  1:567879/147‑1 (MQ=255)
                                                                                                                                  accTTC‑CCGCCCCTAGCCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTc                               <  1:279154/143‑1 (MQ=255)
                                                                                                                                      tt‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACgg                     >  1:749132/1‑151 (MQ=255)
                                                                                                                                         ccgccgCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGt                   >  2:1264164/1‑151 (MQ=255)
                                                                                                                                         ccgccgCGACTCGCAAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACggg                    <  2:779821/150‑1 (MQ=255)
                                                                                                                                         ccgccgCAAGCCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACgg                     <  2:278189/149‑1 (MQ=255)
                                                                                                                                             cgcTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGtat                 >  2:944486/1‑149 (MQ=255)
                                                                                                                                             cccTATTCGCGAACCCCCCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATag               <  2:379487/149‑1 (MQ=255)
                                                                                                                                              gcTAGTCGCGAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACaaa                        <  2:1054946/141‑1 (MQ=255)
                                                                                                                                              gcTAGCCGCAAACCCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagg              <  2:1237219/151‑1 (MQ=255)
                                                                                                                                                    cgcgAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATagggcaagg        >  1:1230127/1‑151 (MQ=255)
                                                                                                                                                          accaccGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCATTg  >  1:527652/1‑151 (MQ=255)
                                                                                                                                                          |                                                                                                                                                      
GAACGGTTATCCAGCCCCACGCGCTTCACGTTTGCGGACTCGCCTTTACGGCAAATCAACTGGTG‑CTCGCCAACGTAGGTGTGCGG‑CCCCAACTCCAGCGCGAGGCATAAGCCTTTTT‑GCGTGAGATAACTTT‑CCGCCGCTAGTCGCGAAACCACCGCCATGTCATACACGCCATTAAGCAGACACTCCACGCGAATATCCGCGCCACGCATGTGCGCATAGTAAAAAGGAATGCCATCAAACTGGGCTTTCAATCCGCTCGCCAGGCCTTCGTACAAACGGGTATAGGGCAAGGGCATTG  >  NC_000913/3511773‑3512073

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