BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,732,250	C→T	A301V (GCC→GTC)	xylH →	XylH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,732,250	0	C	T	99.6%	1002.9 / ‑5.2	261	A301V (GCC→GTC)	xylH	XylH
Reads supporting (aligned to +/- strand): ref base C (1/0); new base T (130/130); total (131/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

Only 100 of 261 total aligned reads displayed.

GGCGGCAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGC  >  NC_000913/3732162‑3732338

ggcgtcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATc                                                                                      <  2:5203999/93‑1 (MQ=255)

ggcggcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATc                                                                                      <  1:3826194/93‑1 (MQ=255)

ggcggcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATc                                                                                      <  2:1851154/93‑1 (MQ=255)

ggcggcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATc                                                                                      <  2:1708378/93‑1 (MQ=255)

ggcggcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATc                                                                                      <  1:2459869/93‑1 (MQ=255)

    gcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCg                                                                                  >  1:4732484/1‑93 (MQ=255)

    gcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCg                                                                                  >  2:2598733/1‑93 (MQ=255)

    gcAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCg                                                                                  >  2:2935443/1‑93 (MQ=255)

     cAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCgg                                                                                 >  2:4293497/1‑93 (MQ=255)

       aTCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGAt                                                                               >  2:2091080/1‑93 (MQ=255)

       aTCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGAt                                                                               >  2:5627003/1‑93 (MQ=255)

       aTCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGAt                                                                               >  1:3564796/1‑93 (MQ=255)

         cTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTa                                                                             >  2:3161178/1‑93 (MQ=255)

          tGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTaa                                                                            <  2:1576635/93‑1 (MQ=255)

          tGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTaa                                                                            >  2:6205828/1‑93 (MQ=255)

           ggAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAAt                                                                           >  1:1747967/1‑93 (MQ=255)

            gAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATc                                                                          >  2:1974628/1‑93 (MQ=255)

               gcagcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCtt                                                                       <  1:436524/93‑1 (MQ=255)

               gcagcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCtt                                                                       <  2:3952690/93‑1 (MQ=255)

               gcagcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCtt                                                                       >  2:5062151/1‑93 (MQ=255)

                cagcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTa                                                                      <  2:4829058/93‑1 (MQ=255)

                 agcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCTCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAg                                                                     <  1:5563151/93‑1 (MQ=255)

                 agcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAg                                                                     >  2:2797684/1‑93 (MQ=255)

                 agcaCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAg                                                                     <  2:4017201/93‑1 (MQ=255)

                    aCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTtc                                                                  <  2:852746/93‑1 (MQ=255)

                    aCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTtc                                                                  <  1:1121820/93‑1 (MQ=255)

                    aCGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTtc                                                                  <  1:3521562/93‑1 (MQ=255)

                       tctctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCg                                                               <  1:466925/93‑1 (MQ=255)

                       tctctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCg                                                               <  1:5297702/93‑1 (MQ=255)

                       tctctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCg                                                               >  2:4401793/1‑93 (MQ=255)

                        ctctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGa                                                              >  1:2319484/1‑93 (MQ=255)

                         tctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGAc                                                             <  1:865814/93‑1 (MQ=255)

                         tctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGAc                                                             <  2:5362866/93‑1 (MQ=255)

                         tctcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGAc                                                             >  2:4198752/1‑93 (MQ=255)

                           tcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACtt                                                           <  2:3587844/93‑1 (MQ=255)

                           tcCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACtt                                                           <  1:5273860/93‑1 (MQ=255)

                            ccGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTg                                                          <  1:2943470/93‑1 (MQ=255)

                            ccGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTg                                                          <  1:6121769/93‑1 (MQ=255)

                            ccGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTg                                                          <  1:1965847/93‑1 (MQ=255)

                             cGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACttcg                                                         <  2:6245401/93‑3 (MQ=255)

                             cGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTgg                                                         <  2:2218683/93‑1 (MQ=255)

                             cGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTgg                                                         <  2:634873/93‑1 (MQ=255)

                              gggATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGgc                                                        >  2:3801257/1‑93 (MQ=255)

                               ggATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCTGGATTAATCCTTAGTTCTCGACTTGgcg                                                       >  1:3725102/1‑93 (MQ=255)

                               ggATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGgcg                                                       >  1:1045445/1‑93 (MQ=255)

                                gATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGgcgc                                                      >  1:4508965/1‑93 (MQ=255)

                                 aTTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCt                                                     <  2:4984404/93‑1 (MQ=255)

                                   tAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTgg                                                   >  1:5925737/1‑93 (MQ=255)

                                   tAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTgg                                                   >  1:3685730/1‑93 (MQ=255)

                                    aaCGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGt                                                  <  1:647693/93‑1 (MQ=255)

                                    aaCGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGt                                                  <  1:5876217/93‑1 (MQ=255)

                                    aaCGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGt                                                  >  2:2864410/1‑93 (MQ=255)

                                     aCGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGtt                                                 <  1:4280902/93‑1 (MQ=255)

                                     aCGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGtt                                                 >  2:5288390/1‑93 (MQ=255)

                                     aCGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGtt                                                 >  1:3133579/1‑93 (MQ=255)

                                         tGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCAcc                                             >  1:6317721/1‑93 (MQ=255)

                                         tGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCAcc                                             >  2:344154/1‑93 (MQ=255)

                                         tGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCAcc                                             <  2:5561118/93‑1 (MQ=255)

                                          gAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCt                                            >  2:2039048/1‑93 (MQ=255)

                                           aaCGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCtt                                           <  2:2482048/93‑1 (MQ=255)

                                            aCGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTc                                          >  2:3271441/1‑93 (MQ=255)

                                            aCGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTc                                          <  2:2889036/93‑1 (MQ=255)

                                             cGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCt                                         <  1:1330897/93‑1 (MQ=255)

                                             cGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCt                                         <  1:4493107/93‑1 (MQ=255)

                                              gCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTg                                        >  1:702382/1‑93 (MQ=255)

                                               cACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGc                                       <  1:3220189/93‑1 (MQ=255)

                                                aCCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCg                                      >  2:2780863/1‑93 (MQ=255)

                                                 ccAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCgg                                     >  1:5535762/1‑93 (MQ=255)

                                                 ccAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCgg                                     <  1:730088/93‑1 (MQ=255)

                                                   aaaCTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGa                                   >  1:3312115/1‑93 (MQ=255)

                                                     aCTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGaaa                                 >  2:560050/1‑93 (MQ=255)

                                                      cTTGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAat                                <  1:6043136/93‑1 (MQ=255)

                                                       ttGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAata                               >  1:2235016/1‑93 (MQ=255)

                                                        tGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAatat                              <  2:3030962/93‑1 (MQ=255)

                                                        tGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAatat                              <  1:1637913/93‑1 (MQ=255)

                                                        tGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAatat                              <  1:2099260/93‑1 (MQ=255)

                                                        tGCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAatat                              <  1:2690365/93‑1 (MQ=255)

                                                         gCCGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATc                             >  1:1874516/1‑93 (MQ=255)

                                                          ccGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCg                            >  1:4871558/1‑93 (MQ=255)

                                                           cGTGTTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGc                           <  2:1025406/93‑1 (MQ=255)

                                                              gtTCGCGATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGa                        <  2:365282/93‑1 (MQ=255)

                                                                 cgcgATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACt                     <  2:3754324/93‑1 (MQ=255)

                                                                   cgATTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTgg                   >  2:1284559/1‑93 (MQ=255)

                                                                     aTTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGAc                 >  2:695555/1‑93 (MQ=255)

                                                                     aTTAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGAc                 <  1:1792115/93‑1 (MQ=255)

                                                                       tAACGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGc               <  1:1441402/93‑1 (MQ=255)

                                                                         aCGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCaa             >  2:1347062/1‑93 (MQ=255)

                                                                          cGGATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAAt            <  1:3895551/93‑1 (MQ=255)

                                                                           ggATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAAtt           >  2:5811341/1‑93 (MQ=255)

                                                                           ggATTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAAtt           >  2:1213970/1‑93 (MQ=255)

                                                                             aTTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTgc         >  2:542553/1‑93 (MQ=255)

                                                                             aTTAATGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTgc         <  2:5986708/93‑1 (MQ=255)

                                                                                aaTGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTgcagc      >  2:4598363/1‑93 (MQ=255)

                                                                                aaTGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTgcagc      >  2:1325033/1‑93 (MQ=255)

                                                                                  tGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGGAt    >  1:1564686/1‑93 (MQ=255)

                                                                                  tGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCAt    >  1:597709/1‑93 (MQ=255)

                                                                                  tGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCAt    <  1:2768720/93‑1 (MQ=255)

                                                                                  tGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCAt    <  1:6198019/93‑1 (MQ=255)

                                                                                  tGGTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCAt    <  1:58274/93‑1 (MQ=255)

                                                                                   ggTAGTCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATg   >  2:729701/1‑93 (MQ=255)

GGCGGCAATCTGGAAGCAGCACGTCTCTCCGGGATTAACGTTGAACGCACCAAACTTGCCGTGTTCGCGATTAACGGATTAATGGTAGCCATCGCCGGATTAATCCTTAGTTCTCGACTTGGCGCTGGTTCACCTTCTGCGGGAAATATCGCCGAACTGGACGCAATTGCAGCATGC  >  NC_000913/3732162‑3732338

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