BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	14.5%	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.5%	431.9 / 12.0	230	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base A (44/119); new base C (27/7); total (81/149)
Fisher's exact test for biased strand distribution p-value = 1.61e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 240 total aligned reads displayed.

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229028‑2229326

tACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑Ca                                                                                                                                                        <  2:487972/150‑1 (MQ=255)

   tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAggg                                                                                                                                                     >  1:4603652/1‑150 (MQ=255)

    ggtggtGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAgggg                                                                                                                                                    >  2:3876444/1‑150 (MQ=255)

     gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAggggg                                                                                                                                                   >  2:4737827/1‑150 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑g‑gggggggg                                                                                                                                                  >  2:3077547/1‑149 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggc                                                                                                                                                  <  1:5171431/150‑1 (MQ=255)

      tggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggc                                                                                                                                                  <  1:5348841/150‑1 (MQ=255)

        gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTGCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGggcgg                                                                                                                                                >  1:4093600/1‑150 (MQ=255)

         tgACTGGCGTATTTAAGATGATATTTTAAAATTAATTACTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcggc                                                                                                                                               <  2:5248424/150‑1 (MQ=255)

            cTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GAGGGGGCGGGCaa                                                                                                                                            >  2:4550876/1‑150 (MQ=255)

             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCCGGGGGCGCCCaa                                                                                                                                            >  2:3422669/1‑150 (MQ=255)

             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGcccca                                                                                                                                           >  1:6498369/1‑147 (MQ=255)

             tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑gggggggggggccca                                                                                                                                           >  1:5620965/1‑145 (MQ=255)

              ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCaaac                                                                                                                                          >  1:5304638/1‑148 (MQ=255)

               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CCGGGGGCGCCCaacc                                                                                                                                          >  2:2781019/1‑148 (MQ=255)

               gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcg                                                                                                                                         <  2:1433159/150‑1 (MQ=255)

                 gTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgcc                                                                                                                                       <  2:6045054/150‑1 (MQ=255)

                  tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATcgccg                                                                                                                                      <  2:1930018/150‑1 (MQ=255)

                      tAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGCGGCGGCGATTggccgcccc                                                                                                                                  <  1:911522/150‑3 (MQ=255)

                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTTTTGCGCTCCTGTCTTGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCgccgccgcc                                                                                                                                 <  2:5868372/150‑1 (MQ=255)

                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCgccgccgcc                                                                                                                                 <  1:670352/150‑1 (MQ=255)

                       aaGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCgccgccgcc                                                                                                                                 <  1:4166023/150‑1 (MQ=255)

                         gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GGGGGGGCGGCCCATAAcccccccccc                                                                                                                               >  1:1012044/1‑150 (MQ=255)

                           tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                             <  2:2637409/150‑1 (MQ=255)

                           agatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCt                                                                                                                             <  2:2322045/149‑1 (MQ=255)

                            gatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAAAAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCATCATATAATCACTGGCGCAAGGAGCGC‑G‑CACGGGGCGGCCCAACccccccccccccag                                                                                                                            >  2:5172025/1‑148 (MQ=255)

                              tatTTTAAAATTAATTAATGTCATCAGGTCCGAAGATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCa                                                                                                                          <  2:5490273/150‑1 (MQ=255)

                                ttttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAAAAGCCCCCCCCCCCAGCAcc                                                                                                                        >  2:2509645/1‑150 (MQ=255)

                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCAccc                                                                                                                       <  2:619850/150‑1 (MQ=255)

                                 tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACGCCGCCCCCCCCAGCAccc                                                                                                                       >  2:5362617/1‑150 (MQ=255)

                                  ttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATAGCCGCCGCCCCCGGCAcccc                                                                                                                      >  2:4059460/1‑150 (MQ=255)

                                      aattaattCATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCACCGCCCCCTGCACCCCCggg                                                                                                                  <  2:6944371/150‑1 (MQ=255)

                                      aattaattAATGTCATCAGGTCCGACAATAACGAGACTATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCAGTGCTTCATATAATCAGTGGCGCACGGCGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:4437730/150‑1 (MQ=37)

                                      aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GGGGGGGGGGGCCCGAGCCCCCCCCCCCACCAcccccccg                                                                                                                  >  1:6065187/1‑147 (MQ=255)

                                         taattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGctc                                                                                                               <  2:6691713/150‑1 (MQ=255)

                                             taatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GCGGGGGCGGGCCAGCGGCCCCGCCCCCCGGACCCCCCGGCGGTGGc                                                                                                           >  2:4296728/1‑150 (MQ=37)

                                               atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATGCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GAGGGGGGGGGCTATCGGCGGCGCCCCCTGCACCCCCGGGGTCTGGCgg                                                                                                         >  2:5937062/1‑149 (MQ=255)

                                                tGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGaa                                                                                                        <  2:349834/150‑1 (MQ=255)

                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgg                                                                                                                                                <  2:5073809/108‑1 (MQ=255)

                                                  tcatcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGggcgg                                                                                                                                                >  1:5073809/1‑108 (MQ=255)

                                                   catcaGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGCCCAAAGCCGCCGCCCCCCGCCCCCCCCGGCGCTGGCGAAcca                                                                                                     >  1:6610618/1‑148 (MQ=255)

                                                        ggTCCGAACATAACGAGCATATGTCAGTCTCTCCTCCTGTTGCGCTCCTGTCATGTGCATTTCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CATGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATcg                                                                                                <  2:1772026/150‑1 (MQ=21)

                                                          tCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCCAACGCCCCCCCCCCCCCCACCCCCCCCCGCGCGCGAACaaaaaaccc                                                                                              >  2:3687607/1‑145 (MQ=25)

                                                              aaaaTAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTg                                                                                          <  2:3222614/150‑1 (MQ=255)

                                                                     cGATAATATTTCTGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGCTTTGCTTCATATAACCACGGGCGCATGGCGCGC‑G‑CTGGGGGCGGCCAACCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTg                                                                                   <  1:4016867/150‑1 (MQ=17)

                                                                      gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGc                                                                                  <  2:578461/150‑1 (MQ=255)

                                                                       agaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACGCCGCCCCCCCCAGCACCCCCCGGCACAGGCGAACAAAATACCCGCAGCCCGGGGcc                                                                                 >  2:7184892/1‑150 (MQ=21)

                                                                         aaTATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCt                                                                               <  2:3795713/150‑1 (MQ=255)

                                                                            atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCTTGTGAATTGCTTCATATAATCACTGGCGCAAGGCGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                            <  1:4406410/150‑1 (MQ=255)

                                                                                     tctcATCCTGTTGCGCTCCTGTCTTGTGATTTGCTACATATAACCACGGGCGCATGGTGCGC‑G‑CTGGGGGCGGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑t                                                                  <  1:895348/150‑1 (MQ=14)

                                                                                      ctcATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑tt                                                                 <  2:1667304/150‑1 (MQ=255)

                                                                                         aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAACGCCCCCCCCCCCAGCACCCCCCGGCTCAGGCCAACAAAAAAGCCGCTGCGCCGTGCCCACCGCGTATACC‑taaac                                                              >  2:7205556/1‑146 (MQ=17)

                                                                                            cTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAACCCCCCCCCCCCCCGCCCCCCCGGGCGCTGGCGCACAAAAAACCCCCGGCGCGGGGGCCACGGCGTAAACC‑CTACGGCt                                                           >  2:438846/1‑150 (MQ=12)

                                                                                               ttGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGGGGCCCAAACCCCCCCCCCCCGGCCCCCCCCGGCGATGGGGGACAAAATAGCCCCGGCGCGGTGGCCGCGGCGTTTCCCATTACgctaac                                                         >  2:5620613/1‑145 (MQ=255)

                                                                                                tGCGCTCCTGTCTTGTGCTTTGCTTAATATAACCCGTGGCGCACGTAGCGC‑G‑CGGGGGGCGCCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCg                                                       <  1:3077547/150‑1 (MQ=14)

                                                                                                tGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCg                                                       <  2:3448899/150‑1 (MQ=255)

                                                                                                 gcgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCgg                                                      <  2:2934980/150‑1 (MQ=255)

                                                                                                   gcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGt                                                    <  2:1020920/150‑1 (MQ=255)

                                                                                                     tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGGCCAAAGCCCCCCCCCCCCCCACCCCCCGGCGCAGGGGAAAAAAAACGCCCCGGCCCCGGGGCCACGGCGTATCCC‑CTACGGCGCCCCGGTCg                                                  >  2:3577740/1‑150 (MQ=12)

                                                                                                     tCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCg                                                  <  2:4782207/150‑1 (MQ=255)

                                                                                                     tCCTGTCACGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCAGCACCCCCGGGCTCTGGCGAAAAAAATCGCCGCTGCGCGGTGGCCTCGGCTTATCCC‑TTACGGCTACCGGGTCg                                                  >  2:6020160/1‑150 (MQ=255)

                                                                                                            aTGTGCATTGATTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAACACCGCCCCCCCCCTCAACCCCCCGCACTGGCGAAAAAAAAAGACGCGGCCCGGTGGCCGCGGCGTATACC‑AGAACTCTCCCAGGTCGGGCGCGa                                           >  2:420363/1‑150 (MQ=9)

                                                                                                             tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAg                                          <  1:4469747/150‑1 (MQ=255)

                                                                                                               tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCCAACGCCGCCCCCCCCAGCACCCCCCGGCTCAGGCGAACAAAAACGCCGCGGCGCGGTGCCCACGGCGTATCCC‑CTAACGCTACCCGGGCGGGCGCGCGGt                                        >  2:6623129/1‑150 (MQ=14)

                                                                                                                  aTTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAAc                                     <  1:3653574/150‑1 (MQ=255)

                                                                                                                   ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACa                                    <  2:6389781/150‑1 (MQ=255)

                                                                                                                     gCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATc                                  <  1:3393864/150‑1 (MQ=255)

                                                                                                                     gCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATc                                  <  1:3543023/150‑1 (MQ=255)

                                                                                                                      cTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:7110207/150‑1 (MQ=40)

                                                                                                                       ttCATATAATCCCTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCGGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATccc                                <  2:676627/150‑1 (MQ=255)

                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:6334838/150‑1 (MQ=255)

                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2852/150‑1 (MQ=39)

                                                                                                                        tCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:2990105/150‑1 (MQ=39)

                                                                                                                         cATATAATCACTGGCGCAAGGAGCGC‑G‑CCGGGGGCGGGCAAAAGACCCCCCCCCCCGCACCCCCCCGGTCAGGCGAAAAAAAAAGCCCCGGCCCGGGGCCCACCGCATATAAC‑ATAAGGCTAACCCGTCGGGCGCCCGGGAACaaaacag                              >  2:3864992/1‑144 (MQ=2)

                                                                                                                             taATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaa                          <  1:2275771/150‑1 (MQ=35)

                                                                                                                              aaTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:4571848/150‑1 (MQ=34)

                                                                                                                              aaTCACGGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTaaaa                         <  1:6135332/150‑1 (MQ=21)

                                                                                                                               aTCACTGGCGCATGGAGCGC‑G‑CTGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:2531542/150‑1 (MQ=18)

                                                                                                                               aTCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:5719223/150‑1 (MQ=33)

                                                                                                                               aTCACGGGCGCAAGGAGCGC‑G‑CCGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAc                        <  1:2497774/150‑1 (MQ=18)

                                                                                                                                 cACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgc                      <  2:4650925/150‑1 (MQ=32)

                                                                                                                                  aCTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAAcgcg                     <  1:5004903/150‑1 (MQ=17)

                                                                                                                                   cTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  1:7219586/150‑1 (MQ=32)

                                                                                                                                   cTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACgcgc                    <  2:2255015/150‑1 (MQ=32)

                                                                                                                                    tGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGcc                   <  2:1917551/150‑1 (MQ=32)

                                                                                                                                         cAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:3852149/150‑1 (MQ=32)

                                                                                                                                         cAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:7184970/150‑1 (MQ=32)

                                                                                                                                         cAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  2:1629031/150‑1 (MQ=32)

                                                                                                                                         cAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAg              <  1:1951581/150‑1 (MQ=32)

                                                                                                                                          aaGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGc             <  1:3763066/150‑1 (MQ=32)

                                                                                                                                           aGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  1:1303762/150‑1 (MQ=32)

                                                                                                                                             gAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCca          <  1:5858033/150‑1 (MQ=32)

                                                                                                                                              aGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac         <  1:1595643/150‑1 (MQ=32)

                                                                                                                                                  c‑g‑cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  2:3983672/150‑1 (MQ=32)

                                                                                                                                                  c‑g‑cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     >  2:6081059/1‑150 (MQ=32)

                                                                                                                                                  c‑g‑cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:5870515/150‑1 (MQ=32)

                                                                                                                                                      cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGTTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:5470575/150‑1 (MQ=17)

                                                                                                                                                      cAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAt   <  1:5985453/150‑1 (MQ=32)

                                                                                                                                                       aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:5466789/1‑150 (MQ=32)

                                                                                                                                                       aGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:2854123/1‑150 (MQ=32)

TACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCGAACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCC‑TTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  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: `‑`