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

CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAT  >  NC_000913/2001764‑2002061
                                                                                                                                                     |                                                                                                                                                      
ccTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaa                                                                                                                                                      >  2:251299/1‑150 (MQ=255)
  ttCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaa                                                                                                                                                    >  1:490788/1‑150 (MQ=255)
   tCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaa                                                                                                                                                   <  1:1281217/150‑1 (MQ=255)
    cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGaaaaaaa                                                                                                                                                 <  2:1167315/151‑1 (MQ=255)
     aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAc                                                                                                                                                >  1:225776/1‑151 (MQ=255)
      gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAcc                                                                                                                                               <  2:1481258/151‑1 (MQ=255)
       cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAAccc                                                                                                                                              >  1:548539/1‑151 (MQ=255)
         gTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccg                                                                                                                                            >  2:1186139/1‑151 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                           <  1:796517/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                           <  2:1285553/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                           <  2:122720/150‑1 (MQ=255)
           ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCccgc                                                                                                                                           <  1:135036/150‑1 (MQ=255)
            taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCGGGCAGAAAAAAACCccgcc                                                                                                                                          <  2:600435/150‑1 (MQ=255)
             agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCgg                                                                                                                                        <  1:314830/151‑1 (MQ=255)
                gtaaaTCCACGATAAACAGCCCTGTGTTTTAAGGGTTACGGGCAGGATTATCCGTTTCTGCAGGGTTTTTATTCGGACGATGAGGGGGTGAAATGAGGGGTTATTTGGGGGTTAGAGGAAAATTCCAGAAAAAAAAAACCCccgccg                                                                                                                                         <  2:611623/144‑1 (MQ=255)
                    aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                 >  2:341218/1‑151 (MQ=255)
                       cACGAAAAAAAGCCCTGCTTTATAGGAGTTATCGGCATGATTTTCCGTTTCGGCGGGGTTTTTAATCGGACGATTAGGGGGTAAAAGGGGGGTTTTTTTGGGGTTTACAGGAAAATCCCAGGCAGAAAAAAACCCCGCCGGTGGCgg                                                                                                                                  <  1:1496521/147‑1 (MQ=255)
                        aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                               >  1:336572/1‑149 (MQ=255)
                          gATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                               >  2:1385560/1‑147 (MQ=255)
                           aTAAACAGCCCTGTGTTATATGAGTTATCGGCATGATTTTCCTTTTCTGCAGGTTTTTTAACCGGACGTTTGGTGGGTAAAATGGGGGTTTTTTTGGGGGTTAAAGTTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                          <  2:1101829/151‑1 (MQ=255)
                           aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCAc                                                                                                                          >  1:978156/1‑151 (MQ=255)
                             aaaCAGCCCTGCTTTATATGATTTACCGGCATGTTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGGGGGTGAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                          <  1:1106381/149‑1 (MQ=255)
                                       gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                               >  1:1368979/1‑150 (MQ=255)
                                         gTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                            >  2:428499/1‑151 (MQ=255)
                                          ttATATGATTTATCGGCATGATTATCCGTTTCTGCAGGTTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                            <  1:109474/150‑1 (MQ=255)
                                          ttATATGAGTTACCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGTTTATTTGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAtt                                                                                                            <  1:129091/150‑1 (MQ=255)
                                             tatGAGTTATCGGCATGATTACCCGTTTCTGCAGGTTTTTTAATGGGACGATTAGGGGGTAAAAGGAGGGGTTATTGGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                          <  1:82726/149‑1 (MQ=255)
                                              atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                       >  2:1238759/1‑151 (MQ=255)
                                                  gTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTAt                                                                                                   >  2:1014633/1‑151 (MQ=255)
                                                   ttATCGGCATTATTATCCGTTTCTGCAGGGTTTTTAATCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTa                                                                                                    <  1:644159/149‑1 (MQ=255)
                                                   ttAGCGGCATGATTACCGGTTTCGGCGGGGTTTTTAATCGGACGATTGGTGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATTCCAGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgcgc                                                                                                        <  1:270725/145‑1 (MQ=255)
                                                       cgtcttgtttttCCTTTTCTGCGGGTTTTTTAATGGAACTTTTAGGGGTTAAAAGGGGGGTTTATTGGGGTTTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                <  1:958946/138‑1 (MQ=255)
                                                        ggCATATTTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGGGAAATGGGGGGTTATTTGGGGGTTACAGGAAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGc                                                                                              <  1:846194/150‑1 (MQ=255)
                                                         gCATGATTTCCCGTTTCTGCAGGTTTTTTAATCGGAGGATTAGGGGGTAAAATGAGGGTTTATTGGGGGGTTACAGGTAAATCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                            <  1:359620/151‑1 (MQ=255)
                                                           aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                           >  2:387281/1‑150 (MQ=255)
                                                               ttATCCGTTTCGGCAGGTTTTTAAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTGGGGGTTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATg                                                                                        <  1:269271/149‑1 (MQ=255)
                                                               ttACCCTTTTCTGCAGGTTTTTAAATCGGACGTTTAGTGGGGAAAAGGGGGGTTTATTGGGGGGTTACAGGAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                       <  1:1094035/150‑1 (MQ=255)
                                                                 aTCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                    >  2:371781/1‑151 (MQ=255)
                                                                  tCCTTTTCGGCAGGTTTTTTAAGCGGACGATTGGGGGGTAAAAGGAGGGTTTTTTTGGGGGTAACAGGAAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggc                                                                                    <  2:1249036/150‑1 (MQ=255)
                                                                   ccGTTTCTGCGGGTTTTTTAAGCGGACGATTGGGGGGTGAAATGGGGGTTTATTGGGGGGTTACAGATAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcgg                                                                                     <  1:340381/148‑1 (MQ=255)
                                                                    cGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCCTTATGTTGCCGGATGCGGCGt                                                                                  >  2:1268469/1‑150 (MQ=255)
                                                                      tttCTGGAGGGTTTTTAATCGGACGTTTAGTGGGTAAATTGGGGGTTTTTTTGGGGGTTACAGAAAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcggcg                                                                                   <  2:1475927/147‑1 (MQ=255)
                                                                      tttCTGCGGGTTTTTAAATCGAGCGATTAGGGGGAAAAATGAGGTTTTTTTTGGGGTTTAAAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                <  2:990340/150‑1 (MQ=255)
                                                                      tttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTaa                                                                                >  2:796919/1‑150 (MQ=255)
                                                                        tCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                              >  2:408032/1‑150 (MQ=255)
                                                                           gcggggTTTTTAATCGGACGATTAGTGGGTAAAATGGGGGGTTTTTTGGGGTTTACAGTAAAATTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACg                                                                             <  2:1257190/145‑1 (MQ=255)
                                                                           gagggttttttAATGGGACGATTAGTGGGTAAAAGGGGGGTTTTTTGGGGGGTTACAGGTAAATCCCGGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGc                                                                            <  1:1185186/143‑1 (MQ=255)
                                                                           gCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCt                                                                          >  2:1279577/1‑151 (MQ=255)
                                                                              ggttttttAATCGGACGATTAGTGGGTAAAAGGAGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                        <  1:921534/147‑1 (MQ=255)
                                                                                 tttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                    >  1:902042/1‑151 (MQ=255)
                                                                                 tttttAATCGGACGATTAGTGGGTAAAATGGGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTAt                                                                       <  1:1318495/148‑1 (MQ=255)
                                                                                 tttttAAGCGGACTTTTAGGGGGTAAAAGGAGGGTTTTTTTGGGGTTCACAGGAAATTTCAAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                     <  1:999810/150‑1 (MQ=255)
                                                                                 ttttaAATCGGAGGATTAGTGGGTAAAATAGGGGGTTATTGGGGGTTGACGGATAATTCCCAGCAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTAAACGCCtt                                                                         <  1:1429689/146‑1 (MQ=255)
                                                                                   tttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGc                                                                  >  2:12800/1‑151 (MQ=255)
                                                                                   tttAATCGGACGATTAGTGGGTAAAATGAGGGGTTATTTGGGGGTTACAGGAAAATTCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                   <  1:773757/150‑1 (MQ=255)
                                                                                      aaTCGGACGTTTAGTGGGTAAAATGGGGGTTTTTTTGGGGGTTACAGTAAATTCCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                <  2:1525645/150‑1 (MQ=37)
                                                                                      aaTCGGACGATTAGGGGGTAAAAGGGGGGGTTTTTTGGGGGTTAAAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                               <  2:879909/151‑1 (MQ=38)
                                                                                       aTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTAc                                                              >  2:1191909/1‑151 (MQ=255)
                                                                                          ggACGATTAGTGGGTGAAATGGGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaa                                                           <  2:1152971/151‑1 (MQ=38)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                          >  1:208267/1‑150 (MQ=255)
                                                                                            aCGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaa                                                          >  1:868228/1‑150 (MQ=255)
                                                                                                ttggTGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGGTAAATCCCAGGCAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                             <  2:617009/140‑1 (MQ=37)
                                                                                                   gggggTAAAAGGGGGGGTTTTTTGGGGGTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                    <  2:270133/147‑1 (MQ=37)
                                                                                                     gggTAAAATGAGGGTTTATTGGGGGTTAACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                <  1:147637/151‑1 (MQ=37)
                                                                                                     gggAGAAAGGGGGGTTTTTTGGGGGGTTACAGGTAAATTCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                     <  2:946759/146‑1 (MQ=37)
                                                                                                      ggTAAAATGGGGGTTTATTTGGGGGTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAtt                                               <  2:460501/151‑1 (MQ=37)
                                                                                                         gAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                             >  1:179285/1‑150 (MQ=255)
                                                                                                               aGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      >  1:1162648/1‑151 (MQ=255)
                                                                                                               aGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      >  1:830791/1‑151 (MQ=255)
                                                                                                                ggggTTTTTTGGGGGTTAAAGGTAAATCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      <  2:490769/150‑1 (MQ=37)
                                                                                                                ggggTTATTTGGGGGTTAGAGAAAATTTCCGGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      <  1:716452/150‑1 (MQ=37)
                                                                                                                 gggTTTTTTGGGGGTTTACAGAAAAATCCCGGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                      <  1:41832/150‑1 (MQ=25)
                                                                                                                  ggTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                    >  2:344386/1‑150 (MQ=255)
                                                                                                                     tATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                >  1:1464152/1‑151 (MQ=255)
                                                                                                                       ttgggggtTTACAGGAAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                <  1:436657/146‑1 (MQ=37)
                                                                                                                        ttGGGGTTTACAGGAAACTCCCAGAAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              <  2:426647/150‑1 (MQ=25)
                                                                                                                        ttGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              >  1:812632/1‑150 (MQ=255)
                                                                                                                        ttGGGGGCTACAGGGAAATTCAAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGATGAGAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTCATCCGGCCTACAAAAATGAGCAAATTCAATAAATTGCAATTCa                              >  2:1211672/1‑150 (MQ=37)
                                                                                                                        tggggggTTACAGGTAAATCCCAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              <  1:169189/148‑1 (MQ=255)
                                                                                                                          gggggTTACAGGAAAATCCCAGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                              <  1:1501110/148‑1 (MQ=255)
                                                                                                                            gggTTACAGGTAAATCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTg                         <  2:339809/151‑1 (MQ=255)
                                                                                                                                   cAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGcc                   >  1:1227693/1‑151 (MQ=255)
                                                                                                                                        aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt               >  1:1309265/1‑150 (MQ=255)
                                                                                                                                        aaaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt               >  1:618107/1‑150 (MQ=255)
                                                                                                                                        aaaTTCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa              <  2:225767/151‑1 (MQ=255)
                                                                                                                                        aaaTCCCAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                <  1:1190976/149‑1 (MQ=255)
                                                                                                                                         aaTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa              >  2:217706/1‑150 (MQ=255)
                                                                                                                                          aTTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa              >  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)
                                                                                                                                                     |                                                                                                                                                      
CCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAT  >  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: