BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	915,215	T→G	8.6%	intergenic (‑358/+137)	ybjE ← / ← aqpZ	putative transporter/aquaporin Z

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	915,215	0	T	G	8.6%	370.1 / 6.4	146	intergenic (‑358/+137)	ybjE/aqpZ	putative transporter/aquaporin Z
Reads supporting (aligned to +/- strand): ref base T (87/42); new base G (3/12); total (91/55)
Fisher's exact test for biased strand distribution p-value = 5.28e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 150 total aligned reads displayed.

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364

ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAAGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                          <  1:1342883/150‑1 (MQ=17)

ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAAGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAtt                                                                                                                                                          <  1:1342917/150‑1 (MQ=17)

  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                        <  2:766523/150‑1 (MQ=32)

  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                        <  1:3253455/150‑1 (MQ=32)

  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                        <  2:3116181/150‑1 (MQ=32)

  ggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTgg                                                                                                                                                        <  2:3513704/150‑1 (MQ=32)

   tcAGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑c                                                                                                                                                      >  2:3854088/2‑150 (MQ=17)

   gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑c                                                                                                                                                      <  2:1374942/150‑1 (MQ=32)

    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑cc                                                                                                                                                     >  2:1082782/1‑150 (MQ=32)

    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑cc                                                                                                                                                     <  2:1319924/150‑1 (MQ=32)

    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑cc                                                                                                                                                     >  2:3066511/1‑150 (MQ=32)

    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑cc                                                                                                                                                     >  2:3988026/1‑150 (MQ=32)

    cAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCACGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑cc                                                                                                                                                     <  1:1648113/150‑1 (MQ=17)

     aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccg                                                                                                                                                    <  1:3848483/150‑1 (MQ=32)

     aGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccg                                                                                                                                                    <  1:1735481/150‑1 (MQ=32)

     aGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGGCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTA‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccg                                                                                                                                                    >  1:444402/1‑150 (MQ=255)

     aGCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccg                                                                                                                                                    >  1:2379018/1‑150 (MQ=17)

      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgc                                                                                                                                                   <  1:3962927/150‑1 (MQ=32)

      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgc                                                                                                                                                   <  1:1917119/150‑1 (MQ=32)

      gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgc                                                                                                                                                   <  1:1916871/150‑1 (MQ=32)

       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgcc                                                                                                                                                  <  1:2046604/150‑1 (MQ=32)

       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgcc                                                                                                                                                  <  1:2046446/150‑1 (MQ=32)

       cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGg‑ccgcc                                                                                                                                                  >  2:1367662/1‑150 (MQ=32)

        ccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGccc                                                                                                                                                 >  1:2419014/1‑150 (MQ=17)

        ccACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGAGAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑CCGccc                                                                                                                                                 >  1:3733981/1‑150 (MQ=255)

         cactTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGcccc                                                                                                                                                <  2:2371550/150‑1 (MQ=17)

         cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGcccc                                                                                                                                                <  2:1941510/150‑1 (MQ=32)

          acaTGGATGTGGGCTGTGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGccccc                                                                                                                                               >  1:3585154/1‑150 (MQ=17)

           caTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑t                                                                                                                                             <  2:1184291/150‑1 (MQ=32)

             tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑Agc                                                                                                                                           >  1:3635438/1‑150 (MQ=17)

             tGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑Agc                                                                                                                                           >  1:3953866/1‑150 (MQ=17)

              ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑Tgcg                                                                                                                                          <  2:598601/150‑1 (MQ=32)

              ccaTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑Tgcg                                                                                                                                          >  2:497931/3‑150 (MQ=17)

                 tgtgggctgtgggcGGTTTTTACTGGTTTTTACCCTGGGGCCCGACCGGTTCGCCTTATGTTTTATGTGGTGGGCACCGCGCCGGGGCTTTTTTTTTC‑GCATTAGCCGGGGGGTGCTGGGGGGGGCGGCTTTTGG‑CCGCCCCC‑Tgcgcgc                                                                                                                                       <  2:1186407/150‑1 (MQ=2)

                      gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATAGG‑CCGCCCCC‑TGCGCGCTccat                                                                                                                                  >  1:1665112/1‑148 (MQ=17)

                       cTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCCCCCCC‑GGCGCGCTCCTTg                                                                                                                                 >  1:3696106/1‑150 (MQ=17)

                          aGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTgcgc                                                                                                                              >  2:1575177/1‑150 (MQ=32)

                            ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑CCGCCCCC‑CGCGCGCACCTTGCGccc                                                                                                                            >  1:1159529/1‑149 (MQ=12)

                            ggCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑GGCGCGCTCCGTGCGCCa                                                                                                                            >  1:1159624/1‑150 (MQ=17)

                                 gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑c                                                                                                                      >  2:3911269/1‑150 (MQ=33)

                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGG‑Ca                                                                                                                     >  2:142646/1‑150 (MQ=21)

                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGG‑Ca                                                                                                                     >  2:1521168/1‑150 (MQ=21)

                                  ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGAGTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑Ca                                                                                                                     >  1:3221713/1‑150 (MQ=21)

                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGG‑CAATATGt                                                                                                               >  2:1479815/1‑150 (MQ=255)

                                        aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAGTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCACTGG‑CACTATGt                                                                                                               >  2:1810754/1‑150 (MQ=255)

                                          gggtGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGG‑CCGCCCCC‑AGCGCGCGCCGTGCCCCAGAGG‑CAATATGTTg                                                                                                             >  1:3541053/4‑150 (MQ=14)

                                               ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATGGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAg                                                                                                        >  1:2830193/1‑150 (MQ=255)

                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTc                                                                                                     >  2:303338/1‑150 (MQ=255)

                                                  ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTc                                                                                                     >  2:3511666/1‑150 (MQ=255)

                                                      gcgcCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGa                                                                                                 >  2:675634/1‑150 (MQ=255)

                                                           cGACCGGGTAGCCGTACGGTATAAGTCGTGGGCACGGCGCAGCGGCTATTTTGTTC‑GCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGTTTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGa                                                                                            <  2:2447499/150‑1 (MQ=17)

                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:2834645/1‑150 (MQ=255)

                                                            gACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAg                                                                                           >  2:2452131/1‑150 (MQ=255)

                                                               ccGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGTGGGTGAAGGGGGCGGCCGCGAGTGG‑CCGCCCCC‑AGCGCGCTCCTTGCGCCAGTCG‑CAATATGTTGCTTAGCTCATGAAAGGAgcgc                                                                                        >  2:1535214/1‑150 (MQ=17)

                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                     >  2:1995388/1‑150 (MQ=255)

                                                                  gTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                                     >  2:1709473/1‑150 (MQ=255)

                                                                   tAGCCGTATGGTTTAAGTGGTGGGCACCGCGCAGGGGCTTTTTTGTTC‑GCCTGTGCCCGGGGGTGCTGGGGGGGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGcaaca                                                                                    <  2:679827/150‑1 (MQ=17)

                                                                   tAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGcaaca                                                                                    >  2:1429668/1‑150 (MQ=255)

                                                                    cgCCGTATGGGATACGTCGAGGGCCCCGCGCAGCGGCTTTTTTGTGC‑GCCTGCGCCGGGGGTGGCTGGGGGGGGGGGGGTTTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGcaacaa                                                                                   <  2:2733681/149‑1 (MQ=17)

                                                                      ccGTAAGGGATAAGTCGAGGGCACCGCGAAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                                 >  2:2383511/1‑150 (MQ=255)

                                                                        gTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑AGCGCGCTCCGTGCGCCAGAGG‑CCATATGGTGCGTAGCTAATGACAGGAGCGCAACAAgatg                                                                               >  1:3845554/1‑150 (MQ=14)

                                                                         tatgggATACGTCGTGGGCACCGCGCTGGGGCTTTTTTGTTC‑GCCTGAGCCGGGGGGTGCTGGGGGGGGGGGGGTTTGG‑CCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                              <  2:355390/147‑1 (MQ=255)

                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGGGGCGGCGGGTGG‑CCGCCCCC‑AGCCCGCTCCATGCGCCACTGG‑CAATATTTTTCTTATATGAAGTAATGAGCGCAACAAGATTATTTATAAACaa                                                                   >  1:2892380/1‑149 (MQ=9)

                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:1593803/1‑150 (MQ=255)

                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACat                                                                   >  2:2174919/1‑150 (MQ=255)

                                                                                    tCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCCATTGG‑CCCCCCCC‑AGCGCGCGCCATGCGCCAGTTG‑CAAAAATTTGCGTAGCTCATGAAACGGGCGCAACAAGATGATGAATaaaaaa                                                                   >  1:4132150/1‑147 (MQ=17)

                                                                                        gggCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCACCATGCGCCAGTGG‑CAATATGTTGCTTAGCTCAAGAAAGGAGCGCAACAAGATGATGAATCAACAAATaa                                                               >  2:1921414/1‑150 (MQ=255)

                                                                                        ggccACGGCGCAGCGGCTTTTTTGTTC‑GCCTGAGCCGGGGGGTGCTGGGGGCGGCGGCTTTTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaa                                                               <  2:2690470/147‑1 (MQ=21)

                                                                                          gCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaaca                                                             >  2:654837/1‑150 (MQ=255)

                                                                                            aCCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACATATaacaac                                                           >  2:1360191/1‑150 (MQ=255)

                                                                                              cgcgCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAAAATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAAACAACAGATAACACCAt                                                         >  2:2653899/1‑150 (MQ=255)

                                                                                              cgcgCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAATGG‑GCGCCCCC‑AGCGCCCGCCCTGCGCCAGTGGCCAAAATTTTG‑TTAGCGCAAGGAAGGAGCGCCACAACATTATGAATCAACAAAAAACAACAt                                                         >  1:4137521/1‑150 (MQ=9)

                                                                                                       gggtttttGTTTGGCCAGCGCGGGGGGGGGGGGGGGGGGGGGGCTTTGGG‑CCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaa                                                 <  2:146859/146‑1 (MQ=25)

                                                                                                       gcgATTTTGTTC‑GCCAGAGCCCGGGGGTGCCGGGGGCGGCGGCGCTTGG‑CCCCCCCC‑AGCGCGCGCCGTTCGCCACGGG‑CAATATTTTGCTTATCTCATGAAAGGAGCGCAACAAGAAGATGAAACAAAAAAAAACAACATAATaaaaaa                                                >  1:3085842/1‑150 (MQ=18)

                                                                                                        tttttttGTGC‑GCCCGCCCCGGGGGGTGCGGGGGGGGGGGGCGTGTGC‑CCCCCCCC‑GGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGGAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                               <  1:2797282/147‑1 (MQ=255)

                                                                                                         ttttttGTCC‑GCCTGCGCGGGGGGGTGCTGGGGGGGGGGGGTTTGGG‑CCCCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              <  2:3398175/148‑1 (MQ=25)

                                                                                                         gtttttGTCC‑GCCAGCGCGGGGGGGCGCGGGGGGGGGCGGGGTGGGG‑CCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                              <  2:3398886/148‑1 (MQ=255)

                                                                                                          ttttgGTTC‑GCCAGCGCGGGGGGGTGCTGGGGGGGGGGGCTTTTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAgg                                             <  2:3863140/149‑1 (MQ=25)

                                                                                                           ttttttTC‑GCCAGAGCCGGGGGGTGCTGGGGGGGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                            <  2:2414908/150‑1 (MQ=255)

                                                                                                           ttttGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCTTGG‑CCCCCCCC‑AGCGCCCGCCCTGCGCCACGGG‑CAAAATTTTTCATATCACATGCAAGGAGCGCAACAACAATATTAATCAACAGATAACAACATCATAAAAAAATGc                                            >  1:671596/1‑150 (MQ=12)

                                                                                                           ttttGTGC‑GCCTGCGCGGGGGGGTGGGGGGGGCGGGGGGGTTGGC‑CCCCCCCC‑CGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                            <  1:321489/150‑1 (MQ=25)

                                                                                                                ttC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAACGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                       >  2:3135303/1‑150 (MQ=255)

                                                                                                                ttC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGAc                                       >  2:3990901/1‑150 (MQ=255)

                                                                                                                     ccAGAGCGGGGGGGTGGGGGGGGCGGGGGCGTTGGC‑CCCCCCCC‑CCCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTa                                   <  1:622592/150‑1 (MQ=255)

                                                                                                                        gagCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                                >  2:554259/1‑150 (MQ=255)

                                                                                                                          ccccGGGGGCGCCGGGGGCCGCGCCGATCCC‑CCCCCCCC‑ACCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCa                              <  2:3277497/149‑1 (MQ=255)

                                                                                                                            ccGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                            >  2:3447978/1‑150 (MQ=255)

                                                                                                                              gggggTGCAGGGGGCGGAGGCGATTTG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGAGCAACAAGATGATGAATCAACATATAACAACAGCTTAAAAAAACGCCTGACATTACGCCAGGcc                          >  2:2537611/1‑150 (MQ=37)

                                                                                                                                  gTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGTGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCt                      <  2:2474724/150‑1 (MQ=255)

                                                                                                                                      ggggggCGGGGGGTTTGGG‑CGGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGt                  <  2:2149732/149‑1 (MQ=39)

                                                                                                                                       ggggggggCGG‑GTTTGGCCCCCCCCC‑CCCGCGCGCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt                 <  1:2903619/150‑1 (MQ=255)

                                                                                                                                        ggggCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAATAAGGCATGACATTACGCCAGGCCTTCTGCGtta                >  2:2248823/1‑150 (MQ=255)

                                                                                                                                          ggggggggcgTTGGG‑CCGCCCCC‑GGCGCGCGCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat              <  1:1638010/144‑1 (MQ=255)

                                                                                                                                          ggcggcggcGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaat              >  2:4201251/1‑150 (MQ=255)

                                                                                                                                               cggcgATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATc         >  2:3330949/1‑150 (MQ=255)

                                                                                                                                                ggcgATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  1:3064094/1‑150 (MQ=255)

                                                                                                                                                ggcgATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCa        >  2:2305561/1‑150 (MQ=255)

                                                                                                                                                   gATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATGACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     >  1:4252995/1‑150 (MQ=255)

                                                                                                                                                    tttGG‑CCCCCCCCGGGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGc     <  2:1110402/149‑1 (MQ=255)

                                                                                                                                                      tGG‑CCCCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCttt  <  2:3665064/150‑1 (MQ=255)

GGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAAGTCGAGGGCACCGCGCAGCGGCGATTTTGTTC‑GCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGATTGG‑CCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGG‑CAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATTAATCACGCTTT  >  NC_000913/915066‑915364

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 25 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`