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

CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001768‑2002065
                                                                                                                                                        |                                                                                                                                                   
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                                                                                                                                               <  2:122720/150‑1 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                               <  1:796517/150‑1 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                               <  1:135036/150‑1 (MQ=255)
       ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCccgc                                                                                                                                               <  2:1285553/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:830791/1‑151 (MQ=255)
                                                                                                           aGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          >  1:1162648/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:618107/1‑150 (MQ=255)
                                                                                                                                    aaaTT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAt                   >  1:1309265/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)
                                                                                                                                                  aaaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCAt      <  1:383131/150‑1 (MQ=255)
                                                                                                                                                   aaaaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCa    <  2:387676/151‑1 (MQ=255)
                                                                                                                                                     aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  2:127077/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                   
CAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGG‑TTACAGGTAAATT‑CCAGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001768‑2002065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑