BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,673,594	T→C	Y74C (TAC→TGC)	yphC ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,673,594	0	T	C	99.3%	437.5 / ‑6.6	138	Y74C (TAC→TGC)	yphC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (70/67); minor base A (1/0); total (71/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

AGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATT  >  NC_000913/2673504‑2673683

aGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGcat                                                                                         >  1:3121460/1‑93 (MQ=255)

aGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGcat                                                                                         >  1:2208910/1‑93 (MQ=255)

aGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGcat                                                                                         >  2:4058940/1‑93 (MQ=255)

  cACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATaa                                                                                       >  1:1120085/1‑93 (MQ=255)

   aCGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAAc                                                                                      <  1:2817383/93‑1 (MQ=255)

   aCGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAAc                                                                                      <  1:1308298/93‑1 (MQ=255)

    cgcgGTCGCCCTCTTTAAAATGGCGGAACCCTTGCCCCACCGCCCCAATCTGCCCCAACGGTTCATGACCGTTGATAAAGCCCGGGCATAAcc                                                                                     >  2:4091734/1‑92 (MQ=255)

      cgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGt                                                                                   <  2:4170715/93‑1 (MQ=255)

       ggTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGtt                                                                                  >  2:4145677/1‑93 (MQ=255)

        gTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGttt                                                                                 >  2:686637/1‑93 (MQ=255)

         tCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTa                                                                                <  1:298160/93‑1 (MQ=255)

           gCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATc                                                                              <  2:3873277/93‑1 (MQ=255)

           gCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATc                                                                              <  1:3971382/93‑1 (MQ=255)

             ccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCgg                                                                            >  2:777634/1‑93 (MQ=255)

              ctctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCggg                                                                           >  2:611056/1‑93 (MQ=255)

              ctctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCggg                                                                           >  2:1933785/1‑93 (MQ=255)

               tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          >  1:2131840/1‑93 (MQ=255)

               tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          >  2:30540/1‑93 (MQ=255)

                 tttAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgc                                                                        >  1:802771/1‑93 (MQ=255)

                 tttAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgc                                                                        >  2:2308481/1‑93 (MQ=255)

                 tttAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgc                                                                        >  1:1951149/1‑93 (MQ=255)

                  ttAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgcc                                                                       <  1:766667/93‑1 (MQ=255)

                    aaaaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgccgc                                                                     <  1:2239142/93‑1 (MQ=255)

                      aaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGGTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGAGccgccg                                                                   >  2:439052/1‑93 (MQ=255)

                        tGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCt                                                                 >  2:3867068/1‑93 (MQ=255)

                          gcggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCtgt                                                               <  1:3012152/93‑1 (MQ=255)

                          gcggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCtgt                                                               <  2:3087027/93‑1 (MQ=255)

                           cggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCtgtg                                                              <  1:246066/93‑1 (MQ=255)

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTgg                                                             <  1:2828424/93‑1 (MQ=255)

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTgg                                                             <  1:61529/93‑1 (MQ=255)

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTgg                                                             >  2:1296022/1‑93 (MQ=255)

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTgg                                                             >  2:3576562/1‑93 (MQ=255)

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTgg                                                             <  1:1685198/93‑1 (MQ=255)

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCACTGGCATAACGGTTTATCGGGTGGCGAAGCTGTgg                                                             >  1:4192997/1‑93 (MQ=255)

                             gcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGc                                                            <  1:2440548/93‑1 (MQ=255)

                                gCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACg                                                         >  2:2319495/1‑93 (MQ=255)

                                gCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACg                                                         >  2:1560327/1‑93 (MQ=255)

                                gCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACg                                                         <  1:2676541/93‑1 (MQ=255)

                                  cTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACggg                                                       >  2:4636616/1‑92 (MQ=255)

                                   ttGCCGCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtg                                                      <  2:2040900/93‑1 (MQ=255)

                                      ccccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtg                                                   >  2:892295/1‑93 (MQ=255)

                                       cccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtga                                                  <  1:3355218/93‑1 (MQ=255)

                                       cccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtga                                                  <  1:2695980/93‑1 (MQ=255)

                                        ccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgat                                                 >  1:3429847/1‑93 (MQ=255)

                                        ccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgat                                                 >  2:330439/1‑93 (MQ=255)

                                        ccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgat                                                 <  2:2150694/93‑1 (MQ=255)

                                        ccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgag                                                 >  2:3651168/1‑92 (MQ=255)

                                         cATCGCCACAATCTGCCCGCACGGTTCATGACCGTTTATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatg                                                >  1:4345908/1‑93 (MQ=255)

                                          aTCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatga                                               <  1:3533538/93‑1 (MQ=255)

                                           tCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTgatgat                                              >  1:232836/1‑93 (MQ=255)

                                           tCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTgatgat                                              <  1:967442/93‑1 (MQ=255)

                                            cGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGACgata                                             <  2:3699716/93‑1 (MQ=255)

                                              ccACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgataga                                           <  1:3841642/93‑1 (MQ=255)

                                               cacaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagat                                          <  2:2935648/93‑1 (MQ=255)

                                               cacaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagat                                          <  1:1631098/93‑1 (MQ=255)

                                               cacaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagat                                          >  2:3956448/1‑93 (MQ=255)

                                               cacaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagat                                          >  2:1255923/1‑93 (MQ=255)

                                                acaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagata                                         >  1:4327287/1‑93 (MQ=255)

                                                    tCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGt                                     <  1:4035638/93‑1 (MQ=255)

                                                    tCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGt                                     <  1:4674709/93‑1 (MQ=255)

                                                    tCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGt                                     <  1:789983/93‑1 (MQ=255)

                                                    tCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGt                                     >  2:2590821/1‑93 (MQ=255)

                                                      tGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGCTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGAgg                                   >  1:4180850/1‑93 (MQ=255)

                                                       gCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGa                                  <  2:3305716/93‑1 (MQ=255)

                                                       gCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGa                                  <  2:2948047/93‑1 (MQ=255)

                                                       gCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGa                                  <  1:4531982/93‑1 (MQ=255)

                                                         ccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACa                                <  2:2586140/93‑1 (MQ=255)

                                                         ccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACa                                <  1:3061331/93‑1 (MQ=255)

                                                         ccGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACa                                <  2:1772576/93‑1 (MQ=255)

                                                             aCGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGc                            >  1:4557088/1‑93 (MQ=255)

                                                               ggTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCtt                          >  2:4206016/1‑93 (MQ=255)

                                                                gTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTc                         <  2:4659825/93‑1 (MQ=255)

                                                                 ttCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTcc                        >  1:3620428/1‑93 (MQ=255)

                                                                  tCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCg                       >  1:3265344/1‑93 (MQ=255)

                                                                    aTGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCa                     >  2:3543334/1‑93 (MQ=255)

                                                                     tGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGAGTTGATGATAGATATAGTGGACATCGCTTCCGCaa                    >  1:1767341/1‑93 (MQ=255)

                                                                       aCCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAAt                  >  1:4006437/1‑93 (MQ=255)

                                                                       aCCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAAt                  >  1:209394/1‑93 (MQ=255)

                                                                        ccGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATc                 >  1:3782755/1‑93 (MQ=255)

                                                                        ccGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATc                 >  1:1614426/1‑93 (MQ=255)

                                                                         cGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATAGCTTCCGCAAATcc                >  1:3626570/1‑93 (MQ=255)

                                                                          gTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATccc               >  2:3375601/1‑93 (MQ=255)

                                                                          gTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATccc               <  2:1249747/93‑1 (MQ=255)

                                                                           ttGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCg              <  1:1103764/93‑1 (MQ=255)

                                                                           ttGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCg              >  2:3961228/1‑93 (MQ=255)

                                                                            tGATAAAGCCCTGGAATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgg             >  2:2214325/1‑93 (MQ=255)

                                                                             gATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgga            >  1:1898135/1‑93 (MQ=255)

                                                                             gATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCgga            <  1:1542630/93‑1 (MQ=255)

                                                                              aTAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggag           <  1:355389/93‑1 (MQ=255)

                                                                                aaaGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagga         <  2:2791490/93‑1 (MQ=255)

                                                                                aaaGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagga         <  2:4461217/93‑1 (MQ=255)

                                                                                 aaGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAt        <  2:87173/93‑1 (MQ=255)

                                                                                  aGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAtt       <  2:2745427/93‑1 (MQ=255)

                                                                                    cccTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTc     >  1:1407513/1‑93 (MQ=255)

                                                                                    cccTGGCATAACGGTGTATCGGGTTCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTc     <  1:4420796/93‑1 (MQ=255)

                                                                                     ccTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCa    <  1:4367710/93‑1 (MQ=255)

                                                                                      cTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCAt   >  1:2319085/1‑93 (MQ=255)

                                                                                      cTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCAt   <  1:1122620/93‑1 (MQ=255)

                                                                                      cTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCAt   >  2:1533823/1‑93 (MQ=255)

                                                                                      cTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCAGCAAATCCCGGAGGATTTCAt   >  2:523250/1‑93 (MQ=255)

AGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATT  >  NC_000913/2673504‑2673683

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