BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	488,917	A→G	E1053G (GAG→GGG)	aefA →	putative alpha helix protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	488,917	0	A	G	100.0%	722.5 / NA	221	E1053G (GAG→GGG)	aefA	putative alpha helix protein
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (117/104); total (117/104)

Only 100 of 221 total aligned reads displayed.

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006

   agtcggcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGc                                                                                      <  2:291108/90‑1 (MQ=255)

        gcgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTg                                                                                 <  1:3180541/93‑1 (MQ=255)

         cgACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGt                                                                                <  1:5343303/93‑1 (MQ=255)

           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCCCGTTGGATCATGGGCTGCGTCTGTAt                                                                              <  2:5863262/93‑1 (MQ=255)

           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGCTCATGGGCTGCGTCTGtct                                                                              >  2:822869/1‑91 (MQ=255)

           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              >  1:5692829/1‑93 (MQ=255)

           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              >  2:352122/1‑93 (MQ=255)

           aCTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAt                                                                              <  2:6037014/93‑1 (MQ=255)

             tGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgt                                                                            <  1:5544191/93‑1 (MQ=255)

              gAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTAtgtg                                                                           >  1:3591403/1‑93 (MQ=255)

               aGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGc                                                                          >  1:1294342/1‑93 (MQ=255)

                 cACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGt                                                                        <  2:3238662/93‑1 (MQ=255)

                 cACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGt                                                                        <  1:2394811/93‑1 (MQ=255)

                   cccAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGa                                                                      >  1:568104/1‑93 (MQ=255)

                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  1:3347842/93‑1 (MQ=255)

                     cAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAAc                                                                    <  1:3390640/93‑1 (MQ=255)

                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   >  1:2472721/1‑93 (MQ=255)

                      aaGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACt                                                                   >  2:5483646/1‑93 (MQ=255)

                       aGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTg                                                                  >  2:2431759/1‑93 (MQ=255)

                        gggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGc                                                                 >  1:4797443/1‑93 (MQ=255)

                        gggTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGc                                                                 <  1:188466/93‑1 (MQ=255)

                          gTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGt                                                               <  1:1829648/93‑1 (MQ=255)

                            gATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTgc                                                             >  1:1356879/1‑92 (MQ=255)

                             atccACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGAc                                                            <  1:2518914/90‑1 (MQ=255)

                                cACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGt                                                         <  1:5531526/93‑1 (MQ=255)

                                 aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        >  2:5700028/1‑93 (MQ=255)

                                 aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        <  1:4296381/93‑1 (MQ=255)

                                 aCGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTa                                                        <  2:1822364/93‑1 (MQ=255)

                                  cGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAg                                                       >  2:3923041/1‑93 (MQ=255)

                                   gacccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGt                                                      <  2:28537/90‑1 (MQ=255)

                                    aaCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTc                                                     <  2:1613334/93‑1 (MQ=255)

                                    aaCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTc                                                     >  1:5027072/1‑93 (MQ=255)

                                    aaCCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTc                                                     >  1:5500435/1‑93 (MQ=255)

                                      ccAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGt                                                   >  1:4183430/1‑93 (MQ=255)

                                       cAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTa                                                  >  2:712402/1‑93 (MQ=255)

                                        aaTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 >  1:6297685/1‑93 (MQ=255)

                                        aaTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTAc                                                 <  2:2242800/93‑1 (MQ=255)

                                         aTGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACt                                                <  1:5125566/93‑1 (MQ=255)

                                          tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTg                                               <  2:294141/93‑1 (MQ=255)

                                          tGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTg                                               <  1:2975176/93‑1 (MQ=255)

                                           gCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGt                                              <  2:1544720/93‑1 (MQ=255)

                                             cGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            <  2:2692712/93‑1 (MQ=255)

                                             cGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCg                                            >  1:3912478/1‑93 (MQ=255)

                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           >  1:3417786/1‑93 (MQ=255)

                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           <  2:1711016/93‑1 (MQ=255)

                                              ggAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGa                                           >  2:5008290/1‑93 (MQ=255)

                                               gAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGAt                                          >  1:1261925/1‑93 (MQ=255)

                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         >  2:3770133/1‑93 (MQ=255)

                                                aaGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATg                                         <  2:3558079/93‑1 (MQ=255)

                                                 aGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGa                                        <  1:880698/93‑1 (MQ=255)

                                                  gTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCTTGACCGTAGTCGTACTGTCGATGAg                                       >  2:2821274/1‑93 (MQ=255)

                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                     <  1:4272898/93‑1 (MQ=255)

                                                    cttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCt                                     >  2:2993183/1‑93 (MQ=255)

                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  2:4052918/1‑93 (MQ=255)

                                                     ttcttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTg                                    >  1:2621658/1‑93 (MQ=255)

                                                       cttTACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGaa                                  <  2:1850990/93‑1 (MQ=255)

                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGTGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  2:5254725/1‑93 (MQ=255)

                                                        tttACGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAAc                                 >  2:5125419/1‑93 (MQ=255)

                                                           aCGGCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATTAGCTGAACCGt                              >  2:5314767/1‑93 (MQ=255)

                                                             ggCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTAc                            <  1:4586350/93‑1 (MQ=255)

                                                              gCATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACt                           >  2:5731774/1‑93 (MQ=255)

                                                               cATTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTa                          >  1:914428/1‑93 (MQ=255)

                                                                aTTTGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTAt                         <  1:6092773/93‑1 (MQ=255)

                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        >  1:4098604/1‑93 (MQ=255)

                                                                 tttGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATc                        <  2:4167821/93‑1 (MQ=255)

                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  1:1677198/93‑1 (MQ=255)

                                                                   tGGTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGa                      <  2:4238990/93‑1 (MQ=255)

                                                                    ggTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     <  2:893706/93‑1 (MQ=255)

                                                                    ggTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGAt                     <  1:1610922/93‑1 (MQ=255)

                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  1:2726374/93‑1 (MQ=255)

                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  2:1199151/93‑1 (MQ=255)

                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  2:4802962/93‑1 (MQ=255)

                                                                     gTGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATc                    <  1:5912982/93‑1 (MQ=255)

                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   >  2:6214641/1‑93 (MQ=255)

                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   <  2:600170/93‑1 (MQ=255)

                                                                      tGCCAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCa                   >  2:527552/1‑93 (MQ=255)

                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 >  1:3524709/1‑93 (MQ=255)

                                                                        ccAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGc                 <  2:4369145/93‑1 (MQ=255)

                                                                         cAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCt                >  2:316077/1‑93 (MQ=255)

                                                                         cAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCt                >  1:3069617/1‑93 (MQ=255)

                                                                         cAGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCt                >  1:958876/1‑93 (MQ=255)

                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               <  2:3377404/93‑1 (MQ=255)

                                                                          aGCACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtg               >  1:5017635/1‑93 (MQ=255)

                                                                            cACGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCtgtg             >  1:3608963/1‑93 (MQ=255)

                                                                             aCGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGc            <  2:2599359/93‑1 (MQ=255)

                                                                              cGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGcc           <  2:552677/93‑1 (MQ=255)

                                                                              cGTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGcc           <  1:4769183/93‑1 (MQ=255)

                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          <  1:4345816/93‑1 (MQ=255)

                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  1:4005924/1‑93 (MQ=255)

                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          >  1:698121/1‑93 (MQ=255)

                                                                               gTTGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCg          <  2:2421368/93‑1 (MQ=255)

                                                                                ttGGATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGt         >  1:4827580/1‑93 (MQ=255)

                                                                                  ggATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGa       >  2:4505549/1‑93 (MQ=255)

                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      <  2:4316339/93‑1 (MQ=255)

                                                                                   gATCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaa      >  1:2075794/1‑93 (MQ=255)

                                                                                    aTCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa     <  2:5149763/93‑1 (MQ=255)

                                                                                    aTCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaa     >  2:2260187/1‑93 (MQ=255)

                                                                                     tCATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGaaaa    <  1:2796475/93‑1 (MQ=255)

                                                                                      cATGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAAc   >  1:730525/1‑93 (MQ=255)

                                                                                       aTGGGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACg  >  2:2057734/1‑93 (MQ=255)

TGAAGGCGGCGACTGAGCACCCAAGGGTGATGCACGAACCAATGCCGGAAGTCTTCTTTACGGCATTTGGTGCCAGCACGTTGGATCATGAGCTGCGTCTGTATGTGCGTGAACTGCGTGACCGTAGTCGTACTGTCGATGAGCTGAACCGTACTATCGATCAGCTGTGCCGTGAAAACG  >  NC_000913/488827‑489006

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