BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,229,177	A→C	8.8%	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	8.8%	599.4 / 6.6	287	intergenic (+112/+261)	yohP/dusC	uncharacterized protein/tRNA‑dihydrouridine synthase C
Reads supporting (aligned to +/- strand): ref base A (51/167); new base C (22/9); total (94/193)
Fisher's exact test for biased strand distribution p-value = 3.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 291 total aligned reads displayed.

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326

  tggtggTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑ggg                                                                                                                                                      <  2:6399559/150‑1 (MQ=255)

    gtggtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑Cg‑ggggg                                                                                                                                                    >  1:5266612/1‑150 (MQ=255)

       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGggc‑gg                                                                                                                                                <  1:7908418/150‑1 (MQ=255)

       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGggc‑gg                                                                                                                                                <  2:6908876/150‑1 (MQ=255)

       gtgACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGggc‑gg                                                                                                                                                <  1:2863628/150‑1 (MQ=255)

         gACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGggc‑gggc                                                                                                                                              >  2:1120895/1‑148 (MQ=255)

          aCTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCCa                                                                                                                                             >  2:1514918/1‑150 (MQ=255)

            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CG‑GGGGGG‑GGccca                                                                                                                                            >  2:7807960/1‑148 (MQ=255)

            tGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GC‑GGGGGG‑GGcccaa                                                                                                                                           >  2:4477426/1‑147 (MQ=255)

             ggagTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATc                                                                                                                                          <  1:4151959/147‑1 (MQ=255)

             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCCGCCTGTTGCGCTCCTGTCATGGGCATTGCTGCATATACTCACTGGCGCAAGGAGCGC‑G‑CG‑gggggg‑gggggccc                                                                                                                                          >  1:2624352/1‑144 (MQ=255)

             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATc                                                                                                                                          <  2:7571445/150‑1 (MQ=255)

             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATc                                                                                                                                          <  2:5027420/150‑1 (MQ=255)

             ggCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATc                                                                                                                                          <  1:1325096/150‑1 (MQ=255)

              gCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATcg                                                                                                                                         >  1:4246991/1‑150 (MQ=255)

                 tATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATcgccg                                                                                                                                      <  2:3684106/150‑1 (MQ=255)

                   tttAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCACTagccgcc                                                                                                                                    >  2:750504/1‑150 (MQ=255)

                      aaGATGATATTTTAAAATAAATTAATGTCATCAGGTCCGAAAATAACTAGAATATTTAAGTCTCTCATCATGTTGCGCTCCTGTCATGTGCAATGCTTCAGATACTCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCCAAACgccgccccc                                                                                                                                 >  1:6629440/1‑150 (MQ=255)

                       aGATGATATTTTAAAATTAATTCATGGCATCAGGTCCGAAAATAACGAGAATATTACAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGccc                                                                                                                                <  1:9102740/150‑1 (MQ=255)

                        gatgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGcccc                                                                                                                               <  1:9223225/150‑1 (MQ=255)

                          tgatATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCACCTGTCATGTGCATTGCTTCATAAAATCACTGGCGCAAGGAGCGC‑G‑GA‑GCGGGCGGGCCAACCACcccccccccc                                                                                                                              >  2:4885707/1‑150 (MQ=255)

                           tatatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTg                                                                                                                            <  1:594500/149‑1 (MQ=255)

                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGCGG‑CG‑GGGGGG‑GGCCCA‑CGCccccccccccccc                                                                                                                           >  1:2124092/1‑147 (MQ=255)

                            atatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCCCATCCCCCCCCCCCCCCGc                                                                                                                           >  2:9128263/1‑150 (MQ=38)

                            atatTTTAAAATTAATGAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGc                                                                                                                           <  2:6953545/150‑1 (MQ=255)

                             tatTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCa                                                                                                                          <  1:8387184/150‑1 (MQ=255)

                                tttAAAATTAATTAATGTCATCAGGTCCGATAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCAccc                                                                                                                       <  2:3852495/150‑1 (MQ=255)

                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCTTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCACCCCCTGCAccc                                                                                                                       >  1:10087770/1‑150 (MQ=255)

                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GG‑GGGGGG‑GGGCAAAAACCCCCCCCCCCCCCAccc                                                                                                                       >  1:333236/1‑150 (MQ=255)

                                tttAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCACCCCCAGCAccc                                                                                                                       >  1:302191/1‑150 (MQ=255)

                                     aattaattAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCggg                                                                                                                  <  2:223538/150‑1 (MQ=255)

                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGG‑GGCCCCTAGCCGCCCCCCCCCCCCCCCCCCGGcg                                                                                                                >  2:432490/1‑149 (MQ=255)

                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                <  2:4092700/150‑1 (MQ=255)

                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGct                                                                                                                <  2:7944853/150‑1 (MQ=255)

                                       ttaattaaTGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCAGCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCCAATAGccgcccccccccgcccccccccgcg                                                                                                                >  1:6851257/1‑145 (MQ=255)

                                         aattaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑GCCG‑GGGGGG‑GGCCAAACCACCCCCCCCCCCCCAcccccccccca                                                                                                               >  2:2904461/1‑144 (MQ=255)

                                          attattGTCATCTGGTCCTAAAATAACTAGAATTTTTCAGTCTCTCTTCCTGTTGCGCTCCTTTCATGTGAATTGATTAATATAATCACTGGCGCAAGGAGGGC‑G‑CT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTg                                                                                                             <  2:2610371/150‑1 (MQ=21)

                                           ttaatGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GA‑GGGGGC‑GGCCCAAAGCCCCCGCCCCCACCCCCCCCCGGcacagg                                                                                                            >  1:1712356/1‑145 (MQ=255)

                                            taatGTCATCTGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGc                                                                                                           <  2:9662687/150‑1 (MQ=255)

                                              atGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑g                                                                                                                                                        >  1:3685400/1‑104 (MQ=255)

                                                     caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CG‑GGGGGC‑GGCCCAAAGCCGCCGCCCCCCCCACCCCCCGGCACAGGCG‑CACaaaaa                                                                                                 >  2:6352831/1‑149 (MQ=25)

                                                     caGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCCAACGCCGCCGCCCCCAGCACCCCCCGGCGCGGGCG‑GACaaaaa                                                                                                 >  2:6352938/1‑149 (MQ=25)

                                                        gTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GG‑GGGGGG‑GGGCCCGCCCCCCCCCCCCCCCCACCCCCCCCCCCGGGCG‑AGCaaaaaacc                                                                                              >  1:7982163/1‑146 (MQ=255)

                                                          ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑GG‑GGGGGG‑GGGCCAAAGCCCCCCCCCCCCCCCCCCCCCCCCCCAGGCG‑ACCaaaaacaccc                                                                                            >  1:7226858/1‑144 (MQ=255)

                                                          ccGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATcgccg                                                                                            <  2:1653011/150‑1 (MQ=255)

                                                            gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCCATCGCCGCCCCCCCCAGCCCCCCCCCGCGCTGGCGCAACaaaacgccccg                                                                                           >  2:10124606/1‑143 (MQ=25)

                                                            gAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCCCCCCCGGGCTCTGGCG‑AACAAAAACGCCGCt                                                                                          >  2:8478222/1‑150 (MQ=255)

                                                                   aCGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGt                                                                                   <  2:6695654/150‑1 (MQ=255)

                                                                     gagaATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCAGCAACCCCCGGCTATGGCG‑AACAAAAAAGCCGCGGCGCGGTGc                                                                                 >  2:3436666/1‑150 (MQ=21)

                                                                         ctatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:6720933/149‑1 (MQ=255)

                                                                         atatTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCAGTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTc                                                                             <  2:5581928/150‑1 (MQ=255)

                                                                           atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTCCATATAATCACTGGCGCACGGAGCGC‑G‑CC‑GGGGGC‑GGCCAAGCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:8796637/150‑1 (MQ=255)

                                                                           atTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATATGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCgg                                                                           <  2:6005747/150‑1 (MQ=255)

                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  1:6308374/150‑1 (MQ=255)

                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  2:4219207/150‑1 (MQ=255)

                                                                               cAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTa                                                                       <  2:5242443/150‑1 (MQ=255)

                                                                                 gTCTCTCATCGTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAAACACTGGCGAATGGAGCGC‑G‑CA‑GGGGGC‑GGCAAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  1:10225172/150‑1 (MQ=18)

                                                                                 gTCTCTCAGCCGGGTGCGCTCCGGTCATGTGCATTGCTTCAGCTAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCCGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  2:10024429/150‑1 (MQ=17)

                                                                                 gTCTCTAATCCTGTTGCGCCCCTCTCATGTGAATTGCTTAATATAAACACTGGCCCAAGGAGCGC‑G‑CT‑GGGGGC‑GCCAAATCGCCGCCCCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATc                                                                     <  1:4765437/150‑1 (MQ=17)

                                                                                     ctcATCCTGTGGCGCTCCTGTCTTGTGAATTGCTTCATATAATCACTGGCGCACGGAGCGC‑G‑CT‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCtt                                                                 <  1:9561980/150‑1 (MQ=18)

                                                                                      tcATCCTGCTGCGCTCCTGTCATGTGCATTGCTTCATATAATCAGTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTa                                                                <  1:11463/150‑1 (MQ=255)

                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:8344443/150‑1 (MQ=255)

                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:7380154/150‑1 (MQ=255)

                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:1721545/150‑1 (MQ=255)

                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  2:6890879/150‑1 (MQ=255)

                                                                                        aTCCTGTTGCGCTCCTGTCATGTGCATTGCTACATATAAGCAGTGGCGCACGGCGCGC‑G‑CT‑GGGGGC‑GCCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACg                                                              <  1:4808619/150‑1 (MQ=17)

                                                                                            tGTTGCGCTCCTGTCATGTGCATTGCTTAATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTa                                                          <  2:1962384/150‑1 (MQ=255)

                                                                                             gTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTAc                                                         <  2:7129948/150‑1 (MQ=255)

                                                                                               tGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCg                                                       <  2:1307880/150‑1 (MQ=255)

                                                                                                 cgcTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CC‑GGGGGC‑GGCCAAACGCCCCCCCCCCCAGCACCCCCCGGCGCAGGCG‑CACAAAAACGCCCCGGCGCGGGGCCCCCGGCGTAACCCCTACCGCTAccccg                                                     >  2:6511822/1‑147 (MQ=17)

                                                                                                   cTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTc                                                   <  2:5768213/150‑1 (MQ=255)

                                                                                                     ccTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCgg                                                 <  1:1514918/150‑1 (MQ=255)

                                                                                                      cTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCggg                                                <  2:1823515/150‑1 (MQ=255)

                                                                                                       tGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCGGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgc                                               <  1:8814541/150‑1 (MQ=255)

                                                                                                        gTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGgcg                                              <  1:2142951/150‑1 (MQ=255)

                                                                                                           cggtgCAGTGCGTCAGAGACTCACTGGCGCCCGGAGCGC‑G‑CC‑GGGGGC‑GGCCACTCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑ACCAAAATCGCCGCTGCGCGGGGCCCTCGGCTTCTCCCTTACGGCGACCGGGTCGGGCGCGa                                           <  2:8907746/148‑1 (MQ=11)

                                                                                                           aTGTGCATTGCTTCATATACTCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:9449583/150‑1 (MQ=255)

                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑TGGGGC‑GGCCAAACGCCGCCGCCCCCAGCCCCCCCCGGCTATGGCG‑AACAAAAACCCCGCGGCGCGGTGGCCACGGCATAACCCATAACGCTACCCGGTCGGGcgcgc                                           >  2:2504346/1‑149 (MQ=255)

                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGGCCAACGCCGCCCCCCCCCGCACCCCCCGGCGCAGGGG‑GCCAAAAAAGCCGCGGCGCGGGGCCCACGGGGGATCCCCGAAGGCGACCCGGGGGGGcgcgc                                           >  2:9862290/1‑149 (MQ=12)

                                                                                                           aTGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGa                                           <  2:972547/150‑1 (MQ=255)

                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGG‑GGCCAAACGCCGCCCCCCCCAGCACCCCCGGGCTCTGGCG‑CACAAAATCCCCGCTGCGCCGTGCCCTCGGCTTATCCCATAACGCGACCCGGTCGGGcgcgcg                                          >  2:3289100/1‑148 (MQ=17)

                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑ACCAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:9051928/150‑1 (MQ=255)

                                                                                                            tgtgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAg                                          <  2:191608/150‑1 (MQ=255)

                                                                                                              tgCATTGCTTCATATAATCAGTGGCGAAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                        <  1:4728962/150‑1 (MQ=255)

                                                                                                              tgCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGt                                        <  2:6795639/150‑1 (MQ=255)

                                                                                                               gCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTa                                       <  2:5479274/150‑1 (MQ=255)

                                                                                                                cATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTaa                                      <  2:6634453/150‑1 (MQ=255)

                                                                                                                  ttGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAAACGCCCCCCCCCCCAGCACCCCCCGGCTCAGGCG‑AAAAAAAAAGCCCCGGCGCGGTGCCCACGGCGTGTACCGTAACGCAAACCGGTCGGGCGCGCGGGAACa                                    >  2:10101488/1‑150 (MQ=9)

                                                                                                                     cTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATcc                                 <  1:5672492/150‑1 (MQ=40)

                                                                                                                      ttCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCCCCCCCAGCACCCCCCGGCTCTGGCG‑AACAAAAAACCCGCGGCCCGGTGCCCACGGCGTATCCCATAAGGCTAACCGGTCGGGCGCGAGGTAACATacc                                >  2:10289735/1‑150 (MQ=14)

                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  2:7248947/150‑1 (MQ=39)

                                                                                                                       tCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCt                               <  1:1855658/150‑1 (MQ=39)

                                                                                                                                    ggCGCCTGGAGCGC‑G‑CC‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGccc                  <  1:8626612/150‑1 (MQ=12)

                                                                                                                                      cgcAAGGAGCGC‑G‑CAGGGGGGC‑GCCCAAAC‑CCCCCCCCCCCACCACCCCCCGGCGCAGGCG‑AACAAAAAAGCCCCGGCGCGGGGGCCCCCGCGGAACACAAAAGGCGACCCGGGCGGGCGCGAGGGAACAAACCAGTAAAACCCGccccc                >  2:3232206/1‑148 (MQ=2)

                                                                                                                                          aGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGcc            <  2:4663737/150‑1 (MQ=32)

                                                                                                                                           ggAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCACGGGCTCTGGCG‑AACAAAATCGCCGCTTCGCTGTGCCCTCGGCTTATCCCTTACTGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGccc           <  1:4954180/150‑1 (MQ=255)

                                                                                                                                             aGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCcac         <  1:974349/150‑1 (MQ=32)

                                                                                                                                                 c‑g‑cA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATcc     <  1:5454848/150‑1 (MQ=32)

                                                                                                                                                      a‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:4850231/1‑150 (MQ=32)

                                                                                                                                                      a‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCAtg  >  2:52904/1‑150 (MQ=32)

ACTGGTGGTGACTGGCGTATTTAAGATGATATTTTAAAATTAATTAATGTCATCAGGTCCGAAAATAACGAGAATATTTCAGTCTCTCATCCTGTTGCGCTCCTGTCATGTGCATTGCTTCATATAATCACTGGCGCAAGGAGCGC‑G‑CA‑GGGGGC‑GGCCAATCGCCGCCGCCCCCTGCACCCCCGGGCTCTGGCG‑AACAAAATCGCCGCTGCGCGGTGCCCTCGGCTTATCCCTTACGGCTACCGGGTCGGGCGCGAGGTAACATCCCTGTAAAACGCGCCCTCAGCCCACATCCATG  >  NC_000913/2229029‑2229326

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