Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	42.3%	177.8 / 33.2	97	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (49/7);  new base C (0/41);  total (49/48)
Fisher's exact test for biased strand distribution p-value = 1.82e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTC  >  NC_000913/993764‑994059                                                                                                                                                           |                                                                                                                                                   cACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGc                                                                                                                                                   >  2:1244530/1‑151 (MQ=255)  aCGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGcc                                                                                                                                                  >  2:1001137/1‑151 (MQ=255)   cGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGccc                                                                                                                                                 >  2:886226/1‑151 (MQ=255)    ggTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCcca                                                                                                                                                <  2:372944/151‑1 (MQ=255)         cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCaccacc                                                                                                                                           <  1:1085298/151‑1 (MQ=255)          aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCg                                                                                                                                          >  1:1010221/1‑151 (MQ=255)           gACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGc                                                                                                                                         >  2:381734/1‑151 (MQ=255)            aCGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGcc                                                                                                                                        <  1:1648997/151‑1 (MQ=255)             cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGcc                                                                                                                                        >  2:1278752/1‑150 (MQ=255)              gcgtaaccagAAACATCATTCGCGTGTCTTCCGGAATTGCAGCCAGCGGTGTGGTGCCCGCTAACACACAGCCTGCGGTGGGCGTTACCAG‑ACCTG‑CCAGCAGGCGCCGCGGGGTTCTTTT‑G‑CCACCTCC‑GCTCCGGCCCACCGCCGCca                                                                                                                                       <  2:543763/141‑1 (MQ=255)                gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCcaca                                                                                                                                     >  1:236944/1‑150 (MQ=255)                     cTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTg                                                                                                                                >  1:1504217/1‑150 (MQ=255)                          aCATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                          >  1:451531/1‑151 (MQ=255)                             tcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCgg                                                                                                                       >  2:1467448/1‑151 (MQ=255)                                 tCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAat                                                                                                                    >  1:1023112/1‑150 (MQ=255)                                   ggttgtCTTCCTGATTTCCGCCCAGGGGGGTGGTCCCCGCAAACCCTCCGCCGGGGGTGGGCTTTCCCAG‑CCCGG‑CGAGCAGCCGCGGGGGGTTTTTTTT‑C‑CCCCCACC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAata                                                                                                                   <  2:497321/146‑1 (MQ=255)                                   gcgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                 >  1:141403/1‑151 (MQ=255)                                             cTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                       >  2:1483258/1‑151 (MQ=255)                                              tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                      >  1:268956/1‑151 (MQ=255)                                              tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                      >  1:1583745/1‑151 (MQ=255)                                                aaTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                    >  1:363144/1‑151 (MQ=255)                                                  tttcgcccaggGGTGGGGCCCCCCCTAACCCACCCCGTGCGGTTGGTTTTTCCAG‑CCCCG‑CCAGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                    <  1:1404455/143‑1 (MQ=255)                                                  ttccacccgggggTGGGGTCCCCGATAACCCCCCGCCGGGGTTGGGCTTTCCCGG‑CCCTC‑CCAGCGGGCGCAGGGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                      <  2:1211714/136‑1 (MQ=255)                                                        ccgggggTGGGGTCCCCCCAACCCCACCGCCGGGGGTTGGCTTTCCCAG‑CCCTG‑CCAGCAGCCGCAGCGGGGTTTTTTT‑C‑CCACCCCC‑GCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAg                                                                                               <  1:955520/143‑1 (MQ=255)                                                          agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                           >  1:818009/1‑150 (MQ=255)                                                           gggggggggtcCCCGCTAACCCACCGCCTGCGGTTGGCTTTTCCAG‑CCCCG‑CCAGCGGGCGCAGGGGGGTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                           <  1:898666/142‑1 (MQ=255)                                                             gttgtgGTCCCCCCAAACCCACCCCCTG‑GGTTGGGCTTTCCCG‑GCCCGCCCAGCGGGCGCGGCGGGTTTCTTTT‑G‑CCCCCACC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGata                                                                                        <  2:977494/148‑1 (MQ=255)                                                             ggTGGGGCGCCCTCAAACAC‑GCGCCCGGGGTTTGGCTTTCCAGCCCCCG‑CCAGGGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                       <  1:1682457/151‑1 (MQ=255)                                                              gtgtGGTCCCCGCTAACCCACCGCCGGGGTTTGGCTTTCCCAG‑CCCGC‑CCGGCGGGCGCGGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtt                                                                                      <  1:137687/151‑1 (MQ=255)                                                                gtggtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                    >  1:526859/1‑151 (MQ=255)                                                                  ggtcCCCCCAAACCCACCGCCTGGGGTTGGCTTTCCCAG‑CCCTG‑CCAGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                  <  1:1490134/151‑1 (MQ=255)                                                                   gtgCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcg                                                                                 >  2:1558699/1‑151 (MQ=255)                                                                      cccccaaccccttCGCCTGGGGTGGGCTTTCCCCG‑CCCGG‑CCGGCGGGCCAAGC‑GGGTTTTTTT‑C‑CCCCCCCCGGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                  <  2:1147888/135‑1 (MQ=255)                                                                       ccGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGt                                                                              >  1:1624704/1‑150 (MQ=255)                                                                            accccaccGCCTGGGGTTGGCTTTCCCAG‑CCCTG‑CCAGCAGGCGCAGCAGGTTCTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATg                                                                          <  1:701954/142‑1 (MQ=255)                                                                             acacATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                       >  1:1657392/1‑151 (MQ=255)                                                                              cccgtCGCCGGCTGTGGGCTTTCCCGG‑CCCGG‑CCGGCAGCCGCGGCGGGGTACTTTT‑G‑CCCCCTCC‑GGCGCGCCCCTCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                       <  2:1354270/146‑1 (MQ=255)                                                                                  tCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                   >  2:298328/1‑150 (MQ=255)                                                                                   cGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                 >  1:1354574/1‑151 (MQ=255)                                                                                     ccgggggTTGGCTTTCCCAG‑CCCCC‑CCGGCGGGCGCAGCGGGGTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTa                                                                  <  1:1655629/143‑1 (MQ=255)                                                                                     ccTGGGGTTGGCTTTCCCAA‑CCCGG‑CCAGCGGGCGCGGCGGGGTCCTTTG‑C‑CCCCCCCC‑GCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                               <  1:1686237/151‑1 (MQ=255)                                                                                     ccTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                               >  1:1460661/1‑151 (MQ=255)                                                                                        gggttgggcttccccAG‑CCCTC‑CCAGCGGGCGCAGCGGGTTTTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                   <  2:463629/131‑1 (MQ=255)                                                                                          ggTTGGCTTTTCCCG‑CCCGC‑CCAGCGGGCGCAGC‑GGGTTCTTTTCC‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                             <  1:381751/148‑1 (MQ=255)                                                                                           gtgggCTTTTCCAG‑CCCCG‑CCAGCGGGCGCAGC‑GGGTTTTTTT‑GCCCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                               <  1:1737952/142‑1 (MQ=255)                                                                                            ttGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                          >  1:1581758/1‑149 (MQ=255)                                                                                            ttGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                          >  1:1335189/1‑149 (MQ=255)                                                                                             tGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                       >  2:1041039/1‑151 (MQ=255)                                                                                              ggCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                       >  1:684316/1‑150 (MQ=255)                                                                                                cGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                    >  1:91682/1‑151 (MQ=255)                                                                                                  ttcccAG‑CCCTC‑CCAGCAGGCGCGGCGGGGTATTTTC‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                      <  1:1526524/144‑1 (MQ=255)                                                                                                  ttcccAG‑CCCGG‑CCGGCGGGCGGAGCGGGGTTTTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                     <  1:909416/145‑1 (MQ=255)                                                                                                   ttcccaccccTC‑CCAGCAGGCGCGGC‑GGGTTCTTTTCC‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaa                                                     <  1:690024/140‑1 (MQ=255)                                                                                                    ttcccagcccccccAGCGGGCGCGGCGGGGTTCTTTC‑G‑CCCCCCCC‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                         <  2:947453/132‑1 (MQ=255)                                                                                                    tccccaccccccccAGCAGGCGAAGCAGGGTTTTTTC‑C‑CCCCCCCC‑GCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                    <  2:48047/137‑1 (MQ=255)                                                                                                    cccAG‑CCCTG‑CCAGCGGGCGC‑GGGGGGTTCTTTTCC‑CCCCCCCC‑GCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                         <  2:904175/141‑1 (MQ=255)                                                                                                     tccccgccccgccgccagGCACAGGAGGGTTTTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaacaa                                                  <  2:777929/135‑1 (MQ=255)                                                                                                      ccccccccagcgggcagGCCCAGCGTGTTTCTTTC‑C‑CCACCGCC‑GCCGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                 <  2:1018088/137‑1 (MQ=255)                                                                                                      cccagcccgcccAGCGGGCGCAGCGGGGTTCTTTC‑C‑CCCCCCCC‑GCGCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaaca                                                   <  1:310103/138‑1 (MQ=255)                                                                                                       aG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                             >  2:649592/1‑151 (MQ=255)                                                                                                        g‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                            >  2:58631/1‑151 (MQ=255)                                                                                                          ccccgcccggcggccgcAGCAGGTTCTTTTC‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                              <  2:114815/135‑1 (MQ=255)                                                                                                          aCCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                           >  2:1160115/1‑151 (MQ=255)                                                                                                            ccgcccggcgggCGCGGCGGGGTATTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGc                                              <  1:1586535/137‑1 (MQ=255)                                                                                                            ccccccAGCGGGCGCAGCGGGTTTTTTTC‑C‑CCCCCCCC‑GCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGt                                            <  1:22537/145‑1 (MQ=255)                                                                                                            cTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                         >  2:924514/1‑151 (MQ=255)                                                                                                             tG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGCATATGTAAATCCAGTTGGATCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGccc                                        >  2:208300/1‑151 (MQ=255)                                                                                                                ccAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                      >  2:1166952/1‑151 (MQ=255)                                                                                                                   gcggccgcAGAGGGGTTCTTTG‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                    <  1:41627/145‑1 (MQ=255)                                                                                                                      ggCGCAGCAGGTTACTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTg                                      <  1:660494/145‑1 (MQ=255)                                                                                                                      ggCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                  >  2:750130/1‑149 (MQ=255)                                                                                                                       gcggAGCAGGTTTTTTTC‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                       <  2:68382/143‑1 (MQ=255)                                                                                                                        cgcAGCAGGGTATTTTT‑C‑CCCCCCCC‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                  <  1:1060672/147‑1 (MQ=255)                                                                                                                        cgcAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                              >  1:510596/1‑151 (MQ=255)                                                                                                                          cagcagGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                             >  2:712059/1‑150 (MQ=255)                                                                                                                            gcagGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                          >  2:625114/1‑151 (MQ=255)                                                                                                                              agGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                          >  1:1243525/1‑149 (MQ=255)                                                                                                                               ggtttcTTTT‑C‑CCCCCCCC‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                         <  1:1168309/144‑1 (MQ=255)                                                                                                                                 gTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                     >  1:771013/1‑151 (MQ=255)                                                                                                                                  tACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGGAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                    >  1:1711131/1‑151 (MQ=255)                                                                                                                                      ttt‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                 >  1:1618942/1‑150 (MQ=255)                                                                                                                                      ttt‑C‑CCCCCACC‑GGTGCGGCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                    <  1:1079013/147‑1 (MQ=255)                                                                                                                                       tttcccccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                   <  2:327208/138‑1 (MQ=255)                                                                                                                                       tttcccccccccc‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                   <  1:1016942/138‑1 (MQ=255)                                                                                                                                        ttcccccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                 <  1:475920/140‑1 (MQ=255)                                                                                                                                        t‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc              >  1:614399/1‑151 (MQ=255)                                                                                                                                            cccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt            <  2:1605681/145‑1 (MQ=255)                                                                                                                                            cccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt            <  1:901109/145‑1 (MQ=255)                                                                                                                                            ccaccaCC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt            >  1:422816/1‑151 (MQ=255)                                                                                                                                             cccccgccgcgtcGGCCCTCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                  <  2:360963/133‑1 (MQ=255)                                                                                                                                             ccccccc‑GCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAc               <  2:1749707/142‑1 (MQ=255)                                                                                                                                               ccccc‑GCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGc         <  1:814547/148‑1 (MQ=255)                                                                                                                                                cccc‑GCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc        <  2:302182/149‑1 (MQ=255)                                                                                                                                                cacc‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc        >  1:1693812/1‑151 (MQ=255)                                                                                                                                                cacc‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc        >  1:1591389/1‑151 (MQ=255)                                                                                                                                                 acc‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct       >  1:1719641/1‑151 (MQ=255)                                                                                                                                                        gCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTc  >  2:602701/1‑150 (MQ=255)                                                                                                                                                           |                                                                                                                                                   CACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAG‑ACCTG‑CCAGCAGGCGCAGCAGGGTACTTTT‑G‑CCACCACC‑GCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTC  >  NC_000913/993764‑994059

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