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

GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001770‑2002067
                                                                                                                                                       |                                                                                                                                                      
gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑cc                                                                                                                                                     <  2:7123/151‑1 (MQ=255)
gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAACGGCCTGATTATCCGTTTCTGCAAGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGaaaaaaa                                                                                                                                                        <  1:587184/149‑1 (MQ=255)
   gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCccg                                                                                                                                                  <  1:400280/151‑1 (MQ=255)
     ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCccgcc                                                                                                                                                >  1:639273/1‑151 (MQ=255)
      taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCccgccg                                                                                                                                               >  2:544877/1‑151 (MQ=255)
               ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCgggg                                                                                                                                      >  1:381125/1‑151 (MQ=255)
                 cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGa                                                                                                                                     >  2:451844/1‑150 (MQ=255)
                     aaaaaCAGCCCTGCGTTATATGATTTATGGGCATGATTATCCGTTTCTGCAGGTTTTTA‑AATCGGACGATTGGGGGGTAAAATGGGGGTTTATTGGGGGG‑TTAAA‑GGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAg                                                                                                                                   <  2:721009/146‑1 (MQ=255)
                      aaaccaccccTGCTTTATTTGGGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTT‑AATCGGACGATTAGTGGGTAAAAGGGGGGTTTTTTTGGGGG‑TTACA‑GATAATTCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACg                                                                                                                               <  1:99413/144‑1 (MQ=255)
                       aaaCACCCGTGGTTTATATGAGTTATCGGCATAATTATCCGTTTCTGCGGGGTTTTA‑AATCGGACGATTGGTGGGGGAATTGAGGGGTTTTTGGGGGG‑TTACA‑GGTAAATCCCAGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCAc                                                                                                                                <  2:525268/149‑1 (MQ=255)
                             cccTGCGTTATATGATTTAGCGGCATGTTTTTCCGTTTCTGGGGGTTTTTT‑AATCGGACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGG‑TTACA‑GGTAAATCCCGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                         <  1:334622/150‑1 (MQ=255)
                              cctggtttATATGAGTTACCGGCAGGTTTACCCGTTTCTGGGGGTTTTTT‑AATCGGACGATTGGTGGGTAAAAGGGGGGGTTTTTTGGGGT‑TAAGA‑GGAAAATCCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                       <  1:399763/145‑1 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  1:474799/1‑150 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  1:641581/1‑150 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  1:329302/1‑150 (MQ=255)
                                    ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                 >  1:248686/1‑151 (MQ=255)
                                    ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                 >  2:248358/1‑151 (MQ=255)
                                    ttATATGAGTTACCGGAATGATTATCCTTTTCGGCGGGTTTTTT‑AACCGGAGGATTAGGGGGGGAAATGGGGGTTTATTGGGGGT‑TTACA‑GGTAAATCCAGGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      <  2:849180/146‑1 (MQ=255)
                                      atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                               >  1:886348/1‑151 (MQ=255)
                                           aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑ATACA‑GGCAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGCGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                           >  2:163252/1‑150 (MQ=255)
                                             ttATCGGCATGATTATCCTTTTCGGCGGGGTTTTT‑AAGCGGACGTTTAGTGGGTGAAA‑GGGGGGTTTTTTGGGGGTTTACA‑GGTAAATTCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                          <  2:121159/149‑1 (MQ=255)
                                               aTCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                      >  1:924231/1‑151 (MQ=255)
                                                  ggCATGTTTATCCGTTTCTGCGGGTTTTTT‑AACCGGACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGG‑TTACA‑GTTAAATTCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                     <  2:847462/149‑1 (MQ=255)
                                                  ggCATGATTTTCCTTTTCTGCGGGGTTTTT‑AATCGGACGATTGGGGGGTAAAAGGGGGGTTTATTTGGGGT‑TAACA‑GGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                     <  1:537812/149‑1 (MQ=255)
                                                      tGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                                >  2:769578/1‑150 (MQ=255)
                                                         ttctCCGTTTCTGAAGGTTTTTT‑AATCGGACGATTAGTGGATAAAATGAGGGTTTTTTTGGGGT‑TGAAG‑GGTAAGTTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                              <  2:299249/146‑1 (MQ=255)
                                                         ttATCCTTTTCTGCGGGTTTTTA‑AACCGGACGATTGGGGGGTAAAATGAGGGTTTA‑TTGGGGG‑TTAAAGGGTAATTTCAAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                 <  2:650846/146‑1 (MQ=255)
                                                         ttATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                            >  1:547821/1‑151 (MQ=255)
                                                         ttATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                            >  1:260163/1‑151 (MQ=255)
                                                         ttACCCTTTTCGGCGGGGTTTTT‑AATCGGACGATTGGGGGGTAAAATGGGGGTTTTTTTGGGGT‑TAACA‑GAAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                  <  1:441995/145‑1 (MQ=255)
                                                           aTCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                          >  2:713774/1‑151 (MQ=255)
                                                           aTCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAAGGAGGGTTTTTTTGGGGG‑TTACA‑GGAAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                           <  1:4220/150‑1 (MQ=255)
                                                           aTCCGTTGCGGGAGTGTTTTG‑AAGAGCAGGATTGGTGGGTAAAAGGAGGGTTTTTTGGGGGT‑TTACA‑GCTAAATCCCCCGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                               <  2:776548/146‑1 (MQ=255)
                                                              cGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                        >  1:300022/1‑150 (MQ=255)
                                                                  tgggtaGGGTTTTG‑AATGGTCTGTTTAGGGGGTAAAACG‑GGGGTTTTTTTGGGG‑TTAAA‑ATTGAAATCCCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                            <  1:966766/137‑1 (MQ=255)
                                                                  tcggcggggTTTTA‑AACCGGGGGATGGGGGGGGGAAAGGGGGGTTTATTTGGGGG‑TGACA‑GGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                    <  2:498467/144‑1 (MQ=255)
                                                                   cTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                  >  2:457132/1‑151 (MQ=255)
                                                                      gcgggTTTTTAAATCGGACGATTGGGGGGTAAAAGGGGGGTTTTTTTGGGGG‑TAACA‑GGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                <  1:609853/149‑1 (MQ=255)
                                                                      cAGGGTTTTT‑AATCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTTGGGGG‑TTACA‑GGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               <  1:317645/151‑1 (MQ=255)
                                                                      cAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  2:367708/1‑151 (MQ=255)
                                                                       aGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  2:830087/1‑150 (MQ=255)
                                                                          tttttt‑AATCGGACTTTTAGTGGGTAAAATGAGGGTTTTTTTGGGGT‑TTAGA‑GATAATTTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                            <  1:599426/149‑1 (MQ=255)
                                                                           ttttt‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                          >  2:821215/1‑151 (MQ=255)
                                                                           ttttt‑AACCGGACGATTGGGGGGTGAAAGGGGGGTTTTTTTGGGGG‑TTAAA‑GGAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               <  1:161038/146‑1 (MQ=255)
                                                                            tttt‑AATCGGACGATTGGTGGGTAAAAGGGGGGTTTTTTTGGGGG‑TTACA‑GTAAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                          <  1:843248/150‑1 (MQ=38)
                                                                             ttt‑AATGGGACGATTAGTGGGTAAAATGAGGGGTTATTGGGGGG‑TTACA‑GGTAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                         <  2:299360/150‑1 (MQ=39)
                                                                             ttt‑AATCGGACGATTAGGGGGTGAAAGGAGGGTTTATTGGGGGG‑TTACA‑GGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                         <  2:757256/150‑1 (MQ=39)
                                                                                 aagcGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGT‑TTACA‑GATAAATTCCAGGAAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                      <  2:678139/147‑1 (MQ=37)
                                                                                  aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                    >  1:627950/1‑151 (MQ=255)
                                                                                   tCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                    >  2:422457/1‑150 (MQ=255)
                                                                                     ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                                  >  1:180947/1‑150 (MQ=255)
                                                                                     ggACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGG‑TTACA‑GGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                                  <  2:404925/150‑1 (MQ=38)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                              >  1:575836/1‑151 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                           >  1:542941/1‑151 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                           >  2:550294/1‑151 (MQ=255)
                                                                                              gTGGGTAAAAGGGGGGTTTTTTGGGGGG‑TTACA‑GGTAAATTCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                          <  2:160128/149‑1 (MQ=37)
                                                                                                     aaaggggggtTTATTGGGGGG‑TTACA‑GGAAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  2:549919/143‑1 (MQ=37)
                                                                                                     aaaTGGGGGGTTTTTTGGGGG‑TTACA‑GATAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  1:137392/149‑1 (MQ=37)
                                                                                                      aaTGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                 >  2:785939/1‑150 (MQ=255)
                                                                                                               ttttttGGGGT‑TTAGA‑TAAAAATTCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                         <  1:161911/146‑1 (MQ=25)
                                                                                                               ttttttGGGGT‑TTAAA‑GATAAATCCCAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                        <  1:431480/147‑1 (MQ=25)
                                                                                                                tATTGGGGGG‑TTAAA‑GGTAAATTCCGGGCAGAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                        <  2:71035/149‑1 (MQ=255)
                                                                                                                     ggggg‑TTAAA‑GGTAAA‑TCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                   <  2:639284/149‑1 (MQ=37)
                                                                                                                      gggg‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                                >  2:318873/1‑151 (MQ=255)
                                                                                                                       ggg‑TTACA‑GGTAAATCCCAGGCAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                               <  2:927969/151‑1 (MQ=255)
                                                                                                                         g‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                              >  2:559404/1‑150 (MQ=255)
                                                                                                                           ttACG‑GGAAAATTCCGGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                              <  2:686138/149‑1 (MQ=37)
                                                                                                                             aCA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                          >  2:182199/1‑151 (MQ=255)
                                                                                                                              cA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                         >  2:122775/1‑151 (MQ=255)
                                                                                                                                      aTTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                    >  1:763056/1‑149 (MQ=255)
                                                                                                                                         ccAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc               >  2:945109/1‑151 (MQ=255)
                                                                                                                                         ccAGGAAAAAAAAAC‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc               <  2:919108/151‑1 (MQ=255)
                                                                                                                                                gAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt        >  2:131376/1‑151 (MQ=255)
                                                                                                                                                  aaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa      <  2:743408/151‑1 (MQ=255)
                                                                                                                                                  aaaaac‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa      <  1:959238/151‑1 (MQ=255)
                                                                                                                                                       a‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa  >  2:382728/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTT‑AATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACA‑GGTAAATTCCAGGCAGAAAAAAA‑CCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001770‑2002067

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