BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	10.1%	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	10.1%	577.5 / 7.0	290	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base A (59/158); new base C (27/8); total (99/191)
Fisher's exact test for biased strand distribution p-value = 2.37e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 294 total aligned reads displayed.

CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229030‑2229326

    tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTGTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑Cggggggg                                                                                                                                                  >  1:7896031/1‑149 (MQ=255)

      gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgg                                                                                                                                                <  1:7137290/150‑1 (MQ=255)

       tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcggc                                                                                                                                               <  1:11884929/150‑1 (MQ=255)

         ccTGGCGTCTGTCAGATGCTCGTTGAACATTCATGACTGTCATCAGGTCCGACAAGCACGAGAATATTTCAGTCTCTCATCCTGTGGCGCTCCTGTCATGTGCATTGCTTCATATCCTCACTGGCGCAAGGAGCGC‑G‑CAGGGggcggcgg                                                                                                                                             <  2:8264578/149‑3 (MQ=25)

         aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCa                                                                                                                                             <  1:7402123/150‑1 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGGGGGCCaa                                                                                                                                            >  1:86487/1‑150 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CGgggggggggccca                                                                                                                                           >  1:10613367/1‑145 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAt                                                                                                                                           >  1:10104905/1‑150 (MQ=255)

           tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAt                                                                                                                                           <  2:10721608/150‑1 (MQ=255)

            ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgggccaaa                                                                                                                                          >  1:1841560/1‑144 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑gggggggggggccccac                                                                                                                                         >  2:11735408/1‑143 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCGGGggcgggcccacc                                                                                                                                         >  1:2737223/1‑143 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCGGGGGGGGcccaaag                                                                                                                                         >  2:1436788/1‑145 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGGGGCCaaaaa                                                                                                                                         >  2:4817633/1‑147 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGcccaacg                                                                                                                                         >  2:1862651/1‑145 (MQ=255)

             gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcg                                                                                                                                         >  2:653982/1‑150 (MQ=255)

               gTATTTAAGATGATATTTTAAACTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgcc                                                                                                                                       <  1:3348805/150‑1 (MQ=255)

               gTATTTAAGATGATATTTTAAAATTAAATAATGTCATCATGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCAGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgcc                                                                                                                                       <  2:7316516/150‑1 (MQ=255)

                 aTTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgccgc                                                                                                                                     <  1:4940763/150‑1 (MQ=255)

                      aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGCGGCGGCGATTGgccgccccct                                                                                                                                <  2:1064852/150‑2 (MQ=39)

                       gatttttttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCGTGTATTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  1:10455460/143‑1 (MQ=38)

                       gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGGGGGCCCA‑CCCCCCCGcccc                                                                                                                               >  1:6265017/1‑150 (MQ=255)

                       gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGGCAATAGcccccccccc                                                                                                                               >  2:3815805/1‑150 (MQ=39)

                       gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGcccc                                                                                                                               <  2:3313230/150‑1 (MQ=255)

                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CATGGGGCGGCCAATCGCCGCCGccccc                                                                                                                              >  2:11972881/1‑150 (MQ=255)

                        atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAACcccccccccccc                                                                                                                              >  2:2283418/1‑150 (MQ=39)

                         tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGATAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                             <  1:3657903/150‑1 (MQ=255)

                          gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCGGGGGGGGGCCAGAGCCGcccccccccg                                                                                                                            >  1:9574390/1‑148 (MQ=255)

                          gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGACATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GAGGGGGCGGGCCATAccccccccccccac                                                                                                                            >  1:7941562/1‑148 (MQ=255)

                              ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATAGCCGCCGCCCCCGGCAcc                                                                                                                        >  2:3081719/1‑150 (MQ=255)

                                ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CGGGGGGGGGCCAA‑AGCCCCCCCCCCCACCAcccc                                                                                                                      >  1:9815084/1‑150 (MQ=255)

                                ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GGGGGGGGGGGCCAACCCCCCCCCCCCCCCCAcccc                                                                                                                      >  2:526980/1‑150 (MQ=255)

                                ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GAGGGGGCGGCCCAAACCCCCCCCCCCCAGcccccc                                                                                                                      >  1:9100537/1‑150 (MQ=255)

                                ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAcccc                                                                                                                      <  2:4010049/150‑1 (MQ=255)

                                      ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑GGGGGGGCGGCCAA‑CACCCCCCCCCCCCGCAcccccccccg                                                                                                                >  1:8362459/1‑145 (MQ=255)

                                      ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGGCCAGCCCCGCCCCCCCCAGCCCCCCCCGGcg                                                                                                                >  2:1631064/1‑149 (MQ=38)

                                        aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGctc‑t                                                                                                             >  2:12202240/1‑150 (MQ=255)

                                        aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCCGCCAAACGCCCCCCCCCCCAGCCCCCCCCGGctca                                                                                                              >  2:9035641/1‑149 (MQ=37)

                                         attaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCGGGGGCGGGCAAACGCCGCCCCCCCCCGCCCCCCCcggggcgg                                                                                                             >  2:2727626/1‑145 (MQ=255)

                                             atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGGGGGCCAA‑CGCCGCCCCCCCCCCCCCCCCCCGGCCC‑AGGCgc                                                                                                        >  2:11619393/1‑149 (MQ=255)

                                              tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCCCCGCCCCCAGCACCCCCCGGCGC‑AGGCGaa                                                                                                       >  2:10025132/1‑150 (MQ=37)

                                               gTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GGGGGGGGGGGCAAAACCCCCCCCCCCCACCACCCCCCCGCGC‑AGGCGacc                                                                                                      >  1:6434031/1‑148 (MQ=255)

                                                   tcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACaaaa                                                                                                  <  2:7987946/150‑1 (MQ=255)

                                                     aGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGGCCAACCCCCCCCCCCCCACCCCCCCCCCGCAC‑AGGCGCCCaaaaaa                                                                                                >  1:8918864/1‑148 (MQ=255)

                                                        tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCAGATAATCACTGGCGCAAGGAGCGCGG‑GCGGGGGGGGCCCA‑CGCCCCCCCCCCCCCCCCCCCCCCGCGC‑GGGCGCACaaaaaaccc                                                                                             >  1:5748811/1‑145 (MQ=255)

                                                            aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAAAGCCGCCCCCCCCCCCACCCCCCGGCTC‑AGGCGAACAAAAAagccgcgg                                                                                         >  2:4991500/1‑148 (MQ=25)

                                                                  aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGt                                                                                   <  2:1011429/150‑1 (MQ=255)

                                                                      gaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCCCCCCCCCCAGCACCCCCCGGCTC‑AGGCGAACAAAAACGCCGCTGCGCGGTGccc                                                                               >  2:10249271/1‑150 (MQ=21)

                                                                        atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCAGTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGGGCCCTc                                                                             <  2:9241056/150‑1 (MQ=255)

                                                                          atTTCAGTCTCTCATCCTGTTGCGATCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCATGGTGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:404106/150‑1 (MQ=255)

                                                                            ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CGGGGGGGGGCCAAAAGCCCCCGCCCCCCGCACCCCCCGGGGC‑AGGCGAACAAAAAACCCCCGGCGCCGGGCCCCCGGcg                                                                         >  2:6481664/1‑149 (MQ=18)

                                                                            ttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCt                                                                         <  2:6345664/150‑1 (MQ=255)

                                                                             tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                        <  2:864781/150‑1 (MQ=255)

                                                                             tCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAAACACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCGGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCtt                                                                        <  2:4571733/150‑1 (MQ=255)

                                                                              cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  2:9622352/150‑1 (MQ=255)

                                                                                   tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAAAGCCGCCCCCCCCAGCACCCCCCGGCGA‑TGGCGAAAAAAAACGCCCCTGCGCCGTGCCCCCGGCATATacca                                                                  >  2:6915669/1‑149 (MQ=18)

                                                                                   tctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CACGGGGCGGCCAAACGCCCCCGCCCCCAGCACCCCCCGGCGC‑TGGCGCACAAAATCGCCGCGGCGCGGTGCCCTCGGCGTATaccg                                                                  >  2:9211100/1‑149 (MQ=18)

                                                                                             ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:6732003/150‑1 (MQ=255)

                                                                                             ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  2:9758077/150‑1 (MQ=255)

                                                                                               gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  2:8897092/150‑1 (MQ=255)

                                                                                                cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCggg                                                     <  2:6273480/150‑1 (MQ=255)

                                                                                                  cTCCTGTCATGTGCATGGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:5162134/150‑1 (MQ=255)

                                                                                                        tCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcgc                                             <  2:7380609/150‑1 (MQ=255)

                                                                                                          aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:8447616/150‑1 (MQ=255)

                                                                                                          aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:10025132/150‑1 (MQ=255)

                                                                                                           tgtgCATTGCTTCATATAATCACTGGCGCACGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:201073/150‑1 (MQ=255)

                                                                                                           tgtgCATTGCTTCATATAATCACTGGCGCAAGGTGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:11250943/150‑1 (MQ=255)

                                                                                                            gtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAgg                                         <  1:860631/150‑1 (MQ=255)

                                                                                                                  tgattCATATAATCACTGGCTCATGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACAt                                   <  1:8037858/147‑1 (MQ=255)

                                                                                                                  tGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGCCAAACTCCCCCCCCCCCAGCACCCCCCCGCGC‑TGGCCCAAAAAAAAGCCCCAGCGCGGTGCCCACGGCGTATACCATAACGCGACCCCGTCGGGCGCGCGGtcaaaa                                   >  2:6118728/1‑145 (MQ=255)

                                                                                                                    cTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATAGCCGCCGCCCCCAGCACCCCCGGGCTC‑TGGCGGACAAAATCGCCGCTGCGCGGTGGCCGCGGCTTATACCATAAGGCTACCCGGTGGGGCGCGAGGTTACAtac                                 >  2:10086456/1‑148 (MQ=14)

                                                                                                                      tCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:5820387/150‑1 (MQ=39)

                                                                                                                      tCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:12103869/150‑1 (MQ=39)

                                                                                                                      tCATATAAGCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:6619326/150‑1 (MQ=255)

                                                                                                                         tataATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTa                            <  1:4419101/150‑1 (MQ=37)

                                                                                                                          ataATCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTac                           <  1:3081719/150‑2 (MQ=255)

                                                                                                                            aaTCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGCCAAAAGCCCCCCCCCCCCGCAACCCCCGGCGCGGGGCGAAAAAAAAC‑CCCCGGCGCCGTGGCCACGGCGTTTACCATAAGGCAACCCGGGCGGGCGGGAGGGCAAAACCCAGTaaaa                         >  2:11009026/1‑150 (MQ=2)

                                                                                                                            aaTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  2:7603015/150‑1 (MQ=34)

                                                                                                                            aaTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATAGCCGGCGCCCCCACCACCCCCCCGCGC‑TGGCCAACAAAAACCCCGCGGCGCCGGGGCCACCGCCTATACCATACCGCGACCCGGGCGGGCGCGGGGTAACCACCCAGTaaaa                         >  2:4876020/1‑150 (MQ=2)

                                                                                                                             aTCACTGGCTCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:1225754/150‑1 (MQ=18)

                                                                                                                             aTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  2:9923747/150‑1 (MQ=33)

                                                                                                                              tCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcg                       <  1:5976505/150‑1 (MQ=32)

                                                                                                                               cACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:258574/150‑1 (MQ=32)

                                                                                                                                aCTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:808600/150‑1 (MQ=32)

                                                                                                                                aCTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:1571145/150‑1 (MQ=32)

                                                                                                                                    gcgcAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCt                 <  1:2515849/150‑1 (MQ=32)

                                                                                                                                       cAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:11299918/150‑1 (MQ=17)

                                                                                                                                       cAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:8179446/150‑1 (MQ=32)

                                                                                                                                        aaGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:2254060/150‑1 (MQ=17)

                                                                                                                                        aaGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  2:10834017/150‑1 (MQ=32)

                                                                                                                                          ggAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  1:5549917/150‑1 (MQ=17)

                                                                                                                                           gAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  2:5837058/150‑1 (MQ=32)

                                                                                                                                            aGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac         <  1:3952541/150‑1 (MQ=32)

                                                                                                                                             gcgc‑g‑cATGGGGCGGCCAATAGCCCCCCCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATAGCCGCTGCGCGGTGCCCACGGCGTATCCCATACCGCTACCCGGTCGGGCGCGCGGTAACAACCCGGTAAAACCCCCCCACAGCCccca        >  2:7835271/1‑150 (MQ=11)

                                                                                                                                              cgc‑g‑cTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAt       <  1:1436788/150‑1 (MQ=17)

                                                                                                                                                c‑g‑cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACAtgg     <  1:6935361/150‑3 (MQ=17)

                                                                                                                                                  g‑cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCa    <  2:1379948/150‑1 (MQ=32)

                                                                                                                                                    cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:3238496/150‑1 (MQ=32)

                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:807272/1‑150 (MQ=32)

                                                                                                                                                     aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:9411138/1‑150 (MQ=32)

CTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTC‑TGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  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: `‑`