Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,004,648	0	T	A	40.8%	113.1 / 42.0	72	N234K (AAT→AAA)	ycbF	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (6/36);  new base A (29/0);  total (36/36)
Fisher's exact test for biased strand distribution p-value = 2.37e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCC  >  NC_000913/1004504‑1004795                                                                                                                                                 |                                                                                                                                                              aaaTGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGGGTGGGAAAAAGGAGCAAAAAA‑TG‑ccc                                                                                                                                                       >  1:505270/1‑151 (MQ=255) aaaTGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGTGGGAAGAAAGAGGAAAAAA‑TG‑ccc                                                                                                                                                       >  2:33657/1‑151 (MQ=255)   aTGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGGGGAAAGAAAGAGGAAAAAA‑TG‑ccccg                                                                                                                                                     >  1:239559/1‑149 (MQ=255)        aaaGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGTGAGGAAAGAAAGAGCAAAAAC‑TT‑CCCT‑G‑aaaa                                                                                                                                               >  2:264241/1‑149 (MQ=255)        aaaGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGGGGAAAAAAAGAGGAAAAAC‑TT‑cccctaaaa                                                                                                                                                 >  2:146804/1‑144 (MQ=255)          aGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑Agg                                                                                                                                           <  2:511973/151‑1 (MQ=255)            cGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGGGGAAAAAAAGAGGAAAAAC‑TT‑CCCT‑G‑aaaaaaggt                                                                                                                                          >  1:1109019/1‑145 (MQ=255)             gTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGTGGAAAGAAAGAGGAAAAAA‑TT‑CCCC‑G‑aaaaaaggt                                                                                                                                          >  2:103699/1‑144 (MQ=255)               tAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGGAAAGAAGGAGGAAAAAA‑TG‑CCCC‑G‑AAAA‑AAGGTCgg                                                                                                                                      >  2:951095/1‑151 (MQ=255)                aaGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGa                                                                                                                                     <  1:118532/151‑1 (MQ=255)                 aGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGaa                                                                                                                                    <  2:986644/151‑1 (MQ=255)                  gTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGaaa                                                                                                                                   <  1:1150056/151‑1 (MQ=255)                  gTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGaaa                                                                                                                                   <  1:923664/151‑1 (MQ=255)                  gTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGaaa                                                                                                                                   <  1:1080714/151‑1 (MQ=255)                   tAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGTGGAAAGAAAGAGGAAAAAA‑TG‑CCCT‑G‑AAAA‑AAGGTCGGgaa                                                                                                                                   >  2:22023/1‑150 (MQ=255)                    aaTGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGGGGGGGAAAAAAGAGGAAAAAC‑TG‑CCCC‑C‑CAAC‑AAGGTCGGgaaa                                                                                                                                  >  1:1209019/1‑150 (MQ=255)                    aaTGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTTAGGGGAAAGAAAGAGCAAAAAC‑TT‑CCCT‑G‑AAAA‑AAGGTCGGgaaa                                                                                                                                  >  2:223543/1‑150 (MQ=255)                     aTGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGGGTGGAAAGAAAGAGGAAAAAA‑TA‑CCCC‑G‑AAAA‑AAGTTCGGgaaa                                                                                                                                  >  1:1224534/1‑149 (MQ=255)                     aTGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑g                                                                                                                               <  2:1313829/151‑1 (MQ=255)                         nggAGTTGTGCCGCCTAAAAGCCAGCGTCNGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATNANNANGGCGGGTTGGGGGGGAAAAAGGANNNAAACCTTC‑CCCC‑A‑AAAA‑AGGTGGGGAAAAG‑Gttt                                                                                                                            >  1:1042871/2‑151 (MQ=255)                           gAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑g                                                                                                                        <  1:820894/151‑1 (MQ=255)                            aGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGCAAAGAAAGGGCAAAAAA‑TG‑CCCTCA‑AAAA‑AGGTTGGGAAACG‑Gtttcgg                                                                                                                         >  1:1335566/1‑148 (MQ=255)                             gTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑Gtt                                                                                                                      <  1:671450/151‑1 (MQ=255)                               tgtgCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGTGTGCTAAAAAGGAGAAAAAAA‑TGTCCCG‑A‑AAAA‑AGGTGGGGAAAAG‑Gttttgtttt                                                                                                                      >  2:552956/1‑145 (MQ=255)                               tgtgCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGTGGGCAAAAAAAGAGAAAAAAA‑TG‑CCCC‑GCAAAA‑AGGTTGGGAAAAG‑GTTTCG‑Gttt                                                                                                                     >  1:1202050/1‑151 (MQ=255)                                gtgCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑Gttttt                                                                                                                   <  1:7852/151‑1 (MQ=255)                                gtgCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑Gtttt                                                                                                                    <  1:646337/150‑1 (MQ=255)                                    cgccTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGTGGAAAGAAAGAGCAAAAAC‑TG‑CCCC‑G‑GAAC‑AGGGTCGGGAAAA‑GGTTTC‑Gttttttttt                                                                                                               >  1:1167181/1‑150 (MQ=255)                                         aaaaGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGGGGGGAAAGAAGGAGGAAAAAC‑TT‑CCCC‑G‑AAAA‑AAGGTCGGGAAACGGTTTGG‑Gtttttttttcgcc                                                                                                           >  2:38047/1‑146 (MQ=255)                                         aaaaGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGTGGAAAGAAAGAGGAAAAAA‑AG‑CCCC‑GAAAAA‑AGGTCGGGAAAAG‑GTTGGG‑Gtttttttttcgcc                                                                                                           >  1:1142096/1‑146 (MQ=255)                                                     cAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGCAAAGAAAGAGAAAAAAT‑TG‑CCCC‑G‑AAAC‑AGGGTCGGGAAAA‑AATTGC‑GGTTTTTTTTTCGCC‑GGGCAAATTccc                                                                                             >  2:647993/1‑151 (MQ=255)                                                      aGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGTGTGGAAAGAAAGAGCAAAAAA‑TT‑CCCC‑C‑AAAA‑AAGATCGGGAAAA‑GCTTTT‑GTTTTTTTTTTCGCC‑GGGCAAATTCCCt                                                                                            >  1:1179967/1‑151 (MQ=255)                                                             cTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGGGCAAAGAAAGAGAAAAAAC‑TG‑CCCT‑G‑AAAA‑AAGATCGGGAAAA‑GGTTGC‑GGTTTTTTTTCCGCA‑GGTC‑AAAACCCTTTTTGTcc                                                                                    >  2:1299739/1‑151 (MQ=255)                                                               ggTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTcc                                                                                    <  1:387426/149‑1 (MQ=255)                                                               ggTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCg                                                                                   <  1:59716/150‑1 (MQ=255)                                                                  gtCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCCAAGAAAGAGCAAAAAC‑TT‑CCCC‑A‑CAAA‑AGGTTTGGGAAAA‑CTTTGA‑GGTTTTTTTTTCGCAGGGTA‑AATTCCTTTTTGTTCGGact                                                                                >  2:1290836/1‑149 (MQ=255)                                                                    cAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCg                                                                             <  1:130494/151‑1 (MQ=255)                                                                      ggCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGca                                                                           <  1:1125994/151‑1 (MQ=255)                                                                      ggCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTTGGTGGGAAAAAGGAGCAAAAAA‑TA‑ACCC‑G‑AAAA‑AAGGTCGGGAAAA‑GTTTTCGGTTTTTTTTTCGCCG‑GGGA‑AATTCCCTTTTGTTCCGactcccc                                                                            >  1:1314055/1‑145 (MQ=255)                                                                       gCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGGGTGGAAAAAAAGAGGAAAAAA‑TT‑CCCC‑G‑AAAA‑AAGGTTGGGAAAC‑GTTTTC‑GGTTTTTTTTTCGCC‑GGGC‑AAATTCCTTTTTGTCCggaatccccc                                                                          >  2:703362/1‑142 (MQ=255)                                                                             gCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACataata                                                                    <  1:797472/151‑1 (MQ=255)                                                                              cACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATAc                                                                   <  1:733013/151‑1 (MQ=255)                                                                               aCCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTTGGTGGAAAAAAAGAGCCAAAAC‑TG‑CCCC‑G‑AAAA‑AAGGTCGGAGAAA‑GGTTGT‑GGTTTTTTTTCCGCC‑GGGG‑AAATCCCCTTTTTTTCGGACACGCCCaaaatatc                                                                  >  1:339487/1‑149 (MQ=255)                                                                                    gtCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTTAGTGGAAAGAAGGAGCAAAAAC‑TT‑CCCT‑G‑AAAC‑AAGGTCGGGAAAA‑GTTTGC‑GTTTTTTTTTTCGCC‑GGGCAAATTCCCTTTT‑GTTCGAACTCCCCCCTAATACG‑ccccc                                                            >  2:912888/1‑151 (MQ=255)                                                                                    gtCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑ccccc                                                            <  2:568280/151‑1 (MQ=255)                                                                                    gtCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑ccccc                                                            <  1:409626/151‑1 (MQ=255)                                                                                      cATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGCCAAGAAGGGGGAAAAAC‑TT‑CCCC‑G‑CAAC‑AAGGTCGGGAAAA‑GGTTGT‑GGTTTTTTTTTCGCG‑GGGCAAATACCCTTTTTG‑CCGGACCCCCACATAATACA‑ccccccg                                                          >  1:37517/1‑149 (MQ=255)                                                                                            gCCTGGCGGGCGATAAATGATTACGGCGGGGTGTGTGCAAAAAAAGAGCAAAATC‑TG‑CCCT‑G‑AAAA‑AAGGTTGGGAAAC‑GGTTGC‑GTTTTTTTTTTCGCA‑GGGC‑AAATTCCTTTTTGTGCGGACTCGCACATAAAAAG‑CCCCCCGGTTgtc                                                    >  1:704905/1‑149 (MQ=255)                                                                                              cTGGCGGGCGATAAATGATTACGGCGGGGTGAGGGCCAAAAAGGAGCAAAAAA‑TG‑CCCC‑G‑AAACAAGGTTGGGAAAA‑‑GGTTGC‑GTTTTTTTTTTCGCA‑GGGC‑AAATACCTTTTTGTTCGAACTCGCACATAAAACC‑CCCCCCGGTTTcacc                                                  >  1:1155153/1‑150 (MQ=255)                                                                                                      cgATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAAt                                          <  2:95822/151‑1 (MQ=255)                                                                                                          aaaTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCa                                       <  1:951090/150‑1 (MQ=255)                                                                                                          aaaTGATTACGGCGGGTTGAGGGCAAAAAAAGGGGAAAAAC‑TG‑CCCC‑G‑AAAA‑AGGGTCGGGAAAA‑GGTTGC‑GTTTTTTTTTTCGCA‑GGGCAAATTCCCTTTT‑GTCCGGAATCGCACATAAAAAGCCCCCCGGTTTG‑ACACCCGGGATcc                                        >  1:1318300/1‑149 (MQ=255)                                                                                                           aaTGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAgg                                     <  2:703454/151‑1 (MQ=255)                                                                                                                  aCGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGtt                               <  2:185485/150‑1 (MQ=255)                                                                                                                    ggcggGTTGAGTGCAAAGAAGGAACAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCAt                            <  2:1119416/151‑1 (MQ=255)                                                                                                                      cggGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGt                          <  2:1204105/151‑1 (MQ=255)                                                                                                                       gggTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGta                         <  2:1092220/151‑1 (MQ=255)                                                                                                                         gTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtact                       >  1:828865/1‑151 (MQ=255)                                                                                                                         gTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtact                       <  2:1137720/151‑1 (MQ=255)                                                                                                                            gAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGtactac                     <  1:624868/150‑1 (MQ=255)                                                                                                                            gAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcc                    <  2:239552/151‑1 (MQ=255)                                                                                                                            gAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcc                    >  1:970657/1‑151 (MQ=255)                                                                                                                             aGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑G‑TTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTAcccc                  >  2:299949/1‑151 (MQ=255)                                                                                                                                gCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCt                 <  1:472588/150‑1 (MQ=255)                                                                                                                                  aaaGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcg              <  2:293603/151‑1 (MQ=255)                                                                                                                                  aaaGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcg              <  1:22024/151‑1 (MQ=255)                                                                                                                                     gaagGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttc           <  1:146811/151‑1 (MQ=255)                                                                                                                                     gaagGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCttcgttc           <  1:33659/151‑1 (MQ=255)                                                                                                                                       agGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTtcgttcgt         <  1:707621/151‑1 (MQ=255)                                                                                                                                          aGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTaaa      <  1:264247/151‑1 (MQ=255)                                                                                                                                            cAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGc    >  1:73543/1‑151 (MQ=255)                                                                                                                                              aaaTC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGccc  >  2:1183006/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                              AAATGGTAAAGCGTTAAGTAATGTGGGAGTTGTGCCGCCTAAAAGCCAGCGTCAGACAAGCTGGTGTCAGGCAATCGCACCGTGTCATGTCGCCTGGCGGGCGATAAATGATTACGGCGGGTTGAGTGCAAAGAAGGAGCAAAATC‑TG‑CCCT‑G‑AAAC‑AGGTTCGGAAAAC‑GTTTGC‑GTTTTTTTTGCCGCA‑GGTC‑AATTCCCTTTTGGTCCGAACTCGCACATAATACG‑CCCCCGGTTTGCACACCGGGAATCCAGGAGAGTTCATGTACTACCCCTTCGTTCGTAAAGCCC  >  NC_000913/1004504‑1004795

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