BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,714,163:1	+T	9.7%	coding (2232/2334 nt)	bisC ←	biotin sulfoxide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,714,163	1	.	T	9.7%	362.2 / 35.2	103	coding (2232/2334 nt)	bisC	biotin sulfoxide reductase
Reads supporting (aligned to +/- strand): ref base . (48/45); new base T (2/8); total (50/53)
Fisher's exact test for biased strand distribution p-value = 9.37e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAG  >  NC_000913/3714014‑3714309

caTTGTCGATTGCGCATGGCAGGCTGAAACCCAACTACCCTCATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGc                                                                                                                                                     <  2:1655060/150‑1 (MQ=255)

caTTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGc                                                                                                                                                     >  1:851072/1‑150 (MQ=255)

    gTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑cgcg                                                                                                                                                >  2:1372160/1‑150 (MQ=255)

    gTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑cgcg                                                                                                                                                >  2:1650427/1‑150 (MQ=255)

    gTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑cgcg                                                                                                                                                >  1:1719403/1‑150 (MQ=255)

      cGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACTGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAg                                                                                                                                              >  2:1500821/1‑150 (MQ=255)

      cGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAg                                                                                                                                              <  2:961672/150‑1 (MQ=255)

      cGATTGCGCATGGCAGGAGGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAg                                                                                                                                              <  2:1154486/150‑1 (MQ=255)

       gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGc                                                                                                                                             <  1:765811/150‑1 (MQ=255)

       gATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGc                                                                                                                                             <  1:1178257/150‑1 (MQ=255)

           gcgcATGGCAGGCTGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTggg                                                                                                                                         <  2:246340/150‑1 (MQ=255)

            cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGCTCGCG‑GCTgggg                                                                                                                                        <  1:1608656/150‑1 (MQ=255)

            cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGCTCGCG‑GCTgggg                                                                                                                                        <  2:1123338/150‑1 (MQ=255)

            cgcATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGCTCGCG‑GCTgggg                                                                                                                                        <  1:1523882/150‑1 (MQ=255)

              cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGgag                                                                                                                                      >  1:238900/1‑150 (MQ=255)

              cATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCGGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGgag                                                                                                                                      <  1:1110456/150‑1 (MQ=255)

                tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGCTCGCG‑GCTGGGgaatg                                                                                                                                    <  2:405044/150‑4 (MQ=255)

                tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGAt                                                                                                                                    >  2:1606633/1‑150 (MQ=255)

                tGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGGGTCAGTTCCTGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGCTCGCG‑GCTGGGgaatg                                                                                                                                    <  2:823633/150‑4 (MQ=255)

                 ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGCTCGCG‑GCTGGGgaatgg                                                                                                                                   <  2:670300/150‑5 (MQ=255)

                 ggCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATc                                                                                                                                   >  1:763152/1‑150 (MQ=255)

                   cAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCtt                                                                                                                                 <  2:1291943/150‑1 (MQ=255)

                     ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                               >  1:1443206/1‑150 (MQ=255)

                     ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                               >  1:1067353/1‑150 (MQ=255)

                     ggATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTg                                                                                                                               >  2:621850/1‑150 (MQ=255)

                        tGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCTGCTATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                                                                            <  2:1443249/150‑1 (MQ=255)

                        tGAAACCCAACTACCCACATGGAGTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTGCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTc                                                                                                                            <  2:1278123/150‑1 (MQ=255)

                         gAAACCCAACTACCCACATGGATTATGAGCTGGCCGGGGGGTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCa                                                                                                                           <  1:412988/150‑1 (MQ=255)

                          aaaCCCAACTTCCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAg                                                                                                                          >  2:364596/1‑150 (MQ=255)

                          aaaCCCAACTACCCACATGGATTATGAGCTGGCCGGGGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAg                                                                                                                          <  1:1538009/150‑1 (MQ=255)

                           aaCCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGc                                                                                                                         >  2:112604/1‑150 (MQ=255)

                             cccAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGAACTTTGGTCAGCAc                                                                                                                       >  2:1353008/1‑150 (MQ=255)

                                 aatACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTc                                                                                                                   <  2:822522/148‑1 (MQ=255)

                                     cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGCCCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCCCGTTCACTg                                                                                                               >  2:1699940/1‑150 (MQ=255)

                                     cccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTg                                                                                                               >  2:1657242/1‑150 (MQ=255)

                                      ccACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgc                                                                                                              >  1:619140/1‑150 (MQ=255)

                                       cacaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                             >  1:137492/1‑150 (MQ=255)

                                       cacaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcg                                                                                                             <  1:1573891/150‑1 (MQ=255)

                                        acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTCAGGGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                            <  2:42493/150‑1 (MQ=255)

                                        acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                            <  1:732327/150‑1 (MQ=255)

                                        acaTGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTAACCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTgcgc                                                                                                            >  2:944585/1‑150 (MQ=255)

                                            ggATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtt                                                                                                        >  1:737067/1‑150 (MQ=255)

                                              aTTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATTCCCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttt                                                                                                       >  2:60332/1‑150 (MQ=255)

                                              aTTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGtttt                                                                                                      >  2:145312/1‑150 (MQ=255)

                                               ttATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGttttt                                                                                                     <  1:1554012/150‑1 (MQ=255)

                                                 aTGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTAc                                                                                                   <  2:1491868/150‑1 (MQ=255)

                                                    aGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACaaa                                                                                                >  1:1046740/1‑150 (MQ=255)

                                                     gCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAAt                                                                                               >  1:1183359/1‑150 (MQ=255)

                                                       tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGTAATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CACCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATa                                                                                              >  1:1290373/1‑150 (MQ=255)

                                                       tGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATAc                                                                                             >  2:1222771/1‑150 (MQ=255)

                                                           cGGTGGTTCTAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCATCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGATCACTGCGCCGTTTTTACAAATACCGTc                                                                                         <  2:727080/150‑1 (MQ=255)

                                                                   cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCCGCGCCGTATAACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTGTACAAAAACCATCAGCgggtt                                                                                 >  2:306427/1‑147 (MQ=255)

                                                                   cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                 >  2:360930/1‑150 (MQ=255)

                                                                   cAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTa                                                                                 <  2:851072/150‑1 (MQ=255)

                                                                       cGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATc                                                                             <  1:400594/150‑1 (MQ=255)

                                                                         cTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCa                                                                           >  2:1358820/1‑150 (MQ=255)

                                                                           gTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                         <  2:475412/150‑1 (MQ=255)

                                                                           gTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGa                                                                         >  1:1017116/1‑150 (MQ=255)

                                                                            tAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGAt                                                                        >  2:342340/1‑150 (MQ=255)

                                                                              aGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCCGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATcc                                                                      <  2:928967/150‑1 (MQ=255)

                                                                                 gtCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGc                                                                   <  2:460130/150‑1 (MQ=255)

                                                                                  tCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑ACCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGcc                                                                  <  2:153857/150‑1 (MQ=255)

                                                                                     gTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGCTCCGGCCATg                                                               >  2:1610718/1‑150 (MQ=255)

                                                                                      ttCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAg                                                                          <  1:30362/138‑1 (MQ=255)

                                                                                       tCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGcc                                                             <  1:364596/150‑1 (MQ=255)

                                                                                          ggACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCt                                                          >  1:1637997/1‑150 (MQ=255)

                                                                                          ggACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCt                                                          >  1:40037/1‑150 (MQ=255)

                                                                                            aCCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTc                                                        <  2:1059953/150‑1 (MQ=255)

                                                                                             ccGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCg                                                       <  1:1372160/150‑1 (MQ=255)

                                                                                              cGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGt                                                      <  1:1606633/150‑1 (MQ=255)

                                                                                               gTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGATGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTg                                                     >  1:344314/1‑150 (MQ=255)

                                                                                                  gT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAAt                                                  <  1:398388/150‑1 (MQ=255)

                                                                                                     aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGc                                                >  1:447563/1‑150 (MQ=255)

                                                                                                     aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGc                                                <  1:773208/150‑1 (MQ=255)

                                                                                                     aTTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTACGTGAATGc                                                >  1:1661291/1‑150 (MQ=255)

                                                                                                        tttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCaaa                                             >  1:1604437/1‑150 (MQ=255)

                                                                                                         ttCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGTACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGCTCCGGCCATGCCCCTTCGTGAATGCAAAt                                            >  2:821740/1‑150 (MQ=255)

                                                                                                           ccAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCCCGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCGTCGTGAATGCAAATcc                                          >  1:1087443/1‑149 (MQ=255)

                                                                                                             aGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACg                                        <  1:621850/150‑1 (MQ=255)

                                                                                                              gCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTTGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCTGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGc                                       <  2:1719403/150‑1 (MQ=255)

                                                                                                                cATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCa                                     >  1:244632/1‑150 (MQ=255)

                                                                                                                cATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCa                                     >  2:38472/1‑150 (MQ=255)

                                                                                                                  tGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGG‑CGCGAGCTGGGGATATCTTAGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAgg                                   <  2:994122/150‑1 (MQ=255)

                                                                                                                     cAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGttt                                >  1:914517/1‑150 (MQ=255)

                                                                                                                         gcCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAAtt                            <  1:1657242/150‑1 (MQ=255)

                                                                                                                            gTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTccc                         >  1:594993/1‑150 (MQ=255)

                                                                                                                                tACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGc                     <  2:594993/150‑1 (MQ=255)

                                                                                                                                 aCCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCt                    <  2:619140/150‑1 (MQ=255)

                                                                                                                                 aCCCGCACAGCAATT‑CCCCTTC‑CTCGAGCTGGGGAGATCTTTGGTCAGCACGTTCTGTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCt                    <  2:1429160/150‑1 (MQ=255)

                                                                                                                                  cccGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTa                   >  2:1606309/1‑150 (MQ=255)

                                                                                                                                    cGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAAt                 <  2:763152/150‑1 (MQ=255)

                                                                                                                                     gCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATg                <  2:137492/150‑1 (MQ=255)

                                                                                                                                        caGCCATT‑CCCTAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAcc             <  2:1730764/150‑1 (MQ=255)

                                                                                                                                        caGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGAcc             <  1:1222771/150‑1 (MQ=255)

                                                                                                                                          gCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGc           >  2:908335/1‑150 (MQ=255)

                                                                                                                                          gCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGc           >  1:186041/1‑150 (MQ=255)

                                                                                                                                              ct‑CCCCAGCTCGCG‑GCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCg       >  1:1034616/2‑150 (MQ=255)

                                                                                                                                               t‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCgg      >  2:1234533/1‑150 (MQ=255)

                                                                                                                                                  cccAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCa    >  2:1719941/1‑150 (MQ=255)

                                                                                                                                                    cAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAg  <  1:944585/150‑1 (MQ=255)

CATTGTCGATTGCGCATGGCAGGATGAAACCCAACTACCCACATGGATTATGAGCTGGCCGGTGGTTCAAACGCTGTAAGTGTCAGTTCCGGACCGTTGT‑ATTTTTCCAGCCATGCCAGCGCCGTATTACCCGCACAGCCATT‑CCCCAGC‑CGCGAGCTGGGGAGATCTTTGGTCAGCACGTTCACTGCGCCGTTTTTACAAATACCGTCAGCGGTTAAATCCAGATCCGGCCATGCCCCTTCGTGAATGCAAATCACGCCAGGTTTAATTCCCTCGCTAATGACCGCTCCGGCAAG  >  NC_000913/3714014‑3714309

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 14 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`