BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,287,122	T→G	5.6%	noncoding (115/171 nt) R19R (AGG→CGG)	rttR ← tpr ←	rtT sRNA, processed from tyrT transcript protamine‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,122	0	T	G	5.6%	734.1 / 11.7	322	noncoding (115/171 nt) R19R (AGG→CGG)	rttR tpr	rtT sRNA, processed from tyrT transcript protamine‑like protein
Reads supporting (aligned to +/- strand): ref base T (146/133); new base G (22/2); total (180/142)
Fisher's exact test for biased strand distribution p-value = 1.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

 ccGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTg                                                                                                                                                      >  1:4648541/1‑150 (MQ=255)

   gAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                                                                                    <  1:612092/150‑1 (MQ=255)

       cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCaa                                                                                                                                                >  1:6338235/1‑150 (MQ=255)

       cTTAGTCGAAGCTTCTCATACTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCggcggcca                                                                                                                                                >  1:2693235/1‑148 (MQ=255)

        ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCGGAGGCAAc                                                                                                                                               >  1:2085298/1‑150 (MQ=255)

        ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCGGCGCGCTTCGCCCTGCGGGTCGTTGCCGGCGGCAAc                                                                                                                                               >  2:1413111/1‑150 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATTGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctc                                                                                                                                           >  2:5638653/1‑150 (MQ=255)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGCGCTGCGCGCTTCGCCCATCGGGTCGTTTCCTGCGGGAAcccgc                                                                                                                                           >  2:355789/1‑146 (MQ=14)

            tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGGGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGcga                                                                                                                                           >  2:4603017/1‑148 (MQ=255)

             cGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctct                                                                                                                                          <  1:6552533/150‑1 (MQ=255)

              gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGATTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGTAGGAGTAGGCTGCGCTGCGGGCTGCGCGCTGGGGGTCGTGTCGTGCGGCGACGcgcgc                                                                                                                                         >  2:672770/1‑146 (MQ=14)

                  cTTCTCATCCTTCCCCGCATGGGCAGAATATTTGAGTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGc                                                                                                                                     <  2:8585682/150‑1 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTGgc                                                                                                                                 >  2:5463237/1‑150 (MQ=255)

                      tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCGTCGGGTCGTTGCCTGCCGCAAAGCTCTCTCGCTGgc                                                                                                                                 >  2:6156714/1‑150 (MQ=255)

                       cATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCAGGGACTTTTGAAAGTGATTGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGgcg                                                                                                                                >  1:5659808/1‑150 (MQ=255)

                       cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACACTCGCTGgcg                                                                                                                                >  1:5259017/1‑150 (MQ=255)

                       cATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTGGCTGCGGCAACGCTCTCACGCTGgcg                                                                                                                                >  2:9084854/1‑150 (MQ=255)

                         tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt                                                                                                                                    <  1:5386379/144‑1 (MQ=255)

                          ccTGCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc                                                                                                                             <  2:6773742/150‑1 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCATCGCTGCGCGCATCGCCATACGTGTCGTTGCCTGAGGCAAAGCACTCTCGCTGGCGCTCGa                                                                                                                           >  2:8409324/1‑150 (MQ=14)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           >  1:8861092/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACGCTCGCTGGCGCTCGa                                                                                                                           >  1:7553462/1‑150 (MQ=255)

                            ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCGCGCTCGCTGGCGCTCGa                                                                                                                           >  2:4933131/1‑150 (MQ=255)

                            gtCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa                                                                                                                           <  2:1694806/149‑1 (MQ=255)

                             tCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGATTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAg                                                                                                                          <  1:4138740/150‑1 (MQ=255)

                              ccccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                         <  1:1688416/150‑1 (MQ=255)

                               cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTc                                                                                                                        >  2:8359306/1‑150 (MQ=255)

                               cccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAAAGCTATCTCGCTGGCGCTCGAgac                                                                                                                        >  1:9505987/1‑148 (MQ=255)

                                ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg                                                                                                                       >  2:5509322/1‑150 (MQ=255)

                                ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg                                                                                                                       >  2:7148792/1‑150 (MQ=255)

                                 cGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGa                                                                                                                      >  1:6050898/1‑150 (MQ=255)

                                    aTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCGCTCACGCGGGCGCGCGAGTCGAAcc                                                                                                                   >  1:8257478/1‑150 (MQ=14)

                                       ggCAGATTCTTTGATTGCGGATTCGTTTGAGTTTTCGTGTTCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTg                                                                                                                <  2:5935659/150‑1 (MQ=11)

                                         cAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGt                                                                                                              >  2:1957858/1‑150 (MQ=255)

                                          aGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAATCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAAAGCTCTCTCGCGGGCGCTCGAGTCGACCCAAGGTc                                                                                                             >  1:2788120/1‑150 (MQ=14)

                                             atatTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                          <  1:1424277/150‑1 (MQ=255)

                                                 ttGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCtt                                                                                                      <  1:4394352/150‑1 (MQ=255)

                                                  tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTtc                                                                                                     >  2:3865757/1‑150 (MQ=255)

                                                  tGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTTCCTGCGGCAAACCACACTCGCTGGCGCACGAGTCGAACCAGTGTCGAAGAtta                                                                                                     >  1:2268803/1‑149 (MQ=14)

                                                   gATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGAGGCAAAGCTCTCACGCGGGCGCACCAGTCGAACCTTGGTAGAAGCttat                                                                                                    >  1:1224278/1‑148 (MQ=14)

                                                   gATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGaa                                                                                                                     >  2:5694670/1‑133 (MQ=255)

                                                                cTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCACACACGCAGGCGCGAGAGACGAACCTTGGTCGAATATTCTCATCCTACCCCGc                                                                                       >  1:5210993/1‑150 (MQ=14)

                                                                 ttGATAGTTCATGGTGTTTTGATAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                      <  1:6733569/150‑1 (MQ=14)

                                                                      agTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCATCCCCGCATGGGc                                                                                 >  2:1301194/1‑150 (MQ=255)

                                                                           aGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACCCAAACACGCTGGCGCTCGAGGAGAACCGTGGTCGAAGCATATCAACCATACCCCCAAGGGAAGAat                                                                            >  1:9161492/1‑150 (MQ=11)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCGAACCATGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAatat                                                                          >  2:555080/1‑150 (MQ=255)

                                                                             ggACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCAAGGGCAGAatat                                                                          >  2:6283428/1‑150 (MQ=255)

                                                                               aCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTGCGACCTTCGCCCTTCGTGTACATGCCTGCCGCAAAGCACTCTCGCTGGCGCACCAGTCGAAACTTGGTAGAAGCTTCTAAAACTTCACCGCATGGGCAGACGAttt                                                                        >  2:3738157/1‑150 (MQ=2)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCAACCTTCCCCGCAAGGTCAGaagagcgac                                                                      >  2:8712912/1‑143 (MQ=11)

                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCAAAACCGGGCCCTCGAGTCGAACCGTGGTCCAACCTTAAAAACCATCCCCCCAAGGGCAGCATATTTaa                                                                      >  1:8025986/1‑150 (MQ=2)

                                                                                   ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAATCTTATCAACCATCCCCCCAAGGGCAGAAAATTTaaat                                                                    >  2:8262271/1‑148 (MQ=11)

                                                                                   ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                    <  1:7047419/150‑1 (MQ=39)

                                                                                   ttGAAAGTGAATGTGGTGGGGGAAGGATTAGTCCGAGCTGCGCGCCACGCCCCTGGGGTCGTTGGCCGCTGTTACGCTCTCTCGCTGCCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAAtt                                                                    <  2:6877117/150‑1 (MQ=11)

                                                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTTGCGCTCGAGTCCAACCTTGGTCGAAGCGTCTCATCCTTCCCCGCAAGGGCAGAATATTTAATTGc                                                                  >  2:2353438/1‑150 (MQ=255)

                                                                                           gATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCGCGAGTCGAACCTTTGTCGCAGCTTATAATACGTCCCCCCAGGGGCCGCAGATTTAAATTCggcggc                                                            >  1:7134469/1‑146 (MQ=11)

                                                                                             tggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGt                                                          >  2:9553216/1‑150 (MQ=37)

                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt                                                                                                                         <  1:6819808/82‑1 (MQ=33)

                                                                                                      gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCATACCCGCAAGGGCAGAATATTTAATTGCGGATTCGTTGGGTAGtt                                                 >  2:7664918/1‑150 (MQ=255)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                          >  2:8863052/1‑92 (MQ=33)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGaa                                                                                                          <  1:8863052/92‑1 (MQ=33)

                                                                                                       ggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCATACCCGCATGGGCAGAATATTTAATTGCGGATTAGTTGGGAAGTTc                                                >  2:7775467/1‑150 (MQ=255)

                                                                                                          aGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAgg                                             <  1:531660/150‑1 (MQ=255)

                                                                                                            gATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAgc                                                                                                                         >  2:5169599/1‑71 (MQ=17)

                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAc                                          <  2:3337837/150‑1 (MQ=255)

                                                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACt                                         <  2:531646/150‑1 (MQ=255)

                                                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGACTCGAACCGTGGTCGAAGCTTCTCAACCATCCCCGCAGGGGGCGAATATTTCAGAGCGGAAGCGTTGGGACGGTCCGGGCCt                                         >  2:4460041/1‑150 (MQ=11)

                                                                                                               tACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                        <  1:7309871/150‑1 (MQ=255)

                                                                                                               tACTCAGCGCGGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACtt                                        <  1:3821250/150‑1 (MQ=255)

                                                                                                                aCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACttt                                       <  1:2431181/150‑1 (MQ=255)

                                                                                                                  tCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGTTGGGGAGTTCGGGGACGTTTg                                     >  2:1081013/1‑150 (MQ=14)

                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   >  2:6102010/1‑150 (MQ=255)

                                                                                                                    aGCGCCGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCGGGCGCTCGAGTCGACCCCTGGTCGAAGCTCCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaa                                   <  2:5734663/150‑1 (MQ=17)

                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGaaa                                  <  1:1876511/150‑1 (MQ=255)

                                                                                                                      cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 <  2:9505987/150‑1 (MQ=255)

                                                                                                                      cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 <  2:2322112/150‑1 (MQ=255)

                                                                                                                      cgcAGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGACGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAg                                 <  2:2145251/150‑1 (MQ=255)

                                                                                                                        cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTg                               <  1:4910336/150‑1 (MQ=255)

                                                                                                                             cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggt                          <  2:6414147/150‑1 (MQ=255)

                                                                                                                              gcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCGTGGGCGAAGCGTCTCATCCGTCCCCGCCTGGGCCGAATCGTTAATTGCGGCTTCGGGGGGACGGTCAGGGCCTGTGGACAGTGAtgggg                         >  2:439017/1‑148 (MQ=14)

                                                                                                                                 tCGCCGTTCGGTTCGTTGCCCGCGTCTTGGGTCTGTAGTTGTGTCTCGAGTAGTTCCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtggtggtg                      <  1:780548/150‑1 (MQ=14)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   <  1:498962/150‑1 (MQ=255)

                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTATCTCGCTGGCGCTCGAGTCGAACCGTGGTCGAAGCGTATCATCCTTCCCCCCAAGGGCAGAAAATTTAAATGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTgggg                   >  2:2903700/1‑150 (MQ=17)

                                                                                                                                     ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAAAATTAAAATGCTGAGCCGTTGTGACCGTGAGGGACGGTTGAAAGTGTtgggggtggggg                  >  2:791759/1‑150 (MQ=14)

                                                                                                                                     ccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAAATTTAAATGCGGAACCGTTGGGAAGGTAAGGGCCTTTTGAAAGTGATTGGGGTggggg                  >  2:8086074/1‑150 (MQ=17)

                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  1:132801/150‑1 (MQ=255)

                                                                                                                                       ttCGGGACGTTGCCAGCGGCCACGCTCTCTCGCTGGCGCTCGCGCCGAACCTTGTTCGAAGCTTCTCTTACTTCCCCGCTTGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGaa                <  2:6466812/150‑1 (MQ=14)

                                                                                                                                          gggTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa             <  1:8359306/150‑1 (MQ=255)

                                                                                                                                          gggTCGTGGCCTGCGGCACCGCTCTCTCGCGGGCGCTCGCGTCGAACCGTGGTCGAAGCTTCTCCTCCGTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGa             <  1:6622496/150‑1 (MQ=17)

                                                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTAGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           <  1:5779332/150‑1 (MQ=255)

                                                                                                                                            agcGTTGCCGCAGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAtt           >  2:2545874/3‑150 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          >  2:751483/1‑150 (MQ=255)

                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTc          <  2:2217402/150‑1 (MQ=255)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                   <  1:4832948/91‑1 (MQ=37)

                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTg                                                                   >  2:4833344/1‑91 (MQ=18)

                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGaa       <  1:891203/150‑1 (MQ=255)

                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                        <  1:8610941/84‑1 (MQ=38)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                        <  1:9373247/81‑1 (MQ=38)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                        >  2:9373247/1‑81 (MQ=38)

                                                                                                                                                    cTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAACCTTACCCCCAAGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGAATTTTGAAAGAGATGGGGGTGGGGGAAGGATTAGAACCtt   >  1:8749346/1‑150 (MQ=17)

                                                                                                                                                     tGCGGCAACGCTCTCTCGCTGGCGCTCAAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGAtg                            >  2:1787240/1‑124 (MQ=255)

GCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTC  >  NC_000913/1286973‑1287271

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