BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,008,990	C→T	R109C (CGT→TGT)	yigM →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,008,990	0	C	T	100.0%	1196.2 / NA	315	R109C (CGT→TGT)	yigM	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base T (149/166); total (149/166)

Only 100 of 320 total aligned reads displayed.

TATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTG  >  NC_000913/4008900‑4009078

tATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGGTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGt                                                                                        >  2:3573032/1‑93 (MQ=255)

 ttcgctACCGGACGGTTTCCGAGCGGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                       <  1:2717144/89‑1 (MQ=255)

 aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                       <  2:2071/93‑1 (MQ=255)

 aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                       >  1:3593770/1‑93 (MQ=255)

 aTCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTc                                                                                       >  1:3187523/1‑93 (MQ=255)

  tctctcCCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCt                                                                                      <  2:347017/93‑1 (MQ=255)

  tctcTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCt                                                                                      >  1:279202/1‑93 (MQ=255)

   ctctACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTg                                                                                     >  2:757994/1‑93 (MQ=255)

    tctACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgc                                                                                    >  2:2751882/1‑93 (MQ=255)

      tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                  >  1:1599916/1‑93 (MQ=255)

      tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                  >  2:232320/1‑93 (MQ=255)

      tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                  >  2:4060330/1‑93 (MQ=255)

      tACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTgcgc                                                                                  <  2:65761/93‑1 (MQ=255)

       tccTGACGGTTCCCGAGCTGCTGCTGTTCGCCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCt                                                                                 <  2:5848718/92‑1 (MQ=255)

       aCCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCt                                                                                 >  1:2835254/1‑93 (MQ=255)

       aCCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCt                                                                                 >  2:5495734/1‑93 (MQ=255)

        ccTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTg                                                                                <  1:5699393/93‑1 (MQ=255)

         cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                               <  1:4459390/93‑1 (MQ=255)

         cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                               <  1:1792802/93‑1 (MQ=255)

         cTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTgg                                                                               <  1:6181575/93‑1 (MQ=255)

          tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                              >  1:4492296/1‑93 (MQ=255)

          tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                              >  2:4870403/1‑93 (MQ=255)

          tGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTggg                                                                              <  2:3988424/93‑1 (MQ=255)

            aCGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGc                                                                            >  2:2183285/1‑93 (MQ=255)

            aCGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGc                                                                            <  2:3696617/93‑1 (MQ=255)

              ggTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTa                                                                          >  1:6358140/1‑93 (MQ=255)

              ggTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTa                                                                          <  1:2276863/93‑1 (MQ=255)

              ggTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTa                                                                          <  1:48106/93‑1 (MQ=255)

               gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACCTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                         >  1:599214/1‑93 (MQ=255)

               gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                         <  2:5931820/93‑1 (MQ=255)

               gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                         <  2:6026149/93‑1 (MQ=255)

               gTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTAAATCACGCTGATTTATGACATCATGCGTAAGCGCTGTCTGCGCTGGGGCTAt                                                                         >  2:3209603/1‑93 (MQ=255)

                  tCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGcc                                                                      <  2:5707048/93‑1 (MQ=255)

                  tCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGcc                                                                      >  1:941748/1‑93 (MQ=255)

                    cGAGCTGCTGCTGTTCACCGTTCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGGCATGCCtt                                                                    <  1:6608478/93‑1 (MQ=255)

                    cGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCtt                                                                    <  2:3849857/93‑1 (MQ=255)

                    cGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCtt                                                                    <  2:3276944/93‑1 (MQ=255)

                     caGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCttt                                                                   <  2:4802374/92‑1 (MQ=255)

                       gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                 <  1:6407770/93‑1 (MQ=255)

                       gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                 >  1:5903957/1‑93 (MQ=255)

                       gctgctgctGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAg                                                                 <  2:5562484/93‑1 (MQ=255)

                           ctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGcc                                                             >  2:5675806/1‑93 (MQ=255)

                           ctgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGcc                                                             <  2:2639727/93‑1 (MQ=255)

                            tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                            <  1:5171559/93‑1 (MQ=255)

                            tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                            <  1:6225358/93‑1 (MQ=255)

                            tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                            <  1:4223072/93‑1 (MQ=255)

                            tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                            >  1:3076738/1‑93 (MQ=255)

                            tgctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCt                                                            <  1:1145490/93‑1 (MQ=255)

                              ctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTg                                                          >  1:1354197/1‑93 (MQ=255)

                              ctgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTg                                                          <  2:257608/93‑1 (MQ=255)

                               tgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGc                                                         <  1:5281537/93‑1 (MQ=255)

                               tgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGc                                                         >  2:4238753/1‑93 (MQ=255)

                               tgTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGc                                                         <  1:2549331/93‑1 (MQ=255)

                                cttCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCt                                                        <  1:5515619/92‑1 (MQ=255)

                                       gagCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTg                                                 >  1:4203560/3‑93 (MQ=255)

                                        tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                >  1:3530250/1‑93 (MQ=255)

                                        tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                <  2:5242138/93‑1 (MQ=255)

                                        tGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGa                                                >  2:6641932/1‑93 (MQ=255)

                                         gCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAt                                               >  1:1279507/1‑93 (MQ=255)

                                         gCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGAt                                               <  1:329100/93‑1 (MQ=255)

                                            gACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTgg                                            >  1:4015140/1‑93 (MQ=255)

                                             aCGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGt                                           <  1:5817612/93‑1 (MQ=255)

                                             aCGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGt                                           <  1:4951038/93‑1 (MQ=255)

                                                   ctctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCggg                                     >  2:1876349/1‑93 (MQ=255)

                                                     ctACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGACGGTGATTGGTGCCGGGat                                   <  1:6540407/93‑1 (MQ=255)

                                                      tACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGatt                                  >  1:2013745/1‑93 (MQ=255)

                                                      tACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGatt                                  >  1:214461/1‑93 (MQ=255)

                                                      tACATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGatt                                  <  1:3023326/93‑1 (MQ=255)

                                                       aCATCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAgcg                                                               <  2:6365555/63‑1 (MQ=255)

                                                         aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                               <  2:1767281/93‑1 (MQ=255)

                                                         aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                               >  2:1713908/1‑93 (MQ=255)

                                                         aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                               <  1:3766947/93‑1 (MQ=255)

                                                         aTCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGattatt                               >  1:3680007/1‑93 (MQ=255)

                                                          tCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTc                              <  2:47136/93‑1 (MQ=255)

                                                          tCACGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTc                              >  1:3458125/1‑93 (MQ=255)

                                                            aCGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGc                            >  1:6126586/1‑93 (MQ=255)

                                                            aCGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGc                            >  1:2948501/1‑93 (MQ=255)

                                                             cGCTGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCt                           >  1:1012407/1‑93 (MQ=255)

                                                                tGATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATg                        <  1:3834273/93‑1 (MQ=255)

                                                                 gATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGa                       >  1:3897997/1‑93 (MQ=255)

                                                                 gATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGa                       >  1:1496794/1‑93 (MQ=255)

                                                                 gATTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGa                       <  2:2504807/93‑1 (MQ=255)

                                                                  aTTTATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGAt                      <  2:1229277/93‑1 (MQ=255)

                                                                   tttATGACATCATGAGTAAGCGCTGGCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATc                     <  1:110932/93‑1 (MQ=255)

                                                                    ttATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCa                    <  2:3793782/93‑1 (MQ=255)

                                                                     tATGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAg                   <  2:3000565/93‑1 (MQ=255)

                                                                      aTGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAgg                  <  1:833984/93‑1 (MQ=255)

                                                                      aTGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAgg                  <  1:3411161/93‑1 (MQ=255)

                                                                       tGACATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGt                 <  2:4670652/93‑1 (MQ=255)

                                                                         aCATCATGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGGGCCGGGATTATTCGCTATGATCAGGtaa               >  1:3419946/1‑91 (MQ=255)

                                                                            tcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCg            >  2:1433995/1‑93 (MQ=255)

                                                                            tcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCg            <  1:4260142/93‑1 (MQ=255)

                                                                            tcatGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCg            >  2:817277/1‑93 (MQ=255)

                                                                             catGAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGa           <  2:3054280/93‑1 (MQ=255)

                                                                                gAGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCa        >  2:1543773/1‑93 (MQ=255)

                                                                                 aGTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAt       <  2:5340227/93‑1 (MQ=255)

                                                                                  gTAAGCGCTGTCTGCGCTGGGTCTATGCCTTTAGCGCCTTGCTGGCGATGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAtt      <  1:6720826/93‑1 (MQ=255)

                                                                                  gTAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAtt      <  2:5431491/93‑1 (MQ=255)

                                                                                   tAAGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCAttt     <  1:5094735/93‑1 (MQ=255)

                                                                                    aaGCGCTGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTc    <  2:5012908/93‑1 (MQ=255)

TATCTCTACCTGACGGTTTCCGAGCTGCTGCTGTTCACCGTGCTGACGCCGCTCTACATCACGCTGATTTATGACATCATGAGTAAGCGCCGTCTGCGCTGGGGCTATGCCTTTAGCGCCTTGCTGGCGGTGATTGGTGCCGGGATTATTCGCTATGATCAGGTCACCGACCATTTCTG  >  NC_000913/4008900‑4009078

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