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

Only 100 of 233 total aligned reads displayed.

GTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973127‑3973305

 tGCCGTGCCATTACGCGGGTGTGGCTTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTc                                                                                       <  1:2651894/93‑1 (MQ=255)

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

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

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

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

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

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

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

      gTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGc                                                                                  <  2:1141386/93‑1 (MQ=255)

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

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

           ttACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGAt                                                                             >  1:1626444/1‑93 (MQ=255)

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

             aCGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGt                                                                           <  2:1283299/93‑1 (MQ=255)

                cgGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGCTGCTGCTCAGGGCGTGATGTccc                                                                        >  2:1823696/1‑92 (MQ=255)

                 gggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCAc                                                                       <  2:193076/93‑1 (MQ=255)

                 gggTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGAAGTCCAc                                                                       >  2:2435544/1‑93 (MQ=255)

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

                    tgtgGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTa                                                                    <  2:1724808/93‑1 (MQ=255)

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

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

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

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

                        gCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACaaa                                                                >  2:836949/1‑93 (MQ=255)

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

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

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

                           tGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             <  1:1979653/93‑1 (MQ=255)

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

                           tGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAggg                                                             <  2:1088985/93‑1 (MQ=255)

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

                              gAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCgg                                                          >  1:2418046/1‑92 (MQ=255)

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

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

                                 aTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGCTGTCCACTTACAAAGGGCGTGCa                                                       >  2:2030493/1‑93 (MQ=255)

                                 aTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCa                                                       <  2:1109584/93‑1 (MQ=255)

                                 aTGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCa                                                       <  1:685209/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                         cATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAAc                                               >  2:1009841/1‑93 (MQ=255)

                                         cATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAAc                                               >  1:1082809/1‑93 (MQ=255)

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

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

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

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

                                           atATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAAGCa                                             <  2:1499512/92‑1 (MQ=255)

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

                                                  gTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCa                                      <  2:574301/93‑1 (MQ=255)

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

                                                     ggCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCatat                                   <  1:1071729/93‑1 (MQ=255)

                                                     ggCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCatat                                   <  2:1261912/93‑1 (MQ=255)

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

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

                                                        gAAAAAGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTgg                                >  1:610157/1‑93 (MQ=255)

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

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

                                                          aaaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGgct                              <  2:4052638/93‑1 (MQ=255)

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

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

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

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

                                                           aaaGCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGGCATATTGgctg                             >  1:3707046/1‑93 (MQ=255)

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

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

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

                                                              gCATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACAGGGAACCATTGGTCATATTGGCTGCtt                          <  2:3211452/93‑1 (MQ=255)

                                                               cATAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCttt                         >  1:3632877/1‑93 (MQ=255)

                                                               cATAATTTGTTTGTGGGAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCttt                         >  2:1435290/1‑93 (MQ=255)

                                                                 tAATTTGTTTGTGGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAg                       >  1:1838079/1‑93 (MQ=255)

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

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

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

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

                                                                       gtttgtGGTAGAAGATGCTGCTCAGGGCGTGGTGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCa                 >  2:2918495/1‑93 (MQ=255)

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

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

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

                                                                            tgGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGaaa            >  1:155007/1‑93 (MQ=255)

                                                                            tgGTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGaaa            >  2:709686/1‑93 (MQ=255)

                                                                             ggTAGAAGTTGTGGCTAAGGGCGTGATGTCCACTTACAAGGGGCGTGCGCT‑GGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAAcc          <  1:3081199/93‑1 (MQ=255)

                                                                             ggTAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGGTTCCATTAAAc           >  1:35492/1‑93 (MQ=255)

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

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

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

                                                                               tAGAAGATGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGGTTTAGCTTCCATGAAACCa         >  1:403777/1‑93 (MQ=255)

                                                                                 gaagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTCGCTTCCATGAAACCaaa       >  1:813971/1‑93 (MQ=255)

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

                                                                                  aagaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaa      <  1:2092158/93‑1 (MQ=255)

                                                                                   agaTGCTGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCaaaaa     >  2:4124174/1‑93 (MQ=255)

                                                                                    ggtgctgcTCAGCGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAAc    <  2:4528094/91‑1 (MQ=255)

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

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

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

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

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

GTGCCGGTCCATTACGCGGGTGTGGCCTGCGAAATGGACACCATTATGGCGTTGGCGAAAAAGCATAATTTGTTTGTGGTAGAAGATGCCGCTCAGGGCGTGATGTCCACTTACAAAGGGCGTGCACTGGGAACCATTGGTCATATTGGCTGCTTTAGCTTCCATGAAACCAAAAACTA  >  NC_000913/3973127‑3973305

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