BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	14.3%	intergenic (+112/+261)	yohP → / ← dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,229,177	0	A	C	14.3%	242.9 / 7.4	143	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base A (41/61); new base C (18/2); total (65/78)
Fisher's exact test for biased strand distribution p-value = 4.10e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229030‑2229326

cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑gg                                                                                                                                                       <  1:2931189/150‑1 (MQ=255)

cTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑gg                                                                                                                                                       <  1:1994805/150‑1 (MQ=255)

 tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCGGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑ggg                                                                                                                                                      >  2:4540032/1‑150 (MQ=255)

   gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGGCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑ggggg                                                                                                                                                    >  1:4494972/1‑150 (MQ=255)

    tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑GGGggc                                                                                                                                                   >  2:489961/1‑150 (MQ=255)

     ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCG‑GGGggc‑g                                                                                                                                                 >  1:4805741/1‑150 (MQ=255)

      gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGggc‑gg                                                                                                                                                <  1:2954342/150‑1 (MQ=255)

      gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGggc‑gg                                                                                                                                                <  1:1644608/150‑1 (MQ=255)

       tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGggc‑ggc                                                                                                                                               >  2:1235140/1‑150 (MQ=255)

       tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGggc‑ggc                                                                                                                                               <  1:4706602/150‑1 (MQ=255)

       tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGggc‑ggc                                                                                                                                               >  2:1283284/1‑150 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATACTCACTGGCGCAAGGAGCGC‑GCA‑CGGGGC‑GGCCaaa                                                                                                                                           >  2:4938449/1‑149 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑GGGGGG‑GGcccca                                                                                                                                           >  2:2350803/1‑147 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGcccag                                                                                                                                           >  1:658005/1‑147 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGGCaaa                                                                                                                                           >  1:1803532/1‑149 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑CGGGGC‑GGcccaa                                                                                                                                           >  2:1070288/1‑147 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCaaa                                                                                                                                           >  2:4567242/1‑149 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCaaa                                                                                                                                           >  2:337804/1‑149 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAAt                                                                                                                                           <  1:5177297/150‑1 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAAt                                                                                                                                           <  1:5177231/150‑1 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGCCAAAcg                                                                                                                                         >  1:2962321/1‑150 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGG‑GGcccaaag                                                                                                                                         >  1:2288397/1‑145 (MQ=255)

                tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCaaaatccc                                                                                                                                      >  1:439890/1‑144 (MQ=255)

                   ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑CGGGGC‑GGCCAATagccgccg                                                                                                                                   >  2:4979999/1‑150 (MQ=255)

                     aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCgccgccgcc                                                                                                                                 >  2:4967885/1‑150 (MQ=255)

                      aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAACATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGccc                                                                                                                                <  2:3497005/150‑1 (MQ=255)

                      aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGT‑GGGGGG‑GGCCAAGCGCCGcccccc                                                                                                                                >  2:3412012/1‑150 (MQ=255)

                      aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGccc                                                                                                                                <  1:3060410/150‑1 (MQ=255)

                      aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGccc                                                                                                                                <  1:3060423/150‑1 (MQ=255)

                       gatgatATTTTAAAATTAATTAATGTCATCAGGTCCTAATATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTATTATAATCACTGTCGCAAGGAGCGA‑GCA‑GGGGGC‑GGCCAATCGCCTCCGcccc                                                                                                                               <  2:735370/150‑1 (MQ=38)

                       gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCCAAACCCGCcccccc                                                                                                                               >  2:4263818/1‑150 (MQ=39)

                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCG‑GGGGGG‑GGCCAAAcccccccccccc                                                                                                                               >  1:4217834/1‑150 (MQ=255)

                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGCGGC‑GGCGATTGgccgccccctgc                                                                                                                              <  1:802070/150‑4 (MQ=39)

                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑CGGGGC‑GGCCAAACcccccccccccc                                                                                                                              >  2:4386761/1‑150 (MQ=39)

                         ttatatTTTAAAATTAATTAATGTCATCTGGTCCTAAAATAACGATAATATTTCAGTCTCTCATCCTTTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCt                                                                                                                             <  1:1869633/148‑1 (MQ=255)

                          tatatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTg                                                                                                                            <  1:1666324/149‑1 (MQ=255)

                               tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCACACGCCGCCCCCCCCCCCAccc                                                                                                                       >  1:2683855/1‑150 (MQ=38)

                                   aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑GGGGGCGGGCCCAAC‑CCCCCCCCCCCAGCAccccccg                                                                                                                   >  1:4412059/1‑148 (MQ=255)

                                   aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGCCCAACGCCGCCGCCCCCAGCAccccccc                                                                                                                   >  1:5466251/1‑148 (MQ=255)

                                    aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGA‑GGGGGC‑GGCCCAACGCCCCCGCCCCCGGCCCCCCCcgg                                                                                                                  >  2:5367771/1‑150 (MQ=255)

                                    aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:2810680/150‑1 (MQ=255)

                                     attaattaATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATGCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGG‑GGGGGG‑GGGCGATCGGGGGCGCCCCCAGGGCCCCCgggg                                                                                                                 >  2:736643/1‑149 (MQ=255)

                                         atttatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                             <  2:4353288/150‑1 (MQ=255)

                                         attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                             <  2:4352985/150‑1 (MQ=255)

                                          ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGGG‑GGGGGG‑GGCCAACCCCCCCCCCCCCCCCCCCCCCCggcgcagg                                                                                                             >  1:5238132/1‑144 (MQ=255)

                                              tGTCATCAGGTCCGAAAATAACTAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGaa                                                                                                       <  2:965143/150‑1 (MQ=255)

                                              tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGCG‑GGGGGG‑GGCCAAAC‑CCCCCCCCCCCAGCAACCCCCGGCGCTGG‑Cgca                                                                                                       >  2:2613633/1‑148 (MQ=255)

                                              tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCCATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGaa                                                                                                       >  2:920197/1‑150 (MQ=255)

                                              tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGAGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑CGGGGC‑GGCCAATCGCCGCCCCCCCCAGGCCCCCCCGGCACTGG‑CCaa                                                                                                       >  2:3534362/1‑150 (MQ=255)

                                                           gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGA‑GGGGGG‑GGGCAATGAGCGCCGCCCCCTGGACCCCCGGGCTCTGG‑GGAACAAAATCGCCGCt                                                                                          >  2:4504000/1‑150 (MQ=25)

                                                              aaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTgcg                                                                                       <  2:5472260/150‑1 (MQ=255)

                                                                tAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCCATCGCCCCCGCCCCCACCACCCCCCGGCGCAGG‑CGAACCAAAACGCCCCGGcgcg                                                                                     >  2:3630756/1‑150 (MQ=21)

                                                                     agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGCCAAAAGCCGCCGCCCCCACCACCCCCCGGCTCTTG‑CGAACAAAAAAGCCCCAGCGCGGTggc                                                                                >  2:2988149/1‑148 (MQ=255)

                                                                     agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCTCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGcc                                                                                >  1:4851348/1‑150 (MQ=255)

                                                                     agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑‑GGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGCGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGccc                                                                               >  1:3498448/1‑150 (MQ=255)

                                                                        atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:4495590/150‑1 (MQ=255)

                                                                          ttttCAGCCTCTCATCCTGTTGGGCTCCTGTCTGGTGCTTTGCTTCATATAATCCGTGGCGCAGGGAGCGC‑GCT‑GGGGGC‑GGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  1:1293866/149‑1 (MQ=18)

                                                                          atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:5191043/150‑1 (MQ=255)

                                                                            ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑CGGGGC‑GGCCAAACGCCGCCCCCCCCAGCACCCCCCGGCTCAGG‑CGCACAAAAACGCCCCAGCGCGGTGcccacccca                                                                         >  2:2780787/1‑144 (MQ=18)

                                                                             tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAAACCCCCCCCCCCCCCGCCCCCCCCGGCCCAGC‑AGCACAAAAACACCCCAGCGCCGGGcccaccccag                                                                        >  2:2153227/1‑143 (MQ=16)

                                                                                gTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:151225/150‑1 (MQ=255)

                                                                                 tctctcATCCTGTTGCGCTCGTGTCTTGTGAATTGCTTCATATAATCACTGGCGCACGGTGTGC‑GCT‑GGGGGC‑GCCAAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  1:2063667/150‑1 (MQ=18)

                                                                                 tctctaATCCTGTTGCGCTCCTGTCTTGTGCTTTGATTAATATAATCACTGGCGCAAGGAGCGC‑GCT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATcc                                                                    <  1:4278234/150‑1 (MQ=17)

                                                                                   tctcATCCTGTGGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                  <  2:4394989/150‑1 (MQ=255)

                                                                                     tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:2463004/150‑1 (MQ=255)

                                                                                     tcATCCTGTTGCGCTCCTGTCATGTGATTTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCT‑GGGGGC‑GGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:4567242/150‑1 (MQ=17)

                                                                                         ccTGTGGCGCTCGTTTCATGTGAATTGCTTAATATAATCACTGGCGCATGGAGCGC‑GCT‑GGGGGC‑GGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGc                                                            <  1:489961/150‑1 (MQ=17)

                                                                                             ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:532514/150‑1 (MQ=255)

                                                                                               gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:4365743/150‑1 (MQ=255)

                                                                                               gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCTGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:4618399/150‑1 (MQ=255)

                                                                                                cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCAGGGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:4267471/150‑1 (MQ=255)

                                                                                                  cTCCTGTCATGTGCTTTGCTTCATATAATCCCTGGCGCAAGGTGCGC‑GCT‑GGGGGC‑GGCCAATCGCGGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAAAAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  1:4572635/150‑1 (MQ=17)

                                                                                                  cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGA‑GGGGGG‑GGCCAAAAGCCCCCCCCCCCAGCACCCCCCGGCTAAGG‑CCAACAAAAAAGCCCCTGCCCCGTGCCCACCGCTAAACCCATACCGCGAccccggc                                                   >  2:3948665/1‑145 (MQ=255)

                                                                                                  cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGCCCAAAGCCGCCCCCCCCAGCACCCCCCGGCTCTGG‑CGCACAAAAACCCCCCGGCGCGGTGCCCACGGCATAACCCATACGGCGACCCGGTc                                                   >  2:3948613/1‑150 (MQ=255)

                                                                                                   tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCg                                                  <  1:3441248/150‑1 (MQ=255)

                                                                                                    ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  2:4805741/150‑1 (MQ=255)

                                                                                                     cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GGC‑CGGGGC‑GGCCCAACCCCCCCCCCCCCCCCAACCCCCCCCTCTGGCCGAAAAAAAAC‑CCCCGGCGCGGGGCCCGCCGCGTAAAACCTAAGGCTAACCGGTCggg                                                >  2:1678642/1‑150 (MQ=255)

                                                                                                       gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                              <  1:4862995/150‑1 (MQ=255)

                                                                                                          aTGTGAATTGATTCATATAATCACTGGCGCATGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:2245859/150‑1 (MQ=255)

                                                                                                           tgtgCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:1550296/150‑1 (MQ=255)

                                                                                                            gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  2:3272551/150‑1 (MQ=255)

                                                                                                            gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAAAAGCCCCCCCCCCCAGCACCCCCCGGCTCAGG‑CGAACAAAAAAGCCCCGGCGCCGTGCCCACGGCGTAAACCCAACCGCAAACCGGTCGGGCGCGAgg                                         >  2:1329064/1‑150 (MQ=9)

                                                                                                                    cTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:398112/150‑1 (MQ=40)

                                                                                                                     ttCATATAATCACTGGCGCAAGGAGCGC‑GCC‑GGGGGC‑GGCCAATAGCCGCCGCCCCCAGCACCCCCCGGCGCATG‑CCAAAAAAAAAGCCCCGGCGCGGTGGCCACGGCGTATACCATAAGGCTACCCGGTCGGGCGCGAGGTaaaaaacc                                >  2:1767527/1‑144 (MQ=255)

                                                                                                                      tCATATAATCACTGGCGCACGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2552209/150‑1 (MQ=255)

                                                                                                                             aTCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:4680500/150‑1 (MQ=33)

                                                                                                                             aTCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:3920074/150‑1 (MQ=33)

                                                                                                                             aTCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCACCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:4504000/150‑1 (MQ=18)

                                                                                                                                aCTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  2:2691983/150‑1 (MQ=32)

                                                                                                                                        aaGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  2:4217834/150‑1 (MQ=32)

                                                                                                                                         aGGAGCGC‑GCT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:1548889/150‑1 (MQ=17)

                                                                                                                                          ggAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  1:1708717/150‑1 (MQ=32)

                                                                                                                                          ggAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  1:1708819/150‑1 (MQ=32)

                                                                                                                                            aGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCTCGAGGTAACATACCTGTAAAACGCGCCCTCAGCCcac         <  2:3549212/150‑1 (MQ=17)

                                                                                                                                              cgc‑gcA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCTCAt       <  1:528360/150‑1 (MQ=17)

                                                                                                                                                c‑gcA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     >  2:2931647/1‑150 (MQ=32)

                                                                                                                                                c‑gcA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     >  2:2931554/1‑150 (MQ=32)

                                                                                                                                                  gcA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa    <  1:766552/150‑1 (MQ=32)

                                                                                                                                                    a‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:3410110/1‑150 (MQ=32)

                                                                                                                                                    a‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:1903356/1‑150 (MQ=32)

CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGG‑CGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229030‑2229326

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