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

AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065
                                                                                                                                                         |                                                                                                                                                   
aGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAAc                                                                                                                                                    >  2:320120/1‑151 (MQ=255)
 gCGGTATAGAGTGAGTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAAcc                                                                                                                                                   <  2:83733/151‑1 (MQ=255)
 gCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAAcc                                                                                                                                                   >  1:242661/1‑151 (MQ=255)
  cGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAAccc                                                                                                                                                  >  2:454038/1‑151 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACTGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAAcccc                                                                                                                                                 >  2:33441/1‑151 (MQ=255)
   ggTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAAccc                                                                                                                                                  >  1:164350/1‑150 (MQ=255)
      ataGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCccgcc                                                                                                                                              >  2:491639/1‑151 (MQ=255)
        agagTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCgg                                                                                                                                            <  3:427047/151‑1 (MQ=255)
          agTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTg                                                                                                                                          >  1:636728/1‑151 (MQ=255)
            tGATTTCACGATAAACAGCCCTGCTTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGGGGGTGAAATGAGGGT‑TTTTTG‑GGGGG‑TTACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTga                                                                                                                                         <  2:8725/150‑2 (MQ=255)
              aaTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCgg                                                                                                                                      >  4:54062/1‑151 (MQ=255)
               aTTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCggg                                                                                                                                     >  3:326397/1‑151 (MQ=255)
                ttCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                    >  3:121823/1‑151 (MQ=255)
                 tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                    >  4:117273/1‑150 (MQ=255)
                 tCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCgggg                                                                                                                                    >  2:462977/1‑150 (MQ=255)
                  cACGATAAACAGCCCTGTTTTATATGATTTATGGGAATGATTATCGTTTTCTGCGGGTTTTTTAAGCGGACGTTTGGGGGGTAAAATGGGGGG‑TTTTTG‑GGGGG‑TTACAGGTAAATCCC‑AGGAAAAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   <  3:393839/150‑1 (MQ=255)
                  cACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  2:1221409/1‑150 (MQ=255)
                   aCGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGa                                                                                                                                   >  2:734683/1‑149 (MQ=255)
                    cGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGc                                                                                                                                >  1:159584/1‑151 (MQ=255)
                      aTAAACAGCCCTGCTTTATATGAGTTATGGGCTTGATTATCGGTTTCTGCAGGGTTTTAAATGGGACGATTAGTGGGTAAAAGGGGGGG‑TTTTTG‑GGGGG‑TTACAGGAAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  2:488676/151‑1 (MQ=255)
                       tAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACg                                                                                                                             >  1:864472/1‑151 (MQ=255)
                        aacccgCCCTGCTTTATAGGAGTTACGGGCATATTTATCCTTTTCTGAGGGGTTTTAAACCGGACATTTGGTGGGT‑AAATGGGGGGTTTTTTT‑GGGGT‑TTACAGGAAATTTCC‑GGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  2:1269292/144‑1 (MQ=255)
                        aaaCAGCCCTGCTTTATATGAGTTATGGGCATGATTATCCTTTTCTGCAGGTTTTTTAAGCGGACGATTAGGGGGTGAAATGGGGGT‑TTATTG‑GGGGG‑TTACAGGAAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACg                                                                                                                             <  1:790020/150‑1 (MQ=255)
                        aaaCAGCCCTGCGTTATAGGAGTTATGGGCATGATTATCCGTTTCTGCGGGGTTTTTAACCGGAAGATTAGTGGGTAAAGGGGGGGG‑TTTTTG‑GGGGT‑TTACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAg                                                                                                                                 <  2:289166/146‑1 (MQ=255)
                        aaaCACCCCTGCTTTATAGGGGTTTTCGGAATGTTTACCCGTTCCGGCGGGTTTTTTAATCGGACGATT‑GGGGGGAAAAAGGGGGG‑TTTTTT‑GGGGGTTTAAAGGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCAc                                                                                                                              <  1:1105434/149‑1 (MQ=255)
                         accAGCCCTGCGTTATATGAGTTATCGGCATTATTACCCTTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGG‑TTTTTG‑GGGGG‑TTACAGTTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                           <  2:1145058/149‑1 (MQ=255)
                          aCAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGtt                                                                                                                           >  4:163361/1‑150 (MQ=255)
                              cccTGCGTTATATGAGTTATCGGCAGGATTATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTAAAAGGAGGGG‑TTTTTG‑GGGGG‑TTAAGGGTAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTg                                                                                                                          <  3:279239/147‑1 (MQ=255)
                                  gCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACaa                                                                                                                   >  1:379289/1‑150 (MQ=255)
                                   cGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAAt                                                                                                                 >  3:414675/1‑151 (MQ=255)
                                        tatGATTTATCGGCATGATTATCCTTTTCTGCAGGTTTTTTAATGGGACGATTGGGGGGTAAAATGGGGGT‑TTTTTG‑GGGGG‑TTACAGGAAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTgc                                                                                                              <  2:851381/149‑1 (MQ=255)
                                         atGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCt                                                                                                           >  4:46006/1‑151 (MQ=255)
                                          tGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCtt                                                                                                          >  2:949535/1‑151 (MQ=255)
                                               tATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGt                                                                                                     >  2:1170014/1‑151 (MQ=255)
                                                   ggcagaaTTATCCGTTTCTGAAGGGTTTTAAATCGGACGATTAGGGGGTAAAATGGGGGG‑TTTTTT‑GGGGG‑TTACAGGAAAATCCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGtt                                                                                                    <  4:432102/142‑1 (MQ=255)
                                                   ggCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                 >  4:297048/1‑151 (MQ=255)
                                                    gCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTAAAATGAGGGT‑TTTTTG‑GGGGG‑TTACAGGTAATTTCC‑AGGAAAAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCg                                                                                                <  1:813858/151‑1 (MQ=255)
                                                      aTGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                               >  3:142984/1‑150 (MQ=255)
                                                         aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                           >  2:464221/1‑151 (MQ=255)
                                                         aTTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgc                                                                                           >  2:126172/1‑151 (MQ=255)
                                                           tATCCTTTTCTGTAGTTTTTTAAAACGGATATTTAGTGGGGAAAAGGAGGGT‑TTATTT‑GGGTG‑TTACAGTAAAATCCC‑AGGAAAAAAAACCCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCgg                                                                                               <  1:130995/145‑1 (MQ=255)
                                                           tATCCGTTTCTGCGGGTTTTTTAATCGGACGATTAGTGGGTGAAATGGGGGG‑TTTTTG‑GGGGG‑TTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATgcg                                                                                          <  4:130378/150‑1 (MQ=255)
                                                              ccGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGt                                                                                      >  2:267696/1‑151 (MQ=255)
                                                                  ttCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAAc                                                                                  >  4:191292/1‑151 (MQ=255)
                                                                         gggTTTTTAATGGGACGATTAGGGGGTAAAATGGGGGT‑TTTTTG‑GGGGG‑TTACAGGTAAATCCC‑GGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTa                                                                            <  1:291057/150‑1 (MQ=255)
                                                                            tttttAATCGGAGGATTAGTGGGTGAAATGGGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATCCC‑GGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                         <  2:943279/150‑1 (MQ=255)
                                                                            tttttAATCGGACGTTTAGTGGGGAAAATGGGGGG‑TTATTG‑GGGGT‑TTACAGG‑AAATTCCAAGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATcc                                                                         <  2:950897/150‑1 (MQ=255)
                                                                            tttttAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCg                                                                        >  3:401032/1‑151 (MQ=255)
                                                                             ttttAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCgg                                                                       >  1:436509/1‑151 (MQ=255)
                                                                                 aaTCGGACGATTAGGGGGTGAAATGAGGGG‑TTTTTT‑GGGGG‑TTACAGGTAAATCCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    <  1:217561/150‑1 (MQ=255)
                                                                                 aaTCGGACAATTAGGGGGTAAAATGAGGGG‑TGATTG‑GGGTG‑TTACGGGAAAATCCC‑GGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCt                                                                    <  3:90963/150‑1 (MQ=37)
                                                                                 aaTCGGAAGATGAGTGGGGAAAAGGGGGGT‑TTATTG‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTa                                                                   <  2:204969/151‑1 (MQ=38)
                                                                                     ggACGATTAGTGGGTAAAATGAGGGT‑TTATTG‑GGGGG‑TTAGGGGTAAATTCC‑GGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                 <  2:501475/149‑1 (MQ=37)
                                                                                     ggACGATTAGGGGGTAAAATGAGGGT‑TTTTTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACa                                                                 <  1:370971/149‑1 (MQ=38)
                                                                                         gATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtg                                                           >  1:30126/1‑151 (MQ=255)
                                                                                            tAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGc                                                        >  1:929075/1‑151 (MQ=255)
                                                                                             agggggTGAAAGGGGGGT‑TTATTG‑GGGGG‑TTAAAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                         <  2:1028407/146‑1 (MQ=37)
                                                                                              gggggTGAAATGGGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATCCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCa                                                       <  3:197569/148‑1 (MQ=255)
                                                                                                gggTGAAATGGGGGT‑TTATTT‑GGGGT‑TTACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAAt                                                    <  1:414842/151‑1 (MQ=38)
                                                                                                gggTGAAATGAGGGT‑TTTTTG‑GGGGG‑TTAAGGGTAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAtgtg                                                         <  2:689527/146‑1 (MQ=37)
                                                                                                    aaaaTGA‑GGG‑TTATTTGGGGGG‑TCAAAGATAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                 <  4:277709/149‑1 (MQ=25)
                                                                                                     aaaggggggt‑TTATTG‑GGGGT‑TTAAAGATAAATCCC‑GGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               <  2:938124/145‑1 (MQ=25)
                                                                                                     aaaTGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               >  3:356800/1‑151 (MQ=255)
                                                                                                     aaaTGAGGGG‑TTATTG‑GGGGG‑TAACAGGTAAATCCC‑TGCCAAAAAAAATCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCa                                                 <  1:199212/149‑1 (MQ=37)
                                                                                                      aaatggggg‑TTTTTT‑GGGGG‑GTAAAAGAAAAATCCAAGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCaaa                                                     <  2:675487/140‑1 (MQ=25)
                                                                                                      aaTGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAAt                                               >  2:813317/1‑150 (MQ=255)
                                                                                                        tGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                             >  1:1053237/1‑150 (MQ=255)
                                                                                                         gAGGGT‑TTTTTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAt                                           <  4:21750/151‑1 (MQ=255)
                                                                                                         gAGGGT‑TTATTT‑GGGGG‑TTACAGGTAATTTCC‑AGGAAGAAAAAACCCCCGCCGGTGGCGGGAAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                             <  3:290190/149‑1 (MQ=255)
                                                                                                           gggt‑TTATTT‑GGGGT‑TTACAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  2:10595/150‑1 (MQ=255)
                                                                                                           gggg‑TTATTT‑GGGGG‑TTACGGATAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                          <  4:132641/150‑1 (MQ=255)
                                                                                                            ggg‑TTTTTT‑GGGGGTTTACAGGTAATTCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATaa                                             <  2:636604/147‑1 (MQ=37)
                                                                                                            ggg‑TTATTG‑GGGGG‑TTACAGTTAACTCCC‑AGACAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGc                                        <  1:1200333/151‑1 (MQ=37)
                                                                                                                 tATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    >  1:987462/1‑151 (MQ=255)
                                                                                                                 tATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGCTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    >  1:60938/1‑151 (MQ=255)
                                                                                                                   ttt‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTc                                   >  2:1008914/1‑150 (MQ=255)
                                                                                                                    ttgggggt‑TTACAGGTAAATTCC‑GGGAAGAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAAtt                                    <  1:879595/146‑1 (MQ=255)
                                                                                                                     tgggggt‑TTAAAGGTAAATCCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCa                                  <  1:364829/148‑1 (MQ=37)
                                                                                                                       ggggg‑TTACAGGTAAATTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAAc                                <  1:399847/150‑1 (MQ=255)
                                                                                                                        gggg‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                              >  1:1139030/1‑151 (MQ=255)
                                                                                                                           gtttACAGGTAAATCCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            <  1:320172/149‑1 (MQ=255)
                                                                                                                           g‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGt                            >  1:317209/1‑150 (MQ=255)
                                                                                                                              tACGGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                          <  1:525409/150‑1 (MQ=255)
                                                                                                                              tACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                         >  4:302268/1‑151 (MQ=255)
                                                                                                                              tACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAgg                         >  2:474633/1‑151 (MQ=255)
                                                                                                                               aCAGGTAAATTCC‑AGGCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGc                        <  2:615173/151‑1 (MQ=255)
                                                                                                                                   gTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTg                     >  1:1128020/1‑150 (MQ=255)
                                                                                                                                     aaaTTCC‑AGGAAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  <  2:676697/151‑1 (MQ=255)
                                                                                                                                      aaTTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  1:987663/1‑150 (MQ=255)
                                                                                                                                       aTTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATa                  >  4:264772/1‑149 (MQ=255)
                                                                                                                                        ttCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagc               >  4:456445/1‑151 (MQ=255)
                                                                                                                                           c‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  1:1273113/1‑150 (MQ=255)
                                                                                                                                           c‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  4:363925/1‑150 (MQ=255)
                                                                                                                                             aGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgc             >  1:1137376/1‑149 (MQ=255)
                                                                                                                                              ggCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcagc          >  4:360711/1‑151 (MQ=255)
                                                                                                                                              ggCAAAAAAAACCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAagcgcag           <  1:974614/150‑1 (MQ=255)
                                                                                                                                                      aaacCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  <  1:523129/151‑1 (MQ=255)
                                                                                                                                                      aaaaCCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAgg  >  2:418192/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
AGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGG‑TTATTT‑GGGGG‑TTACAGGTAAATTCC‑AGGCAGAAAAAAACCCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGG  >  NC_000913/2001769‑2002065

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