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.6%	863.2 / NA	262	Y74C (TAC→TGC)	yphC	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (115/146); total (115/147)
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 272 total aligned reads displayed.

AGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATT  >  NC_000913/2673504‑2673683

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

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

  cACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATaa                                                                                       <  2:1147828/93‑1 (MQ=255)

   aCGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAAc                                                                                      >  2:2661557/1‑93 (MQ=255)

    cgcgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGGCCTGGCATAACg                                                                                     >  1:3400790/1‑93 (MQ=255)

    cgcgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACg                                                                                     <  2:3631360/93‑1 (MQ=255)

    cgcgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACg                                                                                     <  1:4155845/93‑1 (MQ=255)

    cgcgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACg                                                                                     <  2:5838906/93‑1 (MQ=255)

    cgcgGTCGCCCTCTTAAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACg                                                                                     <  2:4572660/93‑1 (MQ=255)

      cgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGt                                                                                   >  1:4485952/1‑93 (MQ=255)

      cgGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGt                                                                                   >  1:4137133/1‑93 (MQ=255)

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

       ggTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGtt                                                                                  <  2:3336358/93‑1 (MQ=255)

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

       ggTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTCATAAAGCCCTAGCATAACgggt                                                                                  >  2:1930898/1‑91 (MQ=255)

        gTCGCCGTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGttt                                                                                 <  2:2833524/93‑1 (MQ=255)

        gTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGttt                                                                                 <  1:5226839/93‑1 (MQ=255)

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

          cGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTAt                                                                               >  2:4272717/1‑93 (MQ=255)

           gCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTCTATc                                                                              >  2:1562200/1‑93 (MQ=255)

            cccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCg                                                                             >  2:2212512/1‑93 (MQ=255)

            cccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCg                                                                             >  1:4068277/1‑93 (MQ=255)

            cccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCg                                                                             >  1:1659333/1‑93 (MQ=255)

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

             ccTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCgg                                                                            <  1:6506358/93‑1 (MQ=255)

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

              ctctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCggg                                                                           <  2:5139385/93‑1 (MQ=255)

              ctctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCggg                                                                           <  2:1020153/93‑1 (MQ=255)

               tctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          <  1:1674873/93‑1 (MQ=255)

               cctTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGt                                                                          <  1:4682348/92‑1 (MQ=255)

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

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

                     aaaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTgccgcc                                                                    <  1:171714/93‑1 (MQ=255)

                     aaaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgccgcc                                                                    >  2:3813742/1‑93 (MQ=255)

                     aaaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTgccgcc                                                                    >  1:2134376/1‑93 (MQ=255)

                      aaTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGccgccg                                                                   >  2:5479556/1‑93 (MQ=255)

                       aTGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCcgccgc                                                                  <  1:5648932/93‑1 (MQ=255)

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

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

                            ggcAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGCTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTgg                                                             >  1:6023939/1‑93 (MQ=255)

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

                              cAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCa                                                           <  1:1998752/93‑1 (MQ=255)

                               aGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCAc                                                          >  1:2372127/1‑93 (MQ=255)

                                  cTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgt                                                       <  1:2312470/93‑1 (MQ=255)

                                  cTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCACTGGCATAACGTTTTATCGGGTGCCGCCTCTGTGGCACGgt                                                       >  2:2712244/1‑93 (MQ=255)

                                   ttGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtg                                                      >  1:671501/1‑93 (MQ=255)

                                    tGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGgtgt                                                     <  1:2113552/93‑1 (MQ=255)

                                     gCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGtt                                                    <  2:204949/93‑1 (MQ=255)

                                      ccccATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtg                                                   <  1:1078913/93‑1 (MQ=255)

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

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

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

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

                                         cATCGCCACAATCTGCCTGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatg                                                <  2:4551735/93‑1 (MQ=255)

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

                                          aTCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTtgatga                                               >  2:5921748/1‑93 (MQ=255)

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

                                            cGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgata                                             <  1:6356472/93‑1 (MQ=255)

                                                acaATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagata                                         <  2:4814675/93‑1 (MQ=255)

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

                                                acaATCTCCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATgatagata                                         <  1:1304917/93‑1 (MQ=255)

                                                  aaTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGAtata                                       <  2:4336266/93‑1 (MQ=255)

                                                  aaTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGAtata                                       <  2:2942278/93‑1 (MQ=255)

                                                  aaTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACAGTGTTGATGATAGAtata                                       <  2:3006497/93‑1 (MQ=255)

                                                   aTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGatatcg                                      >  1:3271695/1‑91 (MQ=255)

                                                   aTCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAg                                      <  2:406538/93‑1 (MQ=255)

                                                      tGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTgg                                   <  2:3126899/93‑1 (MQ=255)

                                                           gCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATc                              <  1:2534871/93‑1 (MQ=255)

                                                           gCACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATc                              >  2:3466348/1‑93 (MQ=255)

                                                            cACGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTTGACATCg                             >  2:6298712/1‑93 (MQ=255)

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

                                                              cGGTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCt                           <  2:6041387/93‑1 (MQ=255)

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

                                                               ggTTCATGACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCtt                          <  1:5298235/93‑1 (MQ=255)

                                                                gTTCATGACCGTTGATAAAGACCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTc                         >  1:6453498/1‑93 (MQ=255)

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

                                                                      gACCGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCaaa                   <  1:3598232/93‑1 (MQ=255)

                                                                         cGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATcc                <  2:4896930/93‑1 (MQ=255)

                                                                         cGTTGATAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATcc                <  2:20198/93‑1 (MQ=255)

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

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

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

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

                                                                              aTAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggag           >  2:1036264/1‑93 (MQ=255)

                                                                               tAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagg          >  1:5502882/1‑93 (MQ=255)

                                                                               tAAAGCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagg          <  1:1375037/93‑1 (MQ=255)

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

                                                                                aaaGCCCTGGCATAACGGTTTATCGGGTCCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCggagga         <  2:874710/93‑1 (MQ=255)

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

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

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

                                                                                   gCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAttt      <  2:353536/93‑1 (MQ=255)

                                                                                   gCCCTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGAttt      <  1:4886975/93‑1 (MQ=255)

                                                                                    cgcTGGCATAATGGTTTATCGGGTTCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTc     <  1:894172/91‑1 (MQ=255)

                                                                                    cccTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTc     <  2:2383561/93‑1 (MQ=255)

                                                                                    cccTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTc     <  2:653811/93‑1 (MQ=255)

                                                                                     ccTGGCATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTTCGCAAATCCCGGAGGATTTCa    >  1:5157403/1‑93 (MQ=255)

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

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

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

AGCACGCGGTCGCCCTCTTTAAAATGGCGGCAGCCTTGCCCCATCGCCACAATCTGCCCGCACGGTTCATGACCGTTGATAAAGCCCTGGTATAACGGTTTATCGGGTGCCGCCGCTGTGGCACGGTGTTGATGATAGATATAGTGGACATCGCTTCCGCAAATCCCGGAGGATTTCATT  >  NC_000913/2673504‑2673683

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