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

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAT  >  NC_000913/2001764‑2002061
                                                                                                                                                     |                                                                                                                                                      
ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGaa                                                                                                                                                      >  2:251299/1‑150 (MQ=255)
  ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGaaaa                                                                                                                                                    >  1:490788/1‑150 (MQ=255)
   tCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGaaaaa                                                                                                                                                   <  1:1281217/150‑1 (MQ=255)
    cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGaaaaaaa                                                                                                                                                 <  2:1167315/151‑1 (MQ=255)
     aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAAc                                                                                                                                                >  1:225776/1‑151 (MQ=255)
      gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAAcc                                                                                                                                               <  2:1481258/151‑1 (MQ=255)
       cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAAccc                                                                                                                                              >  1:548539/1‑151 (MQ=255)
         gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccg                                                                                                                                            >  2:1186139/1‑151 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                           <  1:796517/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                           <  2:1285553/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                           <  2:122720/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                           <  1:135036/150‑1 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCGGGCAGAAAAAAACCccgcc                                                                                                                                          <  2:600435/150‑1 (MQ=255)
             agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                        <  1:314830/151‑1 (MQ=255)
                gtaaaTCCACGATAAACAGCCCTGTGTTTTAAGGGTTACGGGCAGGATTATCCGTTTCTGCAGGGTTTTTATTCGGACGATGAGGGGGTGAAATGAGGGGTTATTTGGGGG‑TTAGAGGAAAATT‑CCAGAAAAAAAAAACCCccgccg                                                                                                                                         <  2:611623/144‑1 (MQ=255)
                    aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                 >  2:341218/1‑151 (MQ=255)
                       cACGAAAAAAAGCCCTGCTTTATAGGAGTTATCGGCATGATTTTCCGTTTCGGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTTGGGGT‑TTACAGGAAAATC‑CCAGGCAGAAAAAAACCCCGCCGGTGGCgg                                                                                                                                  <  1:1496521/147‑1 (MQ=255)
                        aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                               >  1:336572/1‑149 (MQ=255)
                          gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                               >  2:1385560/1‑147 (MQ=255)
                           aTAAACAGCCCTGTGTTATATGAGTTATCGGCATGATTTTCCTTTTCTGCAGGTTTTTTAACCGGACGTTTGGTGGGTAAAATGGGGGTTTTTTTGGGGG‑TTAAAGTTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                          <  2:1101829/151‑1 (MQ=255)
                           aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc                                                                                                                          >  1:978156/1‑151 (MQ=255)
                             aaaCAGCCCTGCTTTATATGATTTACCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTGAAATGGGGGTTTTTTGGGGGG‑TTACAGGTAAATT‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                          <  1:1106381/149‑1 (MQ=255)
                                       gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                               >  1:1368979/1‑150 (MQ=255)
                                         gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                            >  2:428499/1‑151 (MQ=255)
                                          ttATATGATTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTTTTGGGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                            <  1:109474/150‑1 (MQ=255)
                                          ttATATGAGTTACCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGTTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                            <  1:129091/150‑1 (MQ=255)
                                             tatGAGTTATCGGCATGATTACCCGTTTCTGCAGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGAGGGGTTATTGGGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                          <  1:82726/149‑1 (MQ=255)
                                              atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                       >  2:1238759/1‑151 (MQ=255)
                                                  gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                   >  2:1014633/1‑151 (MQ=255)
                                                   ttATCGGCATTATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTGGGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                    <  1:644159/149‑1 (MQ=255)
                                                   ttAGCGGCATGATTACCGGTTTCGGCGGGGTTTTTAATCGGACGATTGGTGGGTAAAATGAGGGTTTATTGGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                        <  1:270725/145‑1 (MQ=255)
                                                       cgtcttgtttttCCTTTTCTGCGGGTTTTTTAATGGAACTTTTAGGGGTTAAAAGGGGGGTTTATTGGGGTT‑TTACAGGAAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                <  1:958946/138‑1 (MQ=255)
                                                        ggCATATTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGGGAAATGGGGGGTTATTTGGGGG‑TTACAGGAAAATC‑CCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                              <  1:846194/150‑1 (MQ=255)
                                                         gCATGATTTCCCGTTTCTGCAGGTTTTTTAATCGGAGGATTAGGGGGTAAAATGAGGGTTTATTGGGGGG‑TTACAGGTAAATC‑CCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                            <  1:359620/151‑1 (MQ=255)
                                                           aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                           >  2:387281/1‑150 (MQ=255)
                                                               ttATCCGTTTCGGCAGGTTTTTAAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTGGGGGT‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                        <  1:269271/149‑1 (MQ=255)
                                                               ttACCCTTTTCTGCAGGTTTTTAAATCGGACGTTTAGTGGGGAAAAGGGGGGTTTATTGGGGGG‑TTACAGGAAAATT‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                       <  1:1094035/150‑1 (MQ=255)
                                                                 aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                    >  2:371781/1‑151 (MQ=255)
                                                                  tCCTTTTCGGCAGGTTTTTTAAGCGGACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGG‑TAACAGGAAAATC‑CCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                    <  2:1249036/150‑1 (MQ=255)
                                                                   ccGTTTCTGCGGGTTTTTTAAGCGGACGATTGGGGGGTGAAATGGGGGTTTATTGGGGGG‑TTACAGATAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                     <  1:340381/148‑1 (MQ=255)
                                                                    cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCCTTATGTTGCCGGATGCGGCGt                                                                                  >  2:1268469/1‑150 (MQ=255)
                                                                      tttCTGGAGGGTTTTTAATCGGACGTTTAGTGGGTAAATTGGGGGTTTTTTTGGGGG‑TTACAGAAAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                   <  2:1475927/147‑1 (MQ=255)
                                                                      tttCTGCGGGTTTTTAAATCGAGCGATTAGGGGGAAAAATGAGGTTTTTTTTGGGGT‑TTAAAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                <  2:990340/150‑1 (MQ=255)
                                                                      tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                >  2:796919/1‑150 (MQ=255)
                                                                        tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                              >  2:408032/1‑150 (MQ=255)
                                                                           gcggggTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTTTTTGGGGT‑TTACAGTAAAATT‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                             <  2:1257190/145‑1 (MQ=255)
                                                                           gagggttttttAATGGGACGATTAGTGGGTAAAAGGGGGGTTTTTTGGGGGG‑TTACAGGTAAATC‑CCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                            <  1:1185186/143‑1 (MQ=255)
                                                                           gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                          >  2:1279577/1‑151 (MQ=255)
                                                                              ggttttttAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTTGGGGG‑TTACAGGTAAATC‑CCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                        <  1:921534/147‑1 (MQ=255)
                                                                                 tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                    >  1:902042/1‑151 (MQ=255)
                                                                                 tttttAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                       <  1:1318495/148‑1 (MQ=255)
                                                                                 tttttAAGCGGACTTTTAGGGGGTAAAAGGAGGGTTTTTTTGGGGT‑TCACAGGAAATTT‑CAAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                     <  1:999810/150‑1 (MQ=255)
                                                                                 ttttaAATCGGAGGATTAGTGGGTAAAATAGGGGGTTATTGGGGGT‑TGACGGATAATTC‑CCAGCAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTAAACGCCtt                                                                         <  1:1429689/146‑1 (MQ=255)
                                                                                   tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                  >  2:12800/1‑151 (MQ=255)
                                                                                   tttAATCGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGG‑TTACAGGAAAATT‑CCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                   <  1:773757/150‑1 (MQ=255)
                                                                                      aaTCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTTGGGGG‑TTACA‑GTAAATTCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                <  2:1525645/150‑1 (MQ=37)
                                                                                      aaTCGGACGATTAGGGGGTAAAAGGGGGGGTTTTTTGGGGG‑TTAAAGGTAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                               <  2:879909/151‑1 (MQ=38)
                                                                                       aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                              >  2:1191909/1‑151 (MQ=255)
                                                                                          ggACGATTAGTGGGTGAAATGGGGGTTTTTTTGGGGG‑TTACAGGTAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                           <  2:1152971/151‑1 (MQ=38)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                          >  1:208267/1‑150 (MQ=255)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                          >  1:868228/1‑150 (MQ=255)
                                                                                                ttggTGGGTAAAATGAGGGTTTTTTTGGGGG‑TTACAGGTAAATC‑CCAGGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                             <  2:617009/140‑1 (MQ=37)
                                                                                                   gggggTAAAAGGGGGGGTTTTTTGGGGG‑TTACAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                    <  2:270133/147‑1 (MQ=37)
                                                                                                     gggTAAAATGAGGGTTTATTGGGGGT‑TAACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                <  1:147637/151‑1 (MQ=37)
                                                                                                     gggAGAAAGGGGGGTTTTTTGGGGGG‑TTACAGGTAAATT‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                     <  2:946759/146‑1 (MQ=37)
                                                                                                      ggTAAAATGGGGGTTTATTTGGGGG‑TTACAGGAAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                               <  2:460501/151‑1 (MQ=37)
                                                                                                         gAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                             >  1:179285/1‑150 (MQ=255)
                                                                                                               aGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      >  1:1162648/1‑151 (MQ=255)
                                                                                                               aGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      >  1:830791/1‑151 (MQ=255)
                                                                                                                ggggTTTTTTGGGGG‑TTAAAGGTAAATC‑CCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      <  2:490769/150‑1 (MQ=37)
                                                                                                                ggggTTATTTGGGGG‑TTAGAGAAAATTT‑CCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      <  1:716452/150‑1 (MQ=37)
                                                                                                                 gggTTTTTTGGGGGTTTACAGAAAAATC‑CCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      <  1:41832/150‑1 (MQ=25)
                                                                                                                  ggTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                    >  2:344386/1‑150 (MQ=255)
                                                                                                                     tATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                >  1:1464152/1‑151 (MQ=255)
                                                                                                                       ttgggggt‑TTACAGGAAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                <  1:436657/146‑1 (MQ=37)
                                                                                                                        ttGGGGT‑TTACAGGAAACTC‑CCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              <  2:426647/150‑1 (MQ=25)
                                                                                                                        ttGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              >  1:812632/1‑150 (MQ=255)
                                                                                                                        ttGGGGG‑CTACAGGGAAATT‑CAAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGATGAGAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTCATCCGGCCTACAAAAATGAGCAAATTCAATAAATTGCAATTCa                              >  2:1211672/1‑150 (MQ=37)
                                                                                                                        tgggggg‑TTACAGGTAAATC‑CCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              <  1:169189/148‑1 (MQ=255)
                                                                                                                          ggggg‑TTACAGGAAAATC‑CCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              <  1:1501110/148‑1 (MQ=255)
                                                                                                                            ggg‑TTACAGGTAAATC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                         <  2:339809/151‑1 (MQ=255)
                                                                                                                                   cAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                   >  1:1227693/1‑151 (MQ=255)
                                                                                                                                        aaaTT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt               >  1:1309265/1‑150 (MQ=255)
                                                                                                                                        aaaTT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt               >  1:618107/1‑150 (MQ=255)
                                                                                                                                        aaaTT‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa              <  2:225767/151‑1 (MQ=255)
                                                                                                                                        aaaTC‑CCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                <  1:1190976/149‑1 (MQ=255)
                                                                                                                                         aaTT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa              >  2:217706/1‑150 (MQ=255)
                                                                                                                                          aTT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa              >  1:818634/1‑149 (MQ=255)
                                                                                                                                              ccAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc         >  1:1130009/1‑151 (MQ=255)
                                                                                                                                               cAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc         >  1:23696/1‑150 (MQ=255)
                                                                                                                                                aGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca        <  1:960007/150‑1 (MQ=255)
                                                                                                                                                 ggCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgca        <  2:1238258/149‑1 (MQ=255)
                                                                                                                                                 ggCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag       <  2:60393/150‑1 (MQ=255)
                                                                                                                                                     gAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt  >  1:781633/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAT  >  NC_000913/2001764‑2002061

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