Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,001,9180AC41.2% 138.0 / 21.5 81intergenic (‑129/+192)fliA/fliCRNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein (flagellin)
Reads supporting (aligned to +/- strand):  ref base A (39/8);  new base C (0/33);  total (39/42)
Fisher's exact test for biased strand distribution p-value = 3.10e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065
                                                                                                                                                              |                                                                                                                                                    
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)
                                                                                                                                           aatgccgagtaaaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   <  2:780662/137‑1 (MQ=25)
                                                                                                                                            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)
                                                                                                                                                      aaaaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAag                <  2:281266/140‑1 (MQ=255)
                                                                                                                                                        aaaaacc‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag           <  1:837322/145‑1 (MQ=255)
                                                                                                                                                        aaaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagcg         <  2:421841/147‑1 (MQ=255)
                                                                                                                                                        aaaaaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc     >  1:356227/1‑151 (MQ=255)
                                                                                                                                                         aaaaac‑CACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa    <  2:416543/151‑1 (MQ=255)
                                                                                                                                                          aaaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATc     >  1:136280/1‑149 (MQ=255)
                                                                                                                                                           aaaa‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  2:835763/1‑151 (MQ=255)
                                                                                                                                                              |                                                                                                                                                    
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTC‑TGCAGGGTTTTTAATCGG‑ACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGG‑TAAATT‑CC‑AGG‑CAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065

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