Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9110GA41.0% 89.6 / 53.5 83intergenic (‑122/+199)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base G (41/8);  new base A (0/34);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.40e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001763‑2002059
                                                                                                                                                             |                                                                                                                                                     
acCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGa                                                                                                                                                      >  1:107296/1‑150 (MQ=255)
     cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGaaaaaa                                                                                                                                                 >  1:71799/1‑150 (MQ=255)
      aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑c                                                                                                                                              >  1:513806/1‑151 (MQ=255)
      aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑c                                                                                                                                              >  1:290403/1‑151 (MQ=255)
        cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑ccc                                                                                                                                            <  1:901876/151‑1 (MQ=255)
           tataGAGTGAATTCACGATAAGCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCccgc                                                                                                                                         >  2:565653/1‑151 (MQ=255)
           tataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCccgc                                                                                                                                         >  2:128919/1‑151 (MQ=255)
            ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGT‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCccgcc                                                                                                                                        >  1:1094218/1‑151 (MQ=255)
             taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCccgccg                                                                                                                                       >  1:563399/1‑151 (MQ=255)
              agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCccgccg                                                                                                                                       >  2:757247/1‑150 (MQ=255)
              agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCgg                                                                                                                                      <  1:1038091/151‑1 (MQ=255)
                  tGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTgg                                                                                                                                   >  2:660315/1‑150 (MQ=255)
                   gAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCg                                                                                                                                 >  1:772945/1‑151 (MQ=255)
                      ttCACGATAAACAGCCCTGTTTTATATGAGTTATCGGCATGTTTATCCTTTTC‑GGCGGGGTTTTTAATCGA‑ACGATGGGGGGTTAAAATGGGGGT‑TTTTTT‑GGGGT‑TTACAGT‑TAAATT‑CC‑AGG‑AAGAAAAAAA‑CCCCGCCGGTg                                                                                                                                    <  2:365247/145‑1 (MQ=255)
                            aTAACCGGCCCTGCTTTATTTGATTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGGGGGTGAAAGGGGGGT‑TTTTTT‑GGGGG‑TAGCAGG‑AAAATC‑CC‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                           <  2:186766/148‑1 (MQ=255)
                               accaccccTGCTTTATATGAGTTATCGGCATGATTTCCCTTTTC‑TGCAGGTTTTTTAACCGG‑ACGATTGGTGGGTGAAAGGGGGGT‑TTTTTT‑GGGGG‑TTAAAGG‑TAATTT‑CA‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                     <  2:781876/146‑1 (MQ=255)
                                aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                     >  2:1148437/1‑150 (MQ=255)
                                aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                    >  1:593920/1‑151 (MQ=255)
                                    cccTGCTTTATATGAGTTATCGGCATGATTTTCCTTTTC‑TGCAGGTTTTTTAATCGG‑ACGATTAGTGGGTAAAATGGGGGT‑TTATTG‑GGGGT‑TTACAGA‑AAAATC‑CC‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                 <  1:628887/150‑1 (MQ=255)
                                    cccTGCGTTATATGAGTTACCGGCATGATTATCCTTTTC‑GGCGGGTTTTTAAATCGG‑ACTATTAGGGGGTGAAATGGGGGT‑TTATTT‑GGGGG‑TTACAGA‑TAAATT‑CC‑AGG‑AAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                 <  1:1142585/150‑1 (MQ=255)
                                        gctttATATTTGTTA‑CGGGATGAATTACCGTTTCGTGGAGGTTTTTAAATCGA‑ATAATAAGGGGGAAAAATGGGGGG‑TT‑TTT‑GGGGG‑TTACAGG‑AAAATTCCC‑GGG‑CAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                             <  1:231270/147‑1 (MQ=255)
                                           ttatagggttgaTCGGATTGTTTATCCGTTTC‑TGAAGGTTTTTAAACGGG‑AGGATTAGGGGGTAAAATGAGGGG‑TTTTTT‑GGGGG‑TTGCAGG‑AAATTC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaaa                                                                                                            <  2:429387/137‑1 (MQ=255)
                                           ttATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                         >  2:66095/1‑151 (MQ=255)
                                            tataTGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACCAATTgc                                                                                                        >  1:631456/1‑151 (MQ=255)
                                                 gAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   >  1:1151604/1‑151 (MQ=255)
                                                  aGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                   >  1:1039149/1‑150 (MQ=255)
                                                    ttATGGGAATGTTTACCCTTTTC‑TGCGGGTTTTTTAATGGG‑ACGTTTAGGGGGTGAAATG‑GGGG‑TTTTTT‑GGGGGTTTAAAGT‑AAAATC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                    <  1:93267/147‑1 (MQ=255)
                                                    ttATCGGCATGATTATCCTTTTC‑TGGGGGTTTTTTAACGGG‑ACGTTGAGGGGGAAAAAGGGGGGT‑TTTTTT‑GGGGT‑TGAAAGA‑AAAATC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                      <  2:815926/145‑1 (MQ=255)
                                                    ttATCGGAATTTTTATCCTTTTC‑TGCGGGTTTTTTAATCGT‑ACGATTAGTGGGTAAAAGGGGGGT‑TTTTTT‑GGGGG‑TAACAGA‑AAAATC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                 <  2:1094147/150‑1 (MQ=255)
                                                            aTGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         >  2:112730/1‑150 (MQ=255)
                                                             tatttATCTTTTTC‑TGCGGGTTTTTAAAGGGG‑ACGATTGGGGGTTAAAAGGGGGGT‑TTATTG‑GGGGT‑TTAAAGT‑AAAATT‑CG‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                          <  1:714440/145‑1 (MQ=255)
                                                             tGATTTTCCTTTTC‑T‑CGGGGTTTTTAAAGGGAACGTTAGGTGGGTAAAAGGAGGGT‑TTTTTT‑GGGGG‑TAACAGG‑TAAATC‑CG‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                         <  2:339574/149‑1 (MQ=255)
                                                              gATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                        >  2:298773/1‑149 (MQ=255)
                                                                ttATCCGTTTC‑TGCGGGGTTTTTAATCGG‑ACGATTAGGGGGTAAAATGGGGGT‑TTTTTT‑GGGGG‑TTACAGT‑AAAATT‑CC‑AGG‑AAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                     <  1:729555/150‑1 (MQ=255)
                                                                 tATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                   >  2:950035/1‑151 (MQ=255)
                                                                      gTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTa                                                                               >  1:748618/1‑150 (MQ=255)
                                                                       tttC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                              >  2:1166372/1‑150 (MQ=255)
                                                                              cAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       >  1:104394/1‑151 (MQ=255)
                                                                               ggggTTTTTAATCGG‑‑CTATTTGGTGGGAAAAAGAGGGGTTTTTTT‑GGGGG‑TAACAGA‑AAATTC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                        <  2:663545/148‑1 (MQ=255)
                                                                               aGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                       >  1:52480/1‑150 (MQ=255)
                                                                                 ggTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                    >  2:251576/1‑151 (MQ=255)
                                                                                       tAATCGG‑ACGATTAGTGGGTAAAATGGGGGT‑TTATTG‑GGGGT‑TTACAGG‑TAAATC‑CC‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                               <  1:480042/150‑1 (MQ=38)
                                                                                        aaTGGG‑AAGATTGGGGGATAAATTGAGGGT‑TTTTTT‑GGGGG‑TTATAGG‑TAATTC‑CC‑AGG‑AAAAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                              <  2:336738/150‑1 (MQ=25)
                                                                                         aTCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                            >  2:792737/1‑151 (MQ=255)
                                                                                             gaatgATTAGTGGGTGAAAGGAGGGT‑TTTTTG‑GGGGT‑TTACAGG‑TAAATC‑CC‑GGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                          <  1:879528/146‑1 (MQ=37)
                                                                                              taCGTTTGGTGGATGAAAAGAGGGG‑TTGTTG‑GGGAG‑TTACGGG‑TATAAC‑CG‑GGG‑AAGAAAAAAT‑CGCCGCCGGGGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                              <  2:1036984/144‑1 (MQ=21)
                                                                                                cGATTAGTGGGTGAAATGAGGGG‑TTATTTGGGGGG‑TTACAGG‑AAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                        >  1:142619/1‑150 (MQ=255)
                                                                                                  aTTAGGGGGTGAAATGAGGGT‑TTTTTT‑GGGGG‑TTAAAAG‑TAAATC‑CG‑AGG‑AAAAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                       <  2:39004/148‑1 (MQ=37)
                                                                                                     gttgggaaaaaTGGGTGG‑TCATTT‑GGGGG‑TAAGAGG‑TAAGTT‑CA‑AGC‑CAAAAAAAAA‑CCCCGCCGTTGGCGGGTAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGTATTCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                     <  2:553066/139‑1 (MQ=25)
                                                                                                     aGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  1:926973/1‑150 (MQ=255)
                                                                                                     aGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                  >  2:852455/1‑150 (MQ=255)
                                                                                                         ggTAAAATG‑GGGG‑TTTTTT‑GGGGG‑TTTAAGGATAAATC‑CC‑GGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                             <  1:10731/151‑1 (MQ=25)
                                                                                                           tGAAATGAGGGT‑TTTTTG‑GGGGG‑TTAAGGT‑TAAATC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                             <  1:1054343/149‑1 (MQ=25)
                                                                                                              aaagggggg‑TTTTTT‑TGGGG‑GTAAAGG‑TAAAAT‑CCAAGG‑CAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                           <  2:862077/144‑1 (MQ=25)
                                                                                                              aaTGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                         >  2:359764/1‑150 (MQ=255)
                                                                                                                  gaggg‑TTTTTT‑TGGGG‑GTAAAGG‑ATAATT‑CC‑CGG‑GAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                       <  2:38993/147‑1 (MQ=25)
                                                                                                                  aGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                    >  1:686584/1‑151 (MQ=255)
                                                                                                                   gggt‑TTTTTT‑GGGGG‑TCACAGT‑AAAATT‑CA‑AGG‑CAAAAAAAAC‑CGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                   <  2:207999/151‑1 (MQ=25)
                                                                                                                    ggg‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  >  1:748139/1‑151 (MQ=255)
                                                                                                                    ggg‑TTATTG‑GGGGT‑TAACAGG‑AAAATC‑CC‑GGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                  <  1:769253/151‑1 (MQ=25)
                                                                                                                        ttATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                               >  2:845684/1‑151 (MQ=255)
                                                                                                                        ttATTG‑GGGGG‑TTACAGG‑AAAATC‑CC‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAATTGTGCAAATTCAATAAATTGCaa                                <  1:876591/150‑1 (MQ=37)
                                                                                                                         tATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                              >  2:327005/1‑151 (MQ=255)
                                                                                                                          tttt‑GGGGGTTTACAGG‑TAAATC‑CC‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAt                               <  1:903370/149‑1 (MQ=37)
                                                                                                                             tgggggg‑TTACAGG‑TAAATC‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                            <  1:650692/148‑1 (MQ=255)
                                                                                                                               ggggg‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                         >  1:1101374/1‑151 (MQ=255)
                                                                                                                                ggggatttaaGG‑TAAATA‑CG‑GGT‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATATTGCCGTATGCGGCGTAAACGCCTTATCCGGCATACAAAAATGTGCAAATTCAATAAATCGCAATTCAACt                         <  2:930663/142‑1 (MQ=25)
                                                                                                                                gggg‑TTACAGA‑TAAATT‑CC‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                        <  2:166202/151‑1 (MQ=255)
                                                                                                                                 gggggtaagag‑TAAAAT‑CC‑AGG‑GAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACt                         <  2:1154782/141‑1 (MQ=25)
                                                                                                                                 ggg‑TTAAGGG‑TAAATC‑CA‑AGG‑CAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                       <  1:630852/151‑1 (MQ=37)
                                                                                                                                      tACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   >  1:25253/1‑151 (MQ=255)
                                                                                                                                       aaaGG‑AAAATC‑CC‑AGA‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   <  1:251584/148‑1 (MQ=37)
                                                                                                                                       aCAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                  >  1:1056044/1‑151 (MQ=255)
                                                                                                                                          aataaattt‑CC‑GGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                     <  2:156992/139‑1 (MQ=37)
                                                                                                                                             tAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa              >  1:16491/1‑150 (MQ=255)
                                                                                                                                              aaaTT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt             >  1:111465/1‑150 (MQ=255)
                                                                                                                                               aaat‑CC‑AGGACAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                   <  1:39979/141‑1 (MQ=255)
                                                                                                                                                aTT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa            >  2:179148/1‑149 (MQ=255)
                                                                                                                                                  tccc‑AGG‑CAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCAGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag          <  2:806482/148‑1 (MQ=255)
                                                                                                                                                    cc‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAGTTGCAATTCAACTTGTAGGCCTGATAagcgc       <  1:950094/151‑1 (MQ=255)
                                                                                                                                                    cc‑AGG‑AAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATaa           <  2:214075/147‑1 (MQ=255)
                                                                                                                                                         ggaaGAAAAACC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag     <  1:560686/147‑1 (MQ=255)
                                                                                                                                                           cAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcgc  >  1:606556/1‑151 (MQ=255)
                                                                                                                                                             |                                                                                                                                                     
ACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGC  >  NC_000913/2001763‑2002059

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