BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	15.7%	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	15.7%	561.2 / 21.7	322	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base A (68/158); new base C (40/12); total (123/199)
Fisher's exact test for biased strand distribution p-value = 8.83e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 341 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

    ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑Agggg                                                                                                                                                      <  2:8879955/150‑1 (MQ=255)

     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑Aggggg                                                                                                                                                     >  2:7801281/1‑150 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑GGGGggc                                                                                                                                                    >  1:10001273/1‑150 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGggc                                                                                                                                                    >  1:3181176/1‑150 (MQ=255)

       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑Cggggggg                                                                                                                                                   >  2:4297102/1‑148 (MQ=255)

       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcg                                                                                                                                                   <  2:5936690/150‑1 (MQ=255)

       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcg                                                                                                                                                   >  2:3494656/1‑150 (MQ=255)

       ggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcg                                                                                                                                                   >  2:8967095/1‑150 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcgg                                                                                                                                                  <  1:170365/150‑1 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcgg                                                                                                                                                  >  2:2795022/1‑150 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcgg                                                                                                                                                  <  1:914314/150‑1 (MQ=255)

         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcggc                                                                                                                                                 <  1:571180/150‑1 (MQ=255)

          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑Gg‑ggggggggggg                                                                                                                                                >  2:3742170/1‑148 (MQ=255)

          gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGc‑‑c                                                                                                                                              >  2:2252619/1‑150 (MQ=255)

           aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑a                                                                                                                                                          >  1:8182739/1‑139 (MQ=255)

            cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑Caa                                                                                                                                            <  2:8143289/150‑1 (MQ=255)

             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGCGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑ACGGGGCGGC‑‑Caaa                                                                                                                                           >  2:605210/1‑149 (MQ=255)

               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑GGGGGGCGGc‑‑ccaaac                                                                                                                                         >  1:6936623/1‑145 (MQ=255)

               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGgccgc‑‑cccacg                                                                                                                                         >  2:9785148/1‑145 (MQ=255)

                     ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAAAcgccgccg                                                                                                                                   >  2:7024574/1‑150 (MQ=255)

                     ttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGCGGCGGC‑‑GATTggccgccc                                                                                                                                   <  1:10860838/150‑2 (MQ=255)

                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATcgccgccg‑c                                                                                                                                 <  2:6229274/150‑1 (MQ=255)

                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCgccgccg‑cc                                                                                                                                <  1:8253093/150‑1 (MQ=255)

                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCgccgccg‑cc                                                                                                                                <  2:4950258/150‑1 (MQ=255)

                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCgccgccg‑cc                                                                                                                                <  1:1070380/150‑1 (MQ=255)

                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑ccc                                                                                                                               <  1:2737160/150‑1 (MQ=255)

                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑ccc                                                                                                                               <  2:9442545/150‑1 (MQ=255)

                        aGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGCGGCGGC‑‑GATTGgccgccc‑cct                                                                                                                               <  1:5979632/150‑2 (MQ=39)

                         tattatatTTTAAATTTAATTAATGTCATCTGGTCCGAAAATAACTAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGAATTGATTAATATAATCACTGGCGCAAGGTGCGC‑GC‑TGGGGGCGGC‑‑CAATCGCCGCCG‑cccc                                                                                                                              <  1:3809061/146‑1 (MQ=38)

                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCCCCGCcc‑cccc                                                                                                                              >  2:6314988/1‑150 (MQ=255)

                          atgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑ccccc                                                                                                                             <  1:3494656/150‑1 (MQ=255)

                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCgccgccg‑cctgcg                                                                                                                            <  2:5371163/150‑5 (MQ=255)

                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGCGCGCGGG‑GGGGGGGGGC‑‑CAAAC‑CCCCCG‑ccccccc                                                                                                                           >  1:2031824/1‑148 (MQ=255)

                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑GGGGGGGGGC‑‑CAA‑CGCCGccc‑ccccccc                                                                                                                           >  1:2405043/1‑148 (MQ=255)

                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGGC‑CGGGGGGGCC‑‑CAAACccccccc‑cccccc                                                                                                                            >  1:4104024/1‑149 (MQ=255)

                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGGGGG‑‑CCAGCGCccccc‑ccccccc                                                                                                                           >  2:10291673/1‑148 (MQ=255)

                             atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATGACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑CCGGGGGGGC‑‑CCAAAGCccccc‑cccccccc                                                                                                                          >  2:9185906/1‑147 (MQ=255)

                               atTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCAc                                                                                                                        <  2:3924750/150‑1 (MQ=255)

                                  ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCGCGGGGGGGGC‑‑CAAACccccccc‑cccccccccccc                                                                                                                      >  2:1892513/1‑143 (MQ=255)

                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑GGGGGGCGGC‑‑CCAACGCCCCCC‑CCCCCAGCcccccc                                                                                                                    >  1:1616732/1‑150 (MQ=255)

                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CCAACGCCACCG‑CCCCCACCAccccc                                                                                                                    >  1:747858/1‑150 (MQ=38)

                                   tAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCAccccc                                                                                                                    <  2:7866532/150‑1 (MQ=255)

                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCgg                                                                                                                  <  2:6856756/150‑1 (MQ=255)

                                     aaaTTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCgg                                                                                                                  <  2:5490235/150‑1 (MQ=255)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcgg                                                                                                                                                  >  2:10259072/1‑120 (MQ=255)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcgg                                                                                                                                                  <  1:10259072/120‑1 (MQ=255)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCggg                                                                                                                 <  2:1627735/150‑1 (MQ=255)

                                            ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGTCGCAAGGAGCGC‑GA‑CCGGGGCGGCAACAAGAACCCCCC‑CCCCCCCCACCCCCGCGctcg                                                                                                             >  2:1338708/1‑149 (MQ=255)

                                             taatGTCATCAGGTCCGAAAATAACGCGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:542468/150‑1 (MQ=255)

                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGc                                                                                                          <  2:10609283/150‑1 (MQ=255)

                                                   catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACaa                                                                                                    <  2:8373629/150‑1 (MQ=255)

                                                    atcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCGCCCCCCTGCACCCCCGGGCTCTGGCGAACaa                                                                                                    >  2:8031965/1‑150 (MQ=255)

                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑CCGGGGCGGC‑‑CCAACCCCCCCC‑CCCCCCCCCCCCCCCCGCCCAGGCCAACcaaaaaacc                                                                                             >  2:5632115/1‑145 (MQ=255)

                                                            cGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑CGGGGGCGGC‑‑CCAAACCCGCCC‑CCCCCCCCCCCCCCCCCCGCCCGCGAACaaaaaaccccc                                                                                           >  2:1441018/1‑143 (MQ=255)

                                                             gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGggcggc                                                                                                                                                 <  1:4069623/98‑1 (MQ=255)

                                                                         aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CCGGGGCGGC‑‑CAAACCCCCCCC‑CCCCCAGCACCCCCCCGCGCAGGCGCACAAAAAAGCCCCCGCGCGGGGccca                                                                              >  2:1704210/1‑149 (MQ=21)

                                                                         aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCt                                                                              <  2:10678051/150‑1 (MQ=255)

                                                                          atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:1229264/150‑1 (MQ=255)

                                                                           tatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCg                                                                            <  2:9622863/150‑1 (MQ=255)

                                                                            atTTCATTCTCTCAGCCTGTTGCGCTCCTGTCTTGTGCATTGCTACATATAATCCCTGGCGCATGGAGCGC‑GC‑TGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  1:5838335/150‑1 (MQ=18)

                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:3104859/150‑1 (MQ=255)

                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:5519033/150‑1 (MQ=255)

                                                                            atTTCAGCCTCTCACCCGGTGGCGCTCCTGTCTTGTGCATTGCTCCATCTAAACACGGGGGCACGGCGCGC‑GC‑CGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  1:539728/150‑1 (MQ=17)

                                                                             tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  2:3816268/150‑1 (MQ=255)

                                                                             tttCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGc                                                                          <  2:10801566/150‑1 (MQ=255)

                                                                                     tctcATCCTGTTGCGCTCCTGTCTTGTGAATTGCTTAATATAATCCCTGGCGCTTGGTGCGG‑GC‑TGGGGGGGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCCAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCt                                                                  <  1:5262257/150‑1 (MQ=18)

                                                                                       tcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  2:10081864/150‑1 (MQ=255)

                                                                                              gTTGCGCTCCTGTCTTGGGCATTGCTTCATATAATCACTGGCGCATGGAGCGC‑GC‑CGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  1:4429191/150‑1 (MQ=255)

                                                                                               ttGCGCTCCTGTCTTGTGCATTGCTTCATATAATCACTGGCGCCCGGTGCGC‑GC‑TGGGGGCGGC‑‑CAAGCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAcc                                                        <  1:10173359/150‑1 (MQ=17)

                                                                                                 gcgcgcCGGTCTGGTGCTTTGCTTCATTTACTCAGGGGCGCATGGTGCGC‑GC‑TGGGGGGGGC‑‑CACTCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:933867/150‑1 (MQ=14)

                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCgg                                                      <  1:6191626/150‑1 (MQ=255)

                                                                                                    cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:8739612/150‑1 (MQ=255)

                                                                                                     tCCTGTCATGTGCATTGCTTCATCTAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGTGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCg                                                  <  2:5859234/150‑1 (MQ=255)

                                                                                                        tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGCGGC‑‑CCAACGCCCCCC‑CCCCCAGCACCCCCCGGCGCGGGCGAACAAAATAGCCCCGGCGCGGGGCCCACGGCGTATACCATACGGCGACCCGGTCCGgc                                               >  2:851324/1‑150 (MQ=255)

                                                                                                         gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑TGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                              <  1:1580131/150‑1 (MQ=255)

                                                                                                           cATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGcgcg                                            <  2:8804859/150‑1 (MQ=255)

                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGCGGC‑‑CCATCGCCGCCG‑CCCCCAGCACCCCCGGGCTCTGGCGAACAAAATCGCCCCTGCGCGGTGCCCACGGCGTATCCCATACGGCTAACGGGTCGGGCGCGa                                           >  2:7231441/1‑150 (MQ=17)

                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:5800870/150‑1 (MQ=255)

                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAAGCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:5632115/150‑1 (MQ=255)

                                                                                                            aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAACCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  1:6506544/150‑1 (MQ=255)

                                                                                                             tgtgCATTGCTTCATATAATCCCTGGCGCAAGGCGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  1:7810051/150‑1 (MQ=255)

                                                                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGCGGC‑‑CAAAACCCCCCC‑CCCCCCGCACCCCCCGGCGCTGGCGAAAAAAAAACCCGCGGCGCCGTGCCCACCGCGTATCCCATAACGCAACCCGGTCGGGCGCGGGGt                                        >  2:2623375/1‑150 (MQ=255)

                                                                                                                gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                       <  1:9902090/150‑1 (MQ=255)

                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGGGGC‑‑CCAACGCCCCCC‑CCCCCCCCACCCCCCCGCGCTGGCGCAAAAAAAACCCCCCGCGCCGTGGCCGCCGCGGATCCCGAAAGGCGACCGGGgcgggggcgggggacca                                    >  2:9664909/1‑145 (MQ=255)

                                                                                                                    tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GG‑CGGGGGGGGC‑‑CCCGAGGCCCCC‑CCCCCATCACCCCCCCCCGCTGGCGAACAAAATCCCCCGAGCGCCGGGCCCGCCGCGTCTCCCGTCCGGCTACCCGGTCGGGGGCGCggggacaa                                   >  2:7563847/1‑144 (MQ=255)

                                                                                                                    tGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGCGGC‑‑CAAAAGCCGCCC‑CCCCCAGCACCCCCGGGCGCAGGCGAACAAAAAACCCGCGGCGCGGTGCCCGCCGCGTATACCAAACGGCGACCCGGGCGGGCGCGCgggcacaa                                   >  2:4574351/1‑144 (MQ=255)

                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:906317/150‑1 (MQ=39)

                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:7352379/150‑1 (MQ=39)

                                                                                                                          atatAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCCCCC‑CCCCCAGCACCCCCCGGCGCTGGCGAACAAAAACCCCGCTGCGCCGGGCCCGCCGCGTATACCGTAAGGCTACCCGGTCGGGCGCGAGGGAACATCCCAGt                             >  2:827116/1‑150 (MQ=11)

                                                                                                                            ataATCACTGGCGCAAGGAGCGC‑GC‑CGGGGGCGGC‑‑CCAAAGCCCCCC‑CCCCCCGCACCCCCCGGCGCGGGCGCCCAAAAAAGCCGCGGCGCGGTGCCCGCCGCGTATCCCAGACCGCTACCCGGGCGGGCGCGCGGGAACAAcccaggaa                           >  2:6138172/1‑145 (MQ=255)

                                                                                                                                  aCTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:4330972/150‑1 (MQ=32)

                                                                                                                                   cTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTACACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTGATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    >  2:198189/1‑150 (MQ=17)

                                                                                                                                        gcAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCa               <  2:3409976/150‑1 (MQ=32)

                                                                                                                                           aGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:2283507/150‑1 (MQ=32)

                                                                                                                                             gAGCGC‑GC‑AGGGGGCGGC‑‑CACTCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:7220060/150‑1 (MQ=17)

                                                                                                                                                 gc‑gc‑AGGGGGCGGC‑‑CCATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCGCGGCGGCTCCCTTACGGCGACCGGGTCGGGCGCGCGGTCACAGCCCGGGAAAACGCGCCCGCCGcccccctc      >  2:2908654/1‑145 (MQ=11)

                                                                                                                                                       aGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:2361547/1‑150 (MQ=32)

                                                                                                                                                       aGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:3165417/1‑150 (MQ=32)

                                                                                                                                                       aGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  <  1:9785148/150‑1 (MQ=32)

                                                                                                                                                       aGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCAGCACCCCCGGGCACTGGCGAACAAAATAGCCGCAGCGCGGTGCCCGCGGCTTATACCATACGGCTACCGGGTAGGGCGCGAGGTAACATCCCTGTAAAAAGCGCCCACAGCCCACAACCaat  >  1:287990/1‑148 (MQ=11)

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GC‑AGGGGGCGGC‑‑CAATCGCCGCCG‑CCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑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: `‑`