BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,726,481	T→C	T174A (ACC→GCC)	xylB ←	xylulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,726,481	0	T	C	85.8%	826.2 / 118.4	352	T174A (ACC→GCC)	xylB	xylulokinase
Reads supporting (aligned to +/- strand): ref base T (22/28); new base C (146/156); total (168/184)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.85e-01

Only 100 of 352 total aligned reads displayed.

CGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGAC  >  NC_000913/3726391‑3726570

 gTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGcc                                                                                        >  2:5492664/1‑93 (MQ=255)

  tataATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCa                                                                                       >  1:5797485/1‑93 (MQ=255)

   ataATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAg                                                                                      >  1:4823350/1‑93 (MQ=255)

   ataATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAg                                                                                      <  2:5330312/93‑1 (MQ=255)

     aaTGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCt                                                                                    <  1:90584/93‑1 (MQ=255)

      aTGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTg                                                                                   >  2:6349576/1‑93 (MQ=255)

       tGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGc                                                                                  >  2:1369120/1‑93 (MQ=255)

       tGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGc                                                                                  >  1:5382500/1‑93 (MQ=255)

       tGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGc                                                                                  <  1:4307653/93‑1 (MQ=255)

         cGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGt                                                                                >  2:6120870/1‑93 (MQ=255)

          gggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTc                                                                               >  2:1602666/1‑93 (MQ=255)

          gggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTc                                                                               >  1:2015613/1‑93 (MQ=255)

           ggCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCa                                                                              >  2:3775749/1‑93 (MQ=255)

             cATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGa                                                                            <  2:2387938/93‑1 (MQ=255)

               tCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACa                                                                          >  2:3201951/1‑93 (MQ=255)

               tCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACa                                                                          >  1:2397714/1‑93 (MQ=255)

                 tGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACata                                                                        >  2:4125586/1‑93 (MQ=255)

                 tGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACata                                                                        <  2:255360/93‑1 (MQ=255)

                 tGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACata                                                                        <  2:4583430/93‑1 (MQ=255)

                   gTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATc                                                                      <  2:6614733/93‑1 (MQ=255)

                     cACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGc                                                                    <  2:3660970/93‑1 (MQ=255)

                     cACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGc                                                                    <  2:7178987/93‑1 (MQ=255)

                      aCGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCt                                                                   >  1:3373585/1‑93 (MQ=255)

                      aCGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCt                                                                   >  1:2277325/1‑93 (MQ=255)

                       cGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTg                                                                  <  1:4458969/93‑1 (MQ=255)

                       cGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTg                                                                  >  1:1321244/1‑93 (MQ=255)

                        gagaTAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTgg                                                                 <  1:4347744/93‑1 (MQ=255)

                          gaaaaGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCa                                                               <  1:2708164/90‑1 (MQ=255)

                          gaTAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCa                                                               >  2:4275017/1‑93 (MQ=255)

                          gaTAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCa                                                               <  2:4009651/93‑1 (MQ=255)

                            tAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCaaa                                                             <  2:6546728/93‑1 (MQ=255)

                            tAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCaaa                                                             >  1:5642649/1‑93 (MQ=255)

                            tAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCAAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCaaa                                                             <  1:1712134/93‑1 (MQ=255)

                             aaGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAAc                                                            <  2:132040/93‑1 (MQ=255)

                             aaGCCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAAc                                                            <  2:4915637/93‑1 (MQ=255)

                               gTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTc                                                          <  1:6418359/93‑1 (MQ=255)

                                tCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTcc                                                         <  2:1931600/93‑1 (MQ=255)

                                tCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTcc                                                         >  1:1042812/1‑93 (MQ=255)

                                  gCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTcccc                                                       <  1:4170098/93‑1 (MQ=255)

                                  gCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTcccc                                                       <  1:3620467/93‑1 (MQ=255)

                                  gCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTcccc                                                       <  1:3510211/93‑1 (MQ=255)

                                    aaGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCg                                                     >  2:5707210/1‑93 (MQ=255)

                                     aGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGt                                                    >  2:3613195/1‑93 (MQ=255)

                                       ccTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCa                                                  >  1:668462/1‑93 (MQ=255)

                                       ccTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCa                                                  <  1:1369627/93‑1 (MQ=255)

                                         tGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATa                                                >  1:5564403/1‑93 (MQ=255)

                                            agcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGc                                             >  1:1523555/1‑93 (MQ=255)

                                             gcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCa                                            <  1:6927374/93‑1 (MQ=255)

                                             gcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCa                                            <  2:813009/93‑1 (MQ=255)

                                             gcaTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCa                                            <  1:7104208/93‑1 (MQ=255)

                                              gaTGGCGTCCCTCCAGTCCCGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAg                                           <  2:4151986/92‑1 (MQ=255)

                                               aTGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGa                                          <  2:5075536/93‑1 (MQ=255)

                                                tGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAgnn                                         >  2:2695047/1‑91 (MQ=255)

                                                tGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGAc                                         >  2:4451371/1‑93 (MQ=255)

                                                 gACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACg                                        >  2:3336247/1‑93 (MQ=255)

                                                    gTCCCTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCaa                                     <  1:28865/93‑1 (MQ=255)

                                                    gTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCaa                                     <  2:47276/93‑1 (MQ=255)

                                                    gTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCaa                                     <  1:4669881/93‑1 (MQ=255)

                                                    gTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCaa                                     <  1:1599985/93‑1 (MQ=255)

                                                       aCTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATCCGCAGACGCAAGTa                                  <  2:5185529/93‑1 (MQ=255)

                                                       aCTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTa                                  >  2:2974374/1‑93 (MQ=255)

                                                       aCTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTa                                  >  2:3568959/1‑93 (MQ=255)

                                                        cTCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTaa                                 <  2:5238380/93‑1 (MQ=255)

                                                         tCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAAt                                >  1:3139815/1‑93 (MQ=255)

                                                         tCCAGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAAt                                <  2:390977/93‑1 (MQ=255)

                                                            aGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCtt                             <  1:413073/93‑1 (MQ=255)

                                                            aGTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCtt                             <  2:766031/93‑1 (MQ=255)

                                                             gTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCttt                            <  2:5830695/93‑1 (MQ=255)

                                                             gTCACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCttt                            >  1:3046159/1‑93 (MQ=255)

                                                               cACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCg                          <  1:6757781/93‑1 (MQ=255)

                                                               cACGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCg                          >  2:6182060/1‑93 (MQ=255)

                                                                aCGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCgg                         >  1:2521914/1‑93 (MQ=255)

                                                                 cGCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGt                        <  2:2816651/93‑1 (MQ=255)

                                                                  gCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGta                       <  2:7309549/93‑1 (MQ=255)

                                                                  gCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGta                       <  1:3232920/93‑1 (MQ=255)

                                                                  gCTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGta                       <  1:7330086/93‑1 (MQ=255)

                                                                   cTTTGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGtaa                      >  2:1815860/1‑93 (MQ=255)

                                                                    tttGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGtaat                     >  1:6889886/1‑93 (MQ=255)

                                                                     ttGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGtaata                    <  2:3435601/93‑1 (MQ=255)

                                                                      tGCGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGtaataa                   <  2:7071305/93‑1 (MQ=255)

                                                                        cGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAataata                 <  2:5102242/93‑1 (MQ=255)

                                                                        cGACATCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAataata                 <  1:597218/93‑1 (MQ=255)

                                                                         gACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATAc                >  2:5575618/1‑93 (MQ=255)

                                                                            aTCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACttt             <  1:3855860/93‑1 (MQ=255)

                                                                             tCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTg            >  1:531321/1‑93 (MQ=255)

                                                                             tCCAGCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTg            >  2:6949383/1‑93 (MQ=255)

                                                                                 gCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGa        >  1:1096230/1‑93 (MQ=255)

                                                                                 gCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGa        <  2:1251215/93‑1 (MQ=255)

                                                                                 gCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGa        >  1:6365607/1‑93 (MQ=255)

                                                                                 gCCACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGa        >  2:3509916/1‑93 (MQ=255)

                                                                                  ccACATGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGAt       <  2:2140258/93‑1 (MQ=255)

                                                                                   cacaTGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGAtt      >  2:6846605/1‑93 (MQ=255)

                                                                                   cacaTGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGAtt      <  1:3990186/93‑1 (MQ=255)

                                                                                    acaTGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGAttt     >  2:4955669/1‑93 (MQ=255)

                                                                                    acaTGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGAttt     >  1:283569/1‑93 (MQ=255)

                                                                                     caTGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTg    >  1:974230/1‑93 (MQ=255)

                                                                                     caTGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTg    <  1:6049510/93‑1 (MQ=255)

                                                                                     caTGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTg    <  1:6155654/93‑1 (MQ=255)

                                                                                       tGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGAc  <  2:6003115/93‑1 (MQ=255)

                                                                                       tGGCGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGAc  >  2:622968/1‑93 (MQ=255)

CGTATAATGCGGGCATCTGGTCACGAGATAAGTCGCAAGCCTGCAGCATGACGTCACTCCAGTCACGCTTTGCGACATCCAGCCACATGGTGCCAGCTGCGTCAGACATATCGCTGGCAAACTCCCCCGTCATACGCAGACGCAAGTAATCTTTCGGTAATAATACTTTGTCGATTTGAC  >  NC_000913/3726391‑3726570

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