BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	915,213	A→T	12.7%	intergenic (‑356/+139)	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,213	0	A	T	12.7%	266.5 / 11.3	123	intergenic (‑356/+139)	ybjE/aqpZ	putative transporter/aquaporin Z
Reads supporting (aligned to +/- strand): ref base A (75/26); new base T (2/15); total (81/42)
Fisher's exact test for biased strand distribution p-value = 1.39e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 123 total aligned reads displayed.

CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATT  >  NC_000913/915065‑915353

cGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt                                                                                                                                               <  2:3771773/150‑1 (MQ=32)

 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAC‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt‑t                                                                                                                                             <  2:1887485/150‑1 (MQ=17)

 ggggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAt‑t                                                                                                                                             <  2:112981/150‑1 (MQ=32)

  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                            <  2:3234682/150‑1 (MQ=32)

  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                            <  1:110212/150‑1 (MQ=32)

  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                            <  2:1910374/150‑1 (MQ=32)

  gggCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑Tg                                                                                                                                            <  2:1845883/150‑1 (MQ=32)

    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgc                                                                                                                                          <  1:217068/150‑1 (MQ=32)

    gCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTCACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgc                                                                                                                                          <  2:3577259/150‑1 (MQ=17)

       gCCCACCTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                       <  2:185200/150‑1 (MQ=17)

       gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                       <  1:3249302/150‑1 (MQ=32)

       gCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                       >  1:372529/1‑150 (MQ=32)

       gCCCACATCCATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                       >  2:3322037/1‑150 (MQ=17)

       gCCCACAGGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGGTACCGCGCGCCCGACCCGGTAGCCGTAAGGGATAC‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgc                                                                                                                                       <  1:1247452/150‑1 (MQ=255)

        cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgcc                                                                                                                                      >  1:1598147/1‑150 (MQ=32)

        cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgcc                                                                                                                                      >  1:2887995/1‑150 (MQ=32)

        cccACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGgccgcc                                                                                                                                      >  1:2888309/1‑150 (MQ=32)

          cacaTGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGcccc                                                                                                                                    <  2:3216345/150‑1 (MQ=32)

          cacaGGGTTGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGcccc                                                                                                                                    <  1:2092602/150‑1 (MQ=17)

               ggATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑Tgcg                                                                                                                              >  2:53672/1‑150 (MQ=32)

               ggATGTGGGCTGAGGGCGCGTTTTAAAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAA‑TGGCCGCCCCC‑Tgcg                                                                                                                              >  2:53595/1‑150 (MQ=17)

                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTccat                                                                                                                      >  2:3012395/1‑148 (MQ=17)

                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCtt                                                                                                                      >  1:2059347/1‑150 (MQ=32)

                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTccgt                                                                                                                      >  1:3756487/1‑148 (MQ=17)

                       gCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCtt                                                                                                                      >  2:1455653/1‑150 (MQ=17)

                         tGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTgc                                                                                                                    >  1:3156206/1‑150 (MQ=32)

                                  tttttACTGGGTTGTTACCTCGCGCCCGACCGGGTAGCGTTATGGTATAT‑GTCGTGGCCACCGCGCTGCGGCTTTTTTTTACGCAAGAGCCCGGGGGTGCGGGGGGCGGCGGCTTT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGc                                                                                                           <  2:579695/149‑1 (MQ=2)

                                  gTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGc                                                                                                           >  2:2579567/1‑150 (MQ=33)

                                   ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑TGCGCGCGCCATGCGCCCGTGGCa                                                                                                          >  1:338800/1‑150 (MQ=255)

                                   ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCa                                                                                                          >  1:2446567/1‑150 (MQ=34)

                                   ttttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCATGCGCCAGTGGCa                                                                                                          >  1:2739283/1‑150 (MQ=21)

                                    tttACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCaa                                                                                                         >  2:775605/1‑150 (MQ=35)

                                        ctggtttGTTCCCTGGCGCCCGACCCGGTCGCCGT‑TGGGTTAACGTGGTGGGCACCGCGCAGCGGCGTTTTTGTTCGCCCGCGCCGGGGGGTGCCGGGGGCGGCGGCGTT‑GGGCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATg                                                                                                     <  2:2526672/144‑1 (MQ=2)

                                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGt                                                                                                    >  1:253428/1‑150 (MQ=39)

                                         aGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTACATGCGCCAGTGGCAATATGt                                                                                                    >  1:1897433/1‑150 (MQ=255)

                                           ggATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGGGGCGTT‑TTGCCGCCCCC‑TGGGCGCTCCTTGCGCCAGTGGCAATATGTTg                                                                                                  >  2:1373214/1‑150 (MQ=255)

                                                ttACCTCGCGCCCGACCGGTTCGCGTTATGGTCAAA‑GTCGGGGGCACGGCGCCGCGGCTTTTTTTTTCGCCAGGGCCGGGGGGTGCCGGGGGGGTGGGCGTT‑TGGCCGCCCCC‑TCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAg                                                                                             <  2:1647926/150‑1 (MQ=9)

                                                ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTTTTCGCCAGAGCCCGGGGGTGCCGGGGGCGGCGGGGAT‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGGCAATAAGTTGCATAg                                                                                             >  1:3402687/1‑150 (MQ=14)

                                                ttACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTAGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TTGCCGCCCCC‑AGCCCGCACCCTGCGCCCGTGGCAATAATTTGCTTAg                                                                                             >  1:3402647/1‑150 (MQ=14)

                                                 tACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGc                                                                                            >  2:3735356/1‑150 (MQ=255)

                                                   ccTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTc                                                                                          >  2:728090/1‑150 (MQ=255)

                                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                         >  2:371296/1‑150 (MQ=255)

                                                    cTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCGTGCGCCAGTGGCAATATGTTGCTTAGCTCa                                                                                         >  2:371253/1‑150 (MQ=255)

                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCag                                                                                        >  1:3879567/1‑149 (MQ=255)

                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTAGCTCAt                                                                                        >  1:3879605/1‑150 (MQ=255)

                                                     tCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCGTGCGCCCGGGGCAATATGTTGCTTAGCTCAt                                                                                        >  1:1812638/1‑150 (MQ=255)

                                                     ccgcgCCCCACCCGGTTGCGGTACGGTTTTC‑GTGGGGGCCACCGCGCTGGGGGTTTTTTTTTCGCCAGCGCCGGGGGGGGCTGGGGGCGGGGGGTTT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCAt                                                                                        <  2:3805790/149‑1 (MQ=12)

                                                      cgcgCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                       >  2:3549649/1‑150 (MQ=255)

                                                      cgcgCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATg                                                                                       >  2:776079/1‑150 (MQ=255)

                                                          cccGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCCGTGGCAATATGTTGCTTTGCGCCTGAAAg                                                                                   >  1:3679929/1‑150 (MQ=14)

                                                            cGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGa                                                                                 >  1:2284225/1‑150 (MQ=255)

                                                             gACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                >  2:3050960/1‑150 (MQ=255)

                                                             gACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                >  2:3817314/1‑150 (MQ=255)

                                                             gACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCACCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAg                                                                                >  1:3775347/1‑150 (MQ=255)

                                                                  ggTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCACCGGAAAGGAGCGcaa                                                                           >  2:3697319/1‑150 (MQ=255)

                                                                   gTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGcaac                                                                          >  2:3581451/1‑150 (MQ=255)

                                                                       ccGTACGGTATAC‑GTGGAGGGCACCGCGACGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAga                                                                      <  2:3604670/150‑1 (MQ=18)

                                                                         gTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatg                                                                    >  2:1924740/1‑150 (MQ=255)

                                                                         gTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCCGCCAT‑TGGCCCCCCCC‑TGCGCGCGCCTTTCGCCACTGGCAATAATTTTCTTATCTCATTCAAGGAGCGCAACAAgaat                                                                    >  1:146865/1‑148 (MQ=17)

                                                                          tAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGGAGGGGGCGGCGGCGTT‑TGGGCGCCCCC‑TGCGCGCTCCTTGCGGCAGTGGCAATTTGTTGGTTTGCTGATGTTAGGTGCGCAACAAgattt                                                                   >  2:105784/1‑148 (MQ=14)

                                                                          tAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCGCCTTGCGCCACTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAgatga                                                                   >  1:3816009/1‑150 (MQ=255)

                                                                              ggATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGaa                                                               >  2:3529613/1‑150 (MQ=255)

                                                                              ggATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCGCCATGCGCCCCTGGCAAAAATTTGCTTAGCTAATGAAAGGAGCGCAACAAAATTATGaa                                                               >  1:2695989/1‑150 (MQ=17)

                                                                              ggATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCG‑CCCC‑TGCGCGCTCCTTGCCCCAGTGGCAAAAAGTTTCTTATCTCATGAAAGGAGCGCAAAAAGAATATGAAt                                                              >  2:2835342/1‑150 (MQ=17)

                                                                                   acgccgtgggCACCGCGCCGCGGCTTTTTTTTTCGCCAGAGCCGGGGGGTGCTGGGGGCGGCGGCGTT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCaa                                                           <  2:1997124/143‑1 (MQ=18)

                                                                                      tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGCT‑TGGCCGCCCCC‑AGCGCCCGCCCTGCGCCCCTGGCAAAATGTTTCATATCTCAGGGAAGGAGCGCAACAAGAAGAATAATAAACaa                                                        >  1:3529964/1‑149 (MQ=14)

                                                                                      tCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCCCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACat                                                        >  2:3015331/1‑150 (MQ=255)

                                                                                              aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATaacaac                                                >  2:813808/1‑150 (MQ=255)

                                                                                              aCCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCATGAGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCCACAAGATGATGAATCAACAAATaacaac                                                >  2:813836/1‑150 (MQ=255)

                                                                                                cgcgCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACAt                                              >  2:1735789/1‑150 (MQ=255)

                                                                                                     cgcggcgTTTTTGTTCGCCAGCGCCGGGGGGTGCTGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaa                                         <  2:3956706/149‑1 (MQ=255)

                                                                                                      ggggcttttttGTTCGCCTGAGCGGGGGGGGGGGGGGGGGGGCGGCTTT‑GGCCCGCCCCC‑CGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaa                                        <  1:3289126/143‑1 (MQ=255)

                                                                                                       gggcttttttttTCGCCAGAGCGGGGGGGTGCTGGGGGCGGCGGCTAT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaa                                       <  2:3868980/144‑1 (MQ=25)

                                                                                                         gcgATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGACTCAACATATAACAACATCTTaaaaaa                                     >  2:928315/1‑150 (MQ=255)

                                                                                                          tttttttGTTCGCCAGTGCGGGGGGGTGCGGGGGGGGGGGGCTTT‑GGCCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                    <  2:2759414/147‑1 (MQ=255)

                                                                                                          cttttttGTTCGCCAGCGCCGGGGGGTGCTGGGGGGGGGGGCGAT‑TGGCGGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTaaaaaaa                                    <  2:2130842/147‑1 (MQ=25)

                                                                                                           gtttttGTTCGCCAGAGCCCGGGGGTGCTGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   <  2:2972112/148‑1 (MQ=255)

                                                                                                           gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   >  2:122623/1‑150 (MQ=255)

                                                                                                           gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   >  2:464054/1‑150 (MQ=255)

                                                                                                           gATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAg                                   >  2:122581/1‑150 (MQ=255)

                                                                                                             ttttttTCGCCTGCGCCGGGGGGTGCGGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                 <  2:3920946/150‑1 (MQ=25)

                                                                                                             ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGTGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGc                                 <  2:883642/150‑1 (MQ=255)

                                                                                                             ttttGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACAGATAACAACATCTTAAAAAAAGGc                                 >  2:2743181/1‑150 (MQ=255)

                                                                                                                     gCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACAtt                         >  2:698874/1‑150 (MQ=255)

                                                                                                                      ccAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCATGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAAAATCTTAAAAAAAGGCCTGAAATTa                        >  2:1469379/1‑150 (MQ=255)

                                                                                                                       cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑GGGCCGCCCCC‑TGCGCGCTCCGTGCGCCAGGGGCAATATGTTGCGTAGCTCATGACAGGAGCGCAACAAGATGAAGCAGCAACATATCACAACATCTTCCAAAAAGGCCTGAccgtcc                       >  2:3519878/1‑145 (MQ=18)

                                                                                                                       cAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAG‑TGGCCGCCCCC‑TGCGCGCGCCTTGCGCCCGGGGCAAGATGGTGCGTCGCGCCGGAAAGGAGCGCACCAAGATGAGGACTCAACCTAGAACACCCGCTTACAAAAAGGCCTGACagtcc                       >  2:3519876/1‑146 (MQ=12)

                                                                                                                         gggggggggggTGGTGGGGGGGGTGGTTTT‑TGGCCGCCGAC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGATCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                     <  2:3402687/144‑1 (MQ=255)

                                                                                                                          ttgcggggggTGTGGGGGGGGGGGGGGTTGTGGCCCCCCCG‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGc                     <  2:3402647/145‑1 (MQ=255)

                                                                                                                          cgcgggggggTGGTGGGGGGGGGGGGTTT‑GGGCCCCCCCC‑TCCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGcc                    <  2:3720223/145‑1 (MQ=255)

                                                                                                                             ccGGGTGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAgg                 >  2:2610340/1‑150 (MQ=255)

                                                                                                                              cGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑AGCGCGCTCCTTGCTCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGAAGATGAATCAACAAATAACAACAACTTAAAAAAAGGCCTGACATTACGCCAGGc                >  2:1436878/1‑150 (MQ=37)

                                                                                                                                 gggTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCtt             >  1:2391853/1‑150 (MQ=255)

                                                                                                                                     gggggggggggggGCGGT‑GGGCCCCCCCCGGGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATGTAACAACATCTTAAAAAAAGGCCTGGCGTAATGTCAGGCCttttt          <  2:1159784/147‑4 (MQ=255)

                                                                                                                                      cAGGGGGCGGCGGCGAG‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCg        >  1:1497719/1‑150 (MQ=255)

                                                                                                                                        ggggggggCGGCGTT‑GGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGtt      <  2:2887995/150‑1 (MQ=255)

                                                                                                                                          gggCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaa    >  1:3331629/1‑150 (MQ=255)

                                                                                                                                            ggggcggccAA‑TCGCCGCCGCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt  >  2:251787/3‑150 (MQ=255)

                                                                                                                                            ggggcggccAA‑TCGCCGCCGCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt  >  1:3530897/3‑150 (MQ=255)

                                                                                                                                            gcgggggcgAT‑TGGCCGCCCCC‑GGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGttaatt  <  2:1288777/150‑1 (MQ=255)

CGGGGCAGCCCACATGGATGTGGGCTGAGGGCGCGTTTTACAGGGATGTTACCTCGCGCCCGACCCGGTAGCCGTAAGGGATAA‑GTCGAGGGCACCGCGCAGCGGCGATTTTGTTCGCCAGAGCCCGGGGGTGCAGGGGGCGGCGGCGAT‑TGGCCGCCCCC‑TGCGCGCTCCTTGCGCCAGTGGCAATATGTTGCTTAGCTCATGAAAGGAGCGCAACAAGATGATGAATCAACATATAACAACATCTTAAAAAAAGGCCTGACATTACGCCAGGCCTTCTGCGTTAATT  >  NC_000913/915065‑915353

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