BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,339,422	T→C	D172D (GAT→GAC)	yrbH →	putative isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,339,422	0	T	C	99.7%	1001.6 / NA	321	D172D (GAT→GAC)	yrbH	putative isomerase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (166/154); minor base A (1/0); total (167/154)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

TGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCGTAAACT  >  NC_000913/3339332‑3339511

tgtgCATCTGTGAGTTAAAGTAGCGAAAGAAGCCTGGCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGc                                                                                         >  2:1762923/1‑93 (MQ=255)

 gtgCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGcc                                                                                        <  1:1730022/93‑1 (MQ=255)

 gtgCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGcc                                                                                        >  1:4885434/1‑93 (MQ=255)

 gtgCAGCTGGGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGcc                                                                                        <  1:3878456/93‑1 (MQ=255)

   gCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCt                                                                                      <  1:1316865/93‑1 (MQ=255)

   gCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCt                                                                                      >  1:6128430/1‑93 (MQ=255)

   gCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCt                                                                                      >  2:6706341/1‑93 (MQ=255)

      tCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGc                                                                                   <  2:1061168/93‑1 (MQ=255)

      tCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGc                                                                                   <  2:2801219/93‑1 (MQ=255)

        tgtgtgTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTg                                                                                 <  1:5562471/93‑1 (MQ=255)

        tgtgtgTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTg                                                                                 <  1:4548153/93‑1 (MQ=255)

         gtgtgtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGt                                                                                >  2:6883859/1‑93 (MQ=255)

          tgtgtTAAAGTAGCGAAAGCAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTc                                                                               <  1:7019173/93‑1 (MQ=255)

          tgtgtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTc                                                                               >  2:1976809/1‑93 (MQ=255)

          tgtgtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTc                                                                               >  2:3282071/1‑93 (MQ=255)

           gtgtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTcg                                                                              <  2:1846577/93‑1 (MQ=255)

           gtgtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTcg                                                                              >  2:7306434/1‑93 (MQ=255)

            tgtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTcgc                                                                             >  1:7232578/1‑93 (MQ=255)

             gtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTcgcg                                                                            >  2:5618662/1‑93 (MQ=255)

             gtTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTcgcg                                                                            >  2:424000/1‑93 (MQ=255)

              ttAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCgcgc                                                                           >  1:7311965/1‑93 (MQ=255)

                 aaGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGt                                                                        <  1:3386787/93‑1 (MQ=255)

                      gCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTaaaa                                                                   <  1:7202304/93‑1 (MQ=255)

                      gCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTaaaa                                                                   <  1:5048986/93‑1 (MQ=255)

                      gCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACCCCCTCGCTGTCGCGCTGTTaaaa                                                                   <  2:5982489/93‑1 (MQ=255)

                         aaaGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCa                                                                <  2:4449291/93‑1 (MQ=255)

                         aaaGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCa                                                                >  1:1579823/1‑93 (MQ=255)

                         aaaGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCa                                                                <  2:7075796/93‑1 (MQ=255)

                          aagaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAc                                                               >  2:121497/1‑93 (MQ=255)

                          aagaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAc                                                               >  1:6075901/1‑93 (MQ=255)

                            gaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgc                                                             >  1:4234101/1‑93 (MQ=255)

                            gaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgc                                                             <  2:403926/93‑1 (MQ=255)

                            gaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgc                                                             <  2:608830/93‑1 (MQ=255)

                            gaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgc                                                             >  2:1689108/1‑93 (MQ=255)

                            gaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgc                                                             <  2:1331851/93‑1 (MQ=255)

                            gaagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgc                                                             <  1:6890525/93‑1 (MQ=255)

                             aagCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCAcgcg                                                            <  2:5476183/93‑1 (MQ=255)

                              agCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCgg                                                           >  1:287356/1‑93 (MQ=255)

                              agCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCgg                                                           >  2:3898209/1‑93 (MQ=255)

                               gCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGc                                                          >  2:3638258/1‑93 (MQ=255)

                                ccTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCt                                                         <  2:2443894/93‑1 (MQ=255)

                                 cTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCtt                                                        >  1:115571/1‑93 (MQ=255)

                                 cTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCtt                                                        <  1:3807971/93‑1 (MQ=255)

                                  tGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCttt                                                       >  2:7902861/1‑93 (MQ=255)

                                   gTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTa                                                      >  2:7556753/1‑93 (MQ=255)

                                     ccGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTAct                                                    <  1:4107245/93‑1 (MQ=255)

                                      cGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTActg                                                   >  2:12582/1‑93 (MQ=255)

                                       gTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTActgc                                                  >  2:7114445/1‑93 (MQ=255)

                                       gTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTActgc                                                  >  1:4084185/1‑93 (MQ=255)

                                        ttAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTActgct                                                 >  2:4712355/1‑93 (MQ=255)

                                           gggCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGaa                                              <  1:6405792/93‑1 (MQ=255)

                                           gggCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGaa                                              <  2:2750157/93‑1 (MQ=255)

                                           ggccTGGCCCCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGaa                                              <  2:5963423/90‑1 (MQ=255)

                                              cTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAt                                           >  2:2075953/1‑93 (MQ=255)

                                              cTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAt                                           <  1:5454617/93‑1 (MQ=255)

                                               tGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAtt                                          >  1:3599028/1‑93 (MQ=255)

                                                ggCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAttt                                         >  2:5715768/1‑93 (MQ=255)

                                                ggCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAttt                                         <  2:5256045/93‑1 (MQ=255)

                                                ggCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAttt                                         <  2:2627599/93‑1 (MQ=255)

                                                ggCACAGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAttt                                         <  2:4628300/93‑1 (MQ=255)

                                                 gCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCTCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAtttt                                        <  1:2805219/93‑1 (MQ=255)

                                                 gCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGAtttt                                        >  1:3461238/1‑93 (MQ=255)

                                                  cACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTg                                       >  1:5927818/1‑93 (MQ=255)

                                                   aCCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTgc                                      <  2:2395732/93‑1 (MQ=255)

                                                     cGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTgcgc                                    >  1:1461521/1‑93 (MQ=255)

                                                      gACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGct                                   >  2:5372068/1‑93 (MQ=255)

                                                      gACCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGct                                   >  2:624886/1‑93 (MQ=255)

                                                       aCCAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGctc                                  >  1:2505519/1‑93 (MQ=255)

                                                        ccAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGctct                                 <  1:5677268/93‑1 (MQ=255)

                                                        ccAGCAGCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGctct                                 <  1:3292588/93‑1 (MQ=255)

                                                         cagcagCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCtctc                                <  1:6454305/93‑1 (MQ=255)

                                                         cagcagCACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCtctc                                >  1:1925251/1‑93 (MQ=255)

                                                          agcagcACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTca                               <  2:6987038/93‑1 (MQ=255)

                                                          agcagcACCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTca                               <  2:3811819/93‑1 (MQ=255)

                                                           gcagcaCCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTcac                              >  1:5137682/1‑93 (MQ=255)

                                                             agcaCCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCacac                            <  2:3477865/93‑1 (MQ=255)

                                                              gcaCCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACAcc                           >  2:2905260/1‑93 (MQ=255)

                                                              gcaCCACCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACAcc                           >  1:7082465/1‑93 (MQ=255)

                                                               caccacCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCCCAccc                          >  1:7445400/1‑93 (MQ=255)

                                                               caccacCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACAccc                          >  2:3076121/1‑93 (MQ=255)

                                                               caccacCGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACAccc                          >  2:650388/1‑93 (MQ=255)

                                                                accaccGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCa                         <  1:7650424/93‑1 (MQ=255)

                                                                accaccGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCa                         >  2:3235395/1‑93 (MQ=255)

                                                                  caccGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAgg                       <  2:3141516/93‑1 (MQ=255)

                                                                   accGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAggc                      >  1:4772936/1‑93 (MQ=255)

                                                                     cGCCACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAggcgg                    >  1:6232735/1‑93 (MQ=255)

                                                                        cACGCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGgcgc                 >  2:4869367/1‑93 (MQ=255)

                                                                         aCGCTGGGTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCa                <  1:6565311/93‑1 (MQ=255)

                                                                           gCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACt              >  2:7646561/1‑93 (MQ=255)

                                                                           gCTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACt              <  1:7830045/93‑1 (MQ=255)

                                                                            cTGGTTATGGGGGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTg             >  2:4145498/1‑93 (MQ=255)

                                                                            cTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTg             <  1:3030732/93‑1 (MQ=255)

                                                                            cTGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTg             <  1:2021032/93‑1 (MQ=255)

                                                                             tGGTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTgg            >  1:951245/1‑93 (MQ=255)

                                                                              ggTTATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTggg           <  2:7210692/93‑1 (MQ=255)

                                                                              ggTTATGGGCGACGCCCTCGCTGTCGCGCTGATAAAAGCACGCGGCTTTACTGCTGCAGATTTTGCGCTCTCACACCCAGGCGGCGCACTggg           >  2:4562084/1‑93 (MQ=255)

                                                                                ttATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTc         >  2:1641931/1‑93 (MQ=255)

                                                                                 tATGGGCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCg        <  1:2401126/93‑1 (MQ=255)

                                                                                     ggCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCGTaaa    <  1:2835339/93‑1 (MQ=255)

                                                                                      gCGACGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCGTAAAc   <  2:1108246/93‑1 (MQ=255)

TGTGCATCTGTGTGTTAAAGTAGCGAAAGAAGCCTGTCCGTTAGGGCTGGCACCGACCAGCAGCACCACCGCCACGCTGGTTATGGGCGATGCCCTCGCTGTCGCGCTGTTAAAAGCACGCGGCTTTACTGCTGAAGATTTTGCGCTCTCACACCCAGGCGGCGCACTGGGTCGTAAACT  >  NC_000913/3339332‑3339511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 20 ≤` `ATCG/ATCG` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`