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

GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  NC_000913/2001770‑2002067
                                                                                                                                                       |                                                                                                                                                      
gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                                     <  2:7123/151‑1 (MQ=255)
gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAACGGCCTGATTATCCGTTTCTGCAAGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                                        <  1:587184/149‑1 (MQ=255)
   gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                                  <  1:400280/151‑1 (MQ=255)
     ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgcc                                                                                                                                                >  1:639273/1‑151 (MQ=255)
      taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgccg                                                                                                                                               >  2:544877/1‑151 (MQ=255)
               ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                      >  1:381125/1‑151 (MQ=255)
                 cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                     >  2:451844/1‑150 (MQ=255)
                     aaaaaCAGCCCTGCGTTATATGATTTATGGGCATGATTATCCGTTTCTGCAGGTTTTTAAATCGGACGATTGGGGGGTAAAATGGGGGTTTATTGGGGGGTTAAAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAg                                                                                                                                   <  2:721009/146‑1 (MQ=255)
                      aaaccaccccTGCTTTATTTGGGTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGATAATTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACg                                                                                                                               <  1:99413/144‑1 (MQ=255)
                       aaaCACCCGTGGTTTATATGAGTTATCGGCATAATTATCCGTTTCTGCGGGGTTTTAAATCGGACGATTGGTGGGGGAATTGAGGGGTTTTTGGGGGGTTACAGGTAAATCCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                                <  2:525268/149‑1 (MQ=255)
                             cccTGCGTTATATGATTTAGCGGCATGTTTTTCCGTTTCTGGGGGTTTTTTAATCGGACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTg                                                                                                                         <  1:334622/150‑1 (MQ=255)
                              cctggtttATATGAGTTACCGGCAGGTTTACCCGTTTCTGGGGGTTTTTTAATCGGACGATTGGTGGGTAAAAGGGGGGGTTTTTTGGGGTTAAGAGGAAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGAc                                                                                                                       <  1:399763/145‑1 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  1:474799/1‑150 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  1:641581/1‑150 (MQ=255)
                                tGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      >  1:329302/1‑150 (MQ=255)
                                    ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                 >  1:248686/1‑151 (MQ=255)
                                    ttATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTg                                                                                                                 >  2:248358/1‑151 (MQ=255)
                                    ttATATGAGTTACCGGAATGATTATCCTTTTCGGCGGGTTTTTTAACCGGAGGATTAGGGGGGGAAATGGGGGTTTATTGGGGGTTTACAGGTAAATCCAGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACa                                                                                                                      <  2:849180/146‑1 (MQ=255)
                                      atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcg                                                                                                               >  1:886348/1‑151 (MQ=255)
                                           aGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGATACAGGCAAATTCCAGGCAGAAAAAAACCCCGCCGGCGGCGGGGAAGCACGTTGCTGACAAATTGCGCttt                                                                                                           >  2:163252/1‑150 (MQ=255)
                                             ttATCGGCATGATTATCCTTTTCGGCGGGGTTTTTAAGCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTGGGGGTTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                          <  2:121159/149‑1 (MQ=255)
                                               aTCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                      >  1:924231/1‑151 (MQ=255)
                                                  ggCATGTTTATCCGTTTCTGCGGGTTTTTTAACCGGACGATTAGGGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGTTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                     <  2:847462/149‑1 (MQ=255)
                                                  ggCATGATTTTCCTTTTCTGCGGGGTTTTTAATCGGACGATTGGGGGGTAAAAGGGGGGTTTATTTGGGGTTAACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTg                                                                                                     <  1:537812/149‑1 (MQ=255)
                                                      tGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                                >  2:769578/1‑150 (MQ=255)
                                                         ttctCCGTTTCTGAAGGTTTTTTAATCGGACGATTAGTGGATAAAATGAGGGTTTTTTTGGGGTTGAAGGGTAAGTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                              <  2:299249/146‑1 (MQ=255)
                                                         ttATCCTTTTCTGCGGGTTTTTAAACCGGACGATTGGGGGGTAAAATGAGGGTTTATTGGGGGTTAAAGGGTAATTTCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                                 <  2:650846/146‑1 (MQ=255)
                                                         ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                            >  1:547821/1‑151 (MQ=255)
                                                         ttATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                            >  1:260163/1‑151 (MQ=255)
                                                         ttACCCTTTTCGGCGGGGTTTTTAATCGGACGATTGGGGGGTAAAATGGGGGTTTTTTTGGGGTTAACAGAAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                  <  1:441995/145‑1 (MQ=255)
                                                           aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                          >  2:713774/1‑151 (MQ=255)
                                                           aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAAGGAGGGTTTTTTTGGGGGTTACAGGAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                           <  1:4220/150‑1 (MQ=255)
                                                           aTCCGTTGCGGGAGTGTTTTGAAGAGCAGGATTGGTGGGTAAAAGGAGGGTTTTTTGGGGGTTTACAGCTAAATCCCCCGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGAt                                                                                               <  2:776548/146‑1 (MQ=255)
                                                              cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                        >  1:300022/1‑150 (MQ=255)
                                                                  tgggtaGGGTTTTGAATGGTCTGTTTAGGGGGTAAAACGGGGGTTTTTTTGGGGTTAAAATTGAAATCCCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                            <  1:966766/137‑1 (MQ=255)
                                                                  tcggcggggTTTTAAACCGGGGGATGGGGGGGGGAAAGGGGGGTTTATTTGGGGGTGACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                    <  2:498467/144‑1 (MQ=255)
                                                                   cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                                  >  2:457132/1‑151 (MQ=255)
                                                                      gcgggTTTTTAAATCGGACGATTGGGGGGTAAAAGGGGGGTTTTTTTGGGGGTAACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                                <  1:609853/149‑1 (MQ=255)
                                                                      cAGGGTTTTTAATCGGACGTTTAGTGGGTGAAAGGGGGGTTTTTTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               <  1:317645/151‑1 (MQ=255)
                                                                      cAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  2:367708/1‑151 (MQ=255)
                                                                       aGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               >  2:830087/1‑150 (MQ=255)
                                                                          ttttttAATCGGACTTTTAGTGGGTAAAATGAGGGTTTTTTTGGGGTTTAGAGATAATTTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATc                                                                            <  1:599426/149‑1 (MQ=255)
                                                                           tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                          >  2:821215/1‑151 (MQ=255)
                                                                           tttttAACCGGACGATTGGGGGGTGAAAGGGGGGTTTTTTTGGGGGTTAAAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCtt                                                                               <  1:161038/146‑1 (MQ=255)
                                                                            ttttAATCGGACGATTGGTGGGTAAAAGGGGGGTTTTTTTGGGGGTTACAGTAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                          <  1:843248/150‑1 (MQ=38)
                                                                             tttAATGGGACGATTAGTGGGTAAAATGAGGGGTTATTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                         <  2:299360/150‑1 (MQ=39)
                                                                             tttAATCGGACGATTAGGGGGTGAAAGGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                         <  2:757256/150‑1 (MQ=39)
                                                                                 aagcGGACGATTAGGGGGTAAAAGGGGGGTTTATTGGGGGTTTACAGATAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                      <  2:678139/147‑1 (MQ=37)
                                                                                  aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                    >  1:627950/1‑151 (MQ=255)
                                                                                   tCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                                    >  2:422457/1‑150 (MQ=255)
                                                                                     ggACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                                  >  1:180947/1‑150 (MQ=255)
                                                                                     ggACGATTAGGGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaa                                                                  <  2:404925/150‑1 (MQ=38)
                                                                                        cGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                              >  1:575836/1‑151 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                           >  1:542941/1‑151 (MQ=255)
                                                                                           ttAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                           >  2:550294/1‑151 (MQ=255)
                                                                                              gTGGGTAAAAGGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                          <  2:160128/149‑1 (MQ=37)
                                                                                                     aaaggggggtTTATTGGGGGGTTACAGGAAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  2:549919/143‑1 (MQ=37)
                                                                                                     aaaTGGGGGGTTTTTTGGGGGTTACAGATAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                   <  1:137392/149‑1 (MQ=37)
                                                                                                      aaTGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                                 >  2:785939/1‑150 (MQ=255)
                                                                                                               ttttttGGGGTTTAGATAAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCa                                         <  1:161911/146‑1 (MQ=25)
                                                                                                               ttttttGGGGTTTAAAGATAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                        <  1:431480/147‑1 (MQ=25)
                                                                                                                tATTGGGGGGTTAAAGGTAAATTCCGGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCaa                                        <  2:71035/149‑1 (MQ=255)
                                                                                                                     gggggTTAAAGGTAAATCCAGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCaa                                   <  2:639284/149‑1 (MQ=37)
                                                                                                                      ggggTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                                >  2:318873/1‑151 (MQ=255)
                                                                                                                       gggTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                               <  2:927969/151‑1 (MQ=255)
                                                                                                                         gTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                              >  2:559404/1‑150 (MQ=255)
                                                                                                                           ttACGGGAAAATTCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                              <  2:686138/149‑1 (MQ=37)
                                                                                                                             aCAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                          >  2:182199/1‑151 (MQ=255)
                                                                                                                              cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                         >  2:122775/1‑151 (MQ=255)
                                                                                                                                      aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                    >  1:763056/1‑149 (MQ=255)
                                                                                                                                         ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc               >  2:945109/1‑151 (MQ=255)
                                                                                                                                         ccAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc               <  2:919108/151‑1 (MQ=255)
                                                                                                                                                gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt        >  2:131376/1‑151 (MQ=255)
                                                                                                                                                  aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa      <  2:743408/151‑1 (MQ=255)
                                                                                                                                                  aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa      <  1:959238/151‑1 (MQ=255)
                                                                                                                                                       aCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCa  >  2:382728/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
GCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCA  >  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: