BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,973,216	C→T	A153A (GCC→GCT)	wecE →	putative regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,973,216	0	C	T	100.0%	878.5 / NA	234	A153A (GCC→GCT)	wecE	putative regulator
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (123/111); total (123/111)

Only 100 of 247 total aligned reads displayed.

CGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973126‑3973305

cGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGAtgctgc                                                                                         <  2:2012210/93‑1 (MQ=255)

   gCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCa                                                                                      <  1:911472/93‑1 (MQ=255)

     cGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAgg                                                                                    <  1:437561/93‑1 (MQ=255)

     cGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAgg                                                                                    <  2:1350889/93‑1 (MQ=255)

     cGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAgg                                                                                    <  2:5959535/93‑1 (MQ=255)

      ggTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAggg                                                                                   >  2:3822009/1‑93 (MQ=255)

      ggTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAggg                                                                                   >  1:4081435/1‑93 (MQ=255)

         ccATTACGCGGGTTTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGt                                                                                <  1:177257/93‑1 (MQ=255)

         ccATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGt                                                                                <  2:2036273/93‑1 (MQ=255)

         ccATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGt                                                                                <  2:1325310/93‑1 (MQ=255)

          cATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTg                                                                               <  2:4622306/93‑1 (MQ=255)

           aTTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGCAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCCGGGCGGGa                                                                              >  1:863039/1‑93 (MQ=255)

           aTTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGTTAGAAGATGCTGCTCAGGGCGTGa                                                                              <  1:2359564/93‑1 (MQ=255)

           aTTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCTAGGCGTGa                                                                              >  1:6030257/1‑93 (MQ=255)

           aTTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGa                                                                              >  1:2448721/1‑93 (MQ=255)

            ttACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGAt                                                                             <  2:581867/93‑1 (MQ=255)

               cgcgGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTc                                                                          <  1:1859037/93‑1 (MQ=255)

                  gggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCAc                                                                       >  1:1710002/1‑93 (MQ=255)

                   ggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACt                                                                      <  1:1920481/93‑1 (MQ=255)

                     tgtgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTa                                                                    >  1:3741438/1‑93 (MQ=255)

                      gtgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTAc                                                                   <  2:564896/93‑1 (MQ=255)

                          ccTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAg                                                               <  2:358391/93‑1 (MQ=255)

                           cTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAgg                                                              >  1:750829/1‑93 (MQ=255)

                           cTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAgg                                                              <  2:1707080/93‑1 (MQ=255)

                            tGCGAAATGGACAGCATTATGGCGTTGTCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             <  2:1488536/93‑1 (MQ=255)

                            tGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             >  1:5001700/1‑93 (MQ=255)

                               gAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGt                                                          <  1:453321/93‑1 (MQ=255)

                               gAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGt                                                          >  2:2404848/1‑93 (MQ=255)

                                 aaTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGc                                                        <  1:5776455/93‑1 (MQ=255)

                                 aaTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGc                                                        >  2:191555/1‑93 (MQ=255)

                                   tGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCAc                                                      <  2:2166037/93‑1 (MQ=255)

                                    ggACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACt                                                     >  1:1401883/1‑93 (MQ=255)

                                    ggACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACt                                                     >  2:5035090/1‑93 (MQ=255)

                                    ggACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACt                                                     >  1:3476137/1‑93 (MQ=255)

                                     gACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTg                                                    <  1:5332270/93‑1 (MQ=255)

                                       cacCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTggg                                                  >  1:3522659/1‑93 (MQ=255)

                                       cacCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTggg                                                  <  1:2973381/93‑1 (MQ=255)

                                         ccATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGaa                                                <  1:4900389/93‑1 (MQ=255)

                                         ccATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGaa                                                >  1:567864/1‑93 (MQ=255)

                                         ccATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGaa                                                <  2:1493496/93‑1 (MQ=255)

                                          cATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACACAGGGCGTGCACTGGGAAc                                               >  1:781551/1‑93 (MQ=255)

                                            ttATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                             >  1:5180635/1‑93 (MQ=255)

                                            ttATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                             >  1:3263232/1‑93 (MQ=255)

                                            ttATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                             >  2:1650218/1‑93 (MQ=255)

                                            ttATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                             >  2:4490406/1‑93 (MQ=255)

                                            ttATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCa                                             <  2:5293738/93‑1 (MQ=255)

                                               tGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCAGTTACAAAGGGCGTGCACTGGGAACCATTg                                          >  2:6191944/1‑93 (MQ=255)

                                               tGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTg                                          >  2:66591/1‑93 (MQ=255)

                                               tGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTg                                          <  2:1001684/93‑1 (MQ=255)

                                                ggCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTgg                                         >  2:3318868/1‑93 (MQ=255)

                                                ggCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTgg                                         >  2:8220/1‑93 (MQ=255)

                                                ggCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTgg                                         >  1:1318174/1‑93 (MQ=255)

                                                 gCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGt                                        <  1:1006213/93‑1 (MQ=255)

                                                  cGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTc                                       >  2:289333/1‑93 (MQ=255)

                                                   gTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCa                                      >  1:2595586/1‑93 (MQ=255)

                                                    ttGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCat                                     >  2:487977/1‑93 (MQ=255)

                                                    ttGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCat                                     <  1:2828392/93‑1 (MQ=255)

                                                    ttGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCat                                     >  2:4345329/1‑93 (MQ=255)

                                                         gAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGTGCGTGCACTGGGAACCATTGGTCATATTgg                                >  1:6257385/1‑93 (MQ=255)

                                                         gAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTgg                                <  2:1109532/93‑1 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               >  1:683063/1‑93 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               >  2:2307693/1‑93 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               >  1:3976555/1‑93 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               >  1:1246038/1‑93 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               >  1:2102991/1‑93 (MQ=255)

                                                            aaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctg                             >  1:478425/1‑93 (MQ=255)

                                                            aaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctg                             <  2:1762652/93‑1 (MQ=255)

                                                               gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCtt                          >  1:1127295/1‑93 (MQ=255)

                                                               gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCtt                          >  1:5887302/1‑93 (MQ=255)

                                                               gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCtt                          >  2:61132/1‑93 (MQ=255)

                                                                cATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCttt                         <  2:2185811/93‑1 (MQ=255)

                                                                cATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCttt                         <  2:5998248/93‑1 (MQ=255)

                                                                cATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCttt                         <  1:1298895/93‑1 (MQ=255)

                                                                 aTAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTa                        >  1:3626895/1‑93 (MQ=255)

                                                                 aTAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTa                        <  1:5983984/93‑1 (MQ=255)

                                                                  tAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAg                       <  2:3432090/93‑1 (MQ=255)

                                                                   aaTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGc                      <  2:1188163/93‑1 (MQ=255)

                                                                   aaTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGc                      <  1:2308434/93‑1 (MQ=255)

                                                                    aTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCt                     >  1:4790083/1‑93 (MQ=255)

                                                                    aTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCt                     >  1:6014769/1‑93 (MQ=255)

                                                                       tgtttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTcc                  <  1:3159895/93‑1 (MQ=255)

                                                                       tgtttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTcc                  >  2:4905420/1‑93 (MQ=255)

                                                                        gtttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCa                 >  1:3052720/1‑93 (MQ=255)

                                                                        gtttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCa                 <  1:59270/93‑1 (MQ=255)

                                                                         tttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCAt                >  2:5226192/1‑93 (MQ=255)

                                                                         tttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCAt                <  2:2374502/93‑1 (MQ=255)

                                                                            gtgGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGaa             >  2:3330679/1‑93 (MQ=255)

                                                                               gTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAcc          <  2:4000248/93‑1 (MQ=255)

                                                                               gTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAcc          >  2:4244350/1‑93 (MQ=255)

                                                                                  gaagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaa       >  2:3918738/1‑93 (MQ=255)

                                                                                  gaagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAACaaa       >  1:5106407/1‑93 (MQ=255)

                                                                                    agaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaaa     <  1:1565522/93‑1 (MQ=255)

                                                                                    agaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaaa     <  2:3148052/93‑1 (MQ=255)

                                                                                     gaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAc    >  2:4992625/1‑93 (MQ=255)

                                                                                     gaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAc    <  2:4664063/93‑1 (MQ=255)

                                                                                     gaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAc    >  2:308907/1‑93 (MQ=255)

                                                                                      aTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACt   >  1:5252368/1‑93 (MQ=255)

                                                                                       tgctgcTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTa  >  2:2810992/1‑93 (MQ=255)

                                                                                       tgctgcTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTa  <  2:2716949/93‑1 (MQ=255)

                                                                                       tgctgcTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTa  <  1:4214215/93‑1 (MQ=255)

CGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973126‑3973305

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