BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,782,932	T→C	Y16H (TAT→CAT)	yibO →	putative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,782,932	0	T	C	100.0%	923.8 / NA	245	Y16H (TAT→CAT)	yibO	putative 2,3‑bisphosphoglycerate‑independent phosphoglycerate mutase
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (129/116); total (129/116)

Only 100 of 276 total aligned reads displayed.

AGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACC  >  NC_000913/3782842‑3783021

agtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccat                                                                                         >  2:1969480/1‑93 (MQ=255)

agtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccat                                                                                         <  2:342172/93‑1 (MQ=255)

agtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccat                                                                                         >  1:3414213/1‑93 (MQ=255)

 gtaAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccatg                                                                                        <  2:896808/93‑1 (MQ=255)

    aaaTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGatggccatggct                                                                                     >  1:6284294/1‑93 (MQ=255)

     aaTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGAtggccatggcta                                                                                    >  2:3276165/1‑93 (MQ=255)

      aTTACGCAAATTTTGGCTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATggccatggctat                                                                                   <  1:5529963/93‑1 (MQ=255)

      aTTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATggccatggctat                                                                                   <  2:5357300/93‑1 (MQ=255)

       ttACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATc                                                                                  <  2:5940371/93‑1 (MQ=255)

          cGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATcgcg                                                                               <  1:428431/93‑1 (MQ=255)

          cGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATcgcg                                                                               <  2:1470327/93‑1 (MQ=255)

               aTTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCgaagaa                                                                          <  1:6509641/93‑1 (MQ=255)

                ttttGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAc                                                                         <  1:1376482/93‑1 (MQ=255)

                 tttGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAca                                                                        <  2:2544486/93‑1 (MQ=255)

                     aCTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAAcagcag                                                                    <  1:171172/93‑1 (MQ=255)

                      ctctTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAgg                                                                   >  2:6461130/1‑93 (MQ=255)

                       tctTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACTTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGa                                                                  <  1:5339959/93‑1 (MQ=255)

                       tctTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGa                                                                  <  1:3777819/93‑1 (MQ=255)

                        ctTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGAt                                                                 <  2:1020454/93‑1 (MQ=255)

                        ctTGAGTATGAGGTTGTCGCAATGATGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGAt                                                                 <  2:1305356/93‑1 (MQ=255)

                         ttGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATa                                                                >  2:4770315/1‑93 (MQ=255)

                         ttGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATa                                                                <  2:4589112/93‑1 (MQ=255)

                         ttGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATa                                                                <  2:76754/93‑1 (MQ=255)

                            aGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACg                                                             >  2:5164163/1‑93 (MQ=255)

                             gTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAGCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGc                                                            <  1:5683489/93‑1 (MQ=255)

                             gTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGc                                                            <  1:3195720/93‑1 (MQ=255)

                              tATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGcc                                                           <  2:3878123/93‑1 (MQ=255)

                              tATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGcc                                                           <  2:5638810/93‑1 (MQ=255)

                              tATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGcc                                                           >  2:1380880/1‑93 (MQ=255)

                              tATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTCATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGcc                                                           <  1:396107/93‑1 (MQ=255)

                               aTGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCa                                                          >  1:5209706/1‑93 (MQ=255)

                                tGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAt                                                         >  1:977325/1‑93 (MQ=255)

                                 gAGGTTGTCGCAATGTTGGTTTCTAAAAATCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAtt                                                        <  1:3449924/93‑1 (MQ=255)

                                  aGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAttt                                                       >  1:739507/1‑93 (MQ=255)

                                   ggTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAtttt                                                      >  2:4109485/1‑93 (MQ=255)

                                   ggTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAtttt                                                      >  1:1919523/1‑93 (MQ=255)

                                   ggTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCAtttt                                                      >  2:4200839/1‑93 (MQ=255)

                                     ttGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    >  2:4296408/1‑93 (MQ=255)

                                     ttGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    <  2:4315865/93‑1 (MQ=255)

                                     ttGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    >  2:4862193/1‑93 (MQ=255)

                                     ttGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    <  2:4967040/93‑1 (MQ=255)

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

                                     ttGTCGCAATGATGGGTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCCTTTCTa                                                    >  2:6215282/1‑93 (MQ=255)

                                     tcgTCGCAATGTTGGTTTCTAAAAACCCTATGGTACTGGTGATTCTGAATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTa                                                    <  1:2795061/91‑1 (MQ=255)

                                      tGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAg                                                   <  1:2924223/93‑1 (MQ=255)

                                        tCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTg                                                 >  1:4701913/1‑93 (MQ=255)

                                          gCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCt                                               >  1:5021085/1‑93 (MQ=255)

                                           cAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTa                                              >  1:3670507/1‑93 (MQ=255)

                                           cAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTa                                              <  2:2765491/93‑1 (MQ=255)

                                            aaTGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTaa                                             >  2:3929806/1‑93 (MQ=255)

                                             aTGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTaaa                                            >  1:2290234/1‑93 (MQ=255)

                                              tGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTaaaa                                           >  2:1890012/1‑93 (MQ=255)

                                               gTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAc                                          >  1:5775686/1‑93 (MQ=255)

                                                ttGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAcc                                         >  1:4221522/1‑93 (MQ=255)

                                                 tGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAccc                                        <  2:5996642/93‑1 (MQ=255)

                                                  ggTTTCTAAAAATCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAcccc                                       <  2:1104193/93‑1 (MQ=255)

                                                  ggTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAAcccc                                       <  2:135719/93‑1 (MQ=255)

                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      <  2:1717611/93‑1 (MQ=255)

                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      <  1:2257745/93‑1 (MQ=255)

                                                   gTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCg                                      >  1:5965733/1‑93 (MQ=255)

                                                    tttCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCgg                                     >  2:23045/1‑93 (MQ=255)

                                                     ttCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGt                                    >  2:6516349/1‑93 (MQ=255)

                                                     ttCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGt                                    <  1:6376458/93‑1 (MQ=255)

                                                      tCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTa                                   <  1:5125710/93‑1 (MQ=255)

                                                      tCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTa                                   <  2:366340/93‑1 (MQ=255)

                                                        tAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAAt                                 <  2:1601862/93‑1 (MQ=255)

                                                        tAAAAAACCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAAt                                 >  2:5022132/1‑93 (MQ=255)

                                                         aaaaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATg                                <  2:4821286/93‑1 (MQ=255)

                                                          aaaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATgg                               >  2:2337536/1‑93 (MQ=255)

                                                          aaaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATgg                               >  2:6610479/1‑93 (MQ=255)

                                                          aaaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATgg                               >  1:513555/1‑93 (MQ=255)

                                                           aaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGa                              <  1:1796716/93‑1 (MQ=255)

                                                           aaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGa                              <  2:3583777/93‑1 (MQ=255)

                                                           aaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGa                              >  2:5892529/1‑93 (MQ=255)

                                                           aaaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGa                              >  2:3157188/1‑93 (MQ=255)

                                                            aaaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGAt                             <  2:4064901/93‑1 (MQ=255)

                                                             aaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATg                            <  1:6047622/93‑1 (MQ=255)

                                                             aaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATg                            >  1:2611553/1‑93 (MQ=255)

                                                             aaCCTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATg                            >  1:1136924/1‑93 (MQ=255)

                                                                cTATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCAc                         >  2:3019122/1‑93 (MQ=255)

                                                                 tATGGTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACt                        >  1:2953845/1‑93 (MQ=255)

                                                                    ggTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACtgtg                     <  1:4321501/93‑1 (MQ=255)

                                                                    ggTACTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACtgtg                     <  1:964567/93‑1 (MQ=255)

                                                                        cTGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGcc                 >  2:5324655/1‑93 (MQ=255)

                                                                         tGGTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCa                >  1:3500805/1‑93 (MQ=255)

                                                                         tGGTCATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTCAGCCa                >  2:3042797/1‑93 (MQ=255)

                                                                          ggTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACACCGGTAATGGATGCACTGTGGGACaa               >  2:3373962/1‑93 (MQ=255)

                                                                          ggTGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGAGAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCaa               >  2:4673519/1‑93 (MQ=255)

                                                                            tGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATc             >  1:2001332/1‑93 (MQ=255)

                                                                            tGATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATc             <  2:6425473/93‑1 (MQ=255)

                                                                             gATTCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCg            <  1:4572308/93‑1 (MQ=255)

                                                                               ttCTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTc          <  1:2532771/93‑1 (MQ=255)

                                                                                 cTGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCg        <  1:5171082/93‑1 (MQ=255)

                                                                                  tGGATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGc       >  1:4514405/1‑93 (MQ=255)

                                                                                   ggATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCa      >  2:3056082/1‑93 (MQ=255)

                                                                                   ggATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCa      >  1:729157/1‑93 (MQ=255)

                                                                                   ggATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCa      <  1:62573/93‑1 (MQ=255)

                                                                                    gATGGCCATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCAt     >  2:5265297/1‑93 (MQ=255)

                                                                                     atggccatggctATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATa    <  1:1201249/93‑1 (MQ=255)

                                                                                     atggccatggctATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATa    >  2:4423448/1‑93 (MQ=255)

AGTAAAATTACGCAAATTTTGACTCTTGAGTATGAGGTTGTCGCAATGTTGGTTTCTAAAAAACCTATGGTACTGGTGATTCTGGATGGCTATGGCTATCGCGAAGAACAGCAGGATAACGCCATTTTTAGTGCTAAAACCCCGGTAATGGATGCACTGTGGGCCAATCGTCCGCATACC  >  NC_000913/3782842‑3783021

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