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%	908.0 / NA	241	A153A (GCC→GCT)	wecE	putative regulator
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (118/123); total (118/123)

Only 100 of 246 total aligned reads displayed.

CGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973126‑3973305

cGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGAtgctgc                                                                                         >  2:879591/1‑93 (MQ=255)

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

     cGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTACTCAgg                                                                                    >  1:3672725/1‑93 (MQ=255)

      ggTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAggg                                                                                   <  2:6072956/93‑1 (MQ=255)

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

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

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

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

        tCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCg                                                                                 <  2:1308710/93‑1 (MQ=255)

         ccATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGt                                                                                >  2:5606006/1‑93 (MQ=255)

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

          cATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTg                                                                               >  1:534020/1‑93 (MQ=255)

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

           aTTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGa                                                                              <  2:4401301/93‑1 (MQ=255)

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

             tACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATg                                                                            <  2:3968949/93‑1 (MQ=255)

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

                gcgGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTcc                                                                         <  1:2090930/93‑1 (MQ=255)

                 cgGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCa                                                                        <  2:1587032/93‑1 (MQ=255)

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

                    gtgtGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACtt                                                                     <  1:13262/93‑1 (MQ=255)

                    gtgtGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACtt                                                                     <  1:176088/93‑1 (MQ=255)

                      gtgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTAc                                                                   >  1:5030334/1‑93 (MQ=255)

                        ggCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACaa                                                                 <  2:1661567/93‑1 (MQ=255)

                         gCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACaaa                                                                <  2:6338778/93‑1 (MQ=255)

                            tgagAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             >  1:3071853/4‑93 (MQ=255)

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

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

                             gCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGc                                                            <  2:2474075/93‑1 (MQ=255)

                             gCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGc                                                            <  1:501992/93‑1 (MQ=255)

                              cGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCg                                                           <  1:524540/93‑1 (MQ=255)

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

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

                                aaaTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTg                                                         >  1:4575388/1‑93 (MQ=255)

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

                                  aTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCa                                                       >  2:361465/1‑93 (MQ=255)

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

                                     gACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTg                                                    >  2:2267573/1‑93 (MQ=255)

                                     gACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTg                                                    >  2:3495549/1‑93 (MQ=255)

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

                                        acCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGa                                                 >  1:2916/1‑93 (MQ=255)

                                           aTTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAAcc                                              >  2:4044809/1‑93 (MQ=255)

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

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

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

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

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            >  1:4670179/1‑93 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  2:2819634/93‑1 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  1:3699820/93‑1 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  1:5658128/93‑1 (MQ=255)

                                             tATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAt                                            <  2:3733593/93‑1 (MQ=255)

                                              aTGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCAtt                                           >  1:4130097/1‑93 (MQ=255)

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

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

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

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

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

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

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

                                                  ggTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAACATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTc                                       <  1:1903256/92‑1 (MQ=255)

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

                                                  cGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTc                                       <  1:1364651/93‑1 (MQ=255)

                                                   gTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGTAACCATTGGTCa                                      <  2:256397/93‑1 (MQ=255)

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

                                                       gCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATAtt                                  <  2:298668/93‑1 (MQ=255)

                                                       gCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATAtt                                  <  1:3158462/93‑1 (MQ=255)

                                                        cGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTg                                 <  2:4133302/93‑1 (MQ=255)

                                                          aaaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgc                               <  1:6329056/93‑1 (MQ=255)

                                                           aaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgct                              >  1:6068628/1‑93 (MQ=255)

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

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

                                                             aaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgc                            <  1:2419584/93‑1 (MQ=255)

                                                             aaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgc                            >  1:3388550/1‑93 (MQ=255)

                                                             aaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgc                            >  1:5204741/1‑93 (MQ=255)

                                                              aGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgctgct                           <  1:2126863/93‑1 (MQ=255)

                                                               gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGAGTGCACTGGGAACCATTGGTCATATTGGCTGCtt                          >  2:1697000/1‑93 (MQ=255)

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

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

                                                                  tAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTaa                       >  1:4818273/1‑92 (MQ=255)

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

                                                                   aaTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGc                      >  2:2091807/1‑93 (MQ=255)

                                                                   aaTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGc                      >  2:2013469/1‑93 (MQ=255)

                                                                    aTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCt                     <  1:1194866/93‑1 (MQ=255)

                                                                    aTTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCt                     <  2:2895185/93‑1 (MQ=255)

                                                                     tttgtttgTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCtt                    >  1:1086964/1‑93 (MQ=255)

                                                                     tttgtttgTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCtt                    <  1:1475013/93‑1 (MQ=255)

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

                                                                          ttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATg               >  2:3789924/1‑93 (MQ=255)

                                                                          ttgtGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATg               >  1:3688898/1‑93 (MQ=255)

                                                                             tgGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGaaa            <  1:276113/93‑1 (MQ=255)

                                                                             tgGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGaaa            <  1:1794379/93‑1 (MQ=255)

                                                                              ggTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAc           <  1:1968956/93‑1 (MQ=255)

                                                                              ggTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAc           <  2:4441883/93‑1 (MQ=255)

                                                                                tAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCa         >  1:2686396/1‑93 (MQ=255)

                                                                                 agaagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaa        <  2:109548/93‑1 (MQ=255)

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

                                                                                   aagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaa      >  2:3377221/1‑93 (MQ=255)

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

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

                                                                                       tgctgcTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAATa  >  2:2177503/1‑93 (MQ=255)

CGTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973126‑3973305

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: `‑`