BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,137,219	C→T	G24R (GGG→AGG)	talC ←	putative transaldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,137,219	0	C	T	98.9%	1331.8 / ‑3.7	374	G24R (GGG→AGG)	talC	putative transaldolase
Reads supporting (aligned to +/- strand): ref base C (4/0); new base T (171/199); total (175/199)
Fisher's exact test for biased strand distribution p-value = 4.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 389 total aligned reads displayed.

CCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTC  >  NC_000913/4137129‑4137308

ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTgg                                                                                         <  2:7525235/93‑1 (MQ=255)

  tcatcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCa                                                                                       >  1:5350866/1‑93 (MQ=255)

  tcatcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCa                                                                                       >  2:4888518/1‑93 (MQ=255)

   catcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCaa                                                                                      >  2:5145702/1‑93 (MQ=255)

      caCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgg                                                                                   <  1:7641824/93‑1 (MQ=255)

        ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATgggg                                                                                 >  2:4164856/1‑93 (MQ=255)

         cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGa                                                                                <  1:4815100/93‑1 (MQ=255)

           aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                              >  1:7065092/1‑93 (MQ=255)

           aTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAat                                                                              <  2:4434559/93‑1 (MQ=255)

            tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                             <  2:2724906/93‑1 (MQ=255)

            tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                             <  2:3450995/93‑1 (MQ=255)

            tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAata                                                                             <  1:836305/93‑1 (MQ=255)

            tCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTAGCAATGGGGAata                                                                             >  1:1325168/1‑93 (MQ=255)

             cGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAatat                                                                            <  1:3795499/93‑1 (MQ=255)

              gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           >  2:6769751/1‑93 (MQ=255)

              gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATg                                                                           <  1:1971711/93‑1 (MQ=255)

               ctctttGCAGACGCGGCAGCACTTCCCATATGGACTCCGTGATGCCAGCGATAATGCTCGGGTTGCTTGTCACCCTGGCAATGGGGAATATGc                                                                          <  2:3764743/90‑1 (MQ=255)

               cTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGc                                                                          <  2:7615952/93‑1 (MQ=255)

                  tttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTg                                                                       <  1:1388284/93‑1 (MQ=255)

                     gCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCa                                                                    <  2:3155302/93‑1 (MQ=255)

                       aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGa                                                                  >  1:6877107/1‑93 (MQ=255)

                       aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGCTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGa                                                                  <  2:2612138/93‑1 (MQ=255)

                        gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                 <  1:2091343/93‑1 (MQ=255)

                        gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGAc                                                                 >  1:1693679/1‑93 (MQ=255)

                         aCGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACg                                                                <  1:7345179/93‑1 (MQ=255)

                          cgcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGt                                                               >  1:5589134/1‑93 (MQ=255)

                           gcgGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGtt                                                              >  1:502709/1‑93 (MQ=255)

                              gcagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGa                                                           >  1:6909774/1‑93 (MQ=255)

                               cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGAc                                                          <  1:4318011/93‑1 (MQ=255)

                                 gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                                        <  2:2440686/93‑1 (MQ=255)

                                 gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACtt                                                        >  1:4424094/1‑93 (MQ=255)

                                  caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                                       <  2:4383608/93‑1 (MQ=255)

                                  caCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTc                                                       <  1:6259206/93‑1 (MQ=255)

                                    cTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                                     >  1:3962284/1‑93 (MQ=255)

                                    cTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTg                                                     <  2:5596410/93‑1 (MQ=255)

                                     ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                    >  2:7725903/1‑93 (MQ=255)

                                     ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                    <  2:742897/93‑1 (MQ=255)

                                     ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                    >  2:7247502/1‑93 (MQ=255)

                                     ttCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGc                                                    >  1:5833731/1‑93 (MQ=255)

                                      tCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCg                                                   >  1:2545300/1‑93 (MQ=255)

                                       cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGa                                                  >  2:5730337/1‑93 (MQ=255)

                                       cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGa                                                  <  1:5679585/93‑1 (MQ=255)

                                           tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                              <  1:2793612/93‑1 (MQ=255)

                                           tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                              <  1:4614355/93‑1 (MQ=255)

                                           tatGGACTCCTTGCTGGCAGCGATAATGCTCGGGCTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGCTCGACTTCTGCGACGtt                                              >  1:3126590/1‑93 (MQ=255)

                                           tatGCACTCCTTGCTGGCAGCGATAATGCGCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGtt                                              <  1:6130891/93‑1 (MQ=255)

                                            atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                             >  1:4724575/1‑93 (MQ=255)

                                            atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTa                                             <  1:1825129/93‑1 (MQ=255)

                                              ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                           >  1:7492049/1‑93 (MQ=255)

                                              ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGc                                           <  2:393248/93‑1 (MQ=255)

                                               gACTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCg                                          >  2:6333007/1‑93 (MQ=255)

                                                aCTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCgg                                         >  1:4089491/1‑93 (MQ=255)

                                                 cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        <  1:2177032/93‑1 (MQ=255)

                                                 cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        <  1:6409970/93‑1 (MQ=255)

                                                 cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        <  1:242617/93‑1 (MQ=255)

                                                 cTCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgt                                        <  1:6752284/93‑1 (MQ=255)

                                                  tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                                       >  1:5905546/1‑93 (MQ=255)

                                                  tCCTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGgtg                                       >  1:3043734/1‑93 (MQ=255)

                                                    cTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTc                                     >  1:3885629/1‑93 (MQ=255)

                                                    cTTGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTc                                     >  2:6517469/1‑93 (MQ=255)

                                                      tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                                   <  1:1426360/93‑1 (MQ=255)

                                                      tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                                   <  2:1915228/93‑1 (MQ=255)

                                                      tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                                   <  1:3357895/93‑1 (MQ=255)

                                                      tGCTGGCAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCa                                   >  2:3267518/1‑93 (MQ=255)

                                                           gCAGCGATAATGCTCGGTTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATAc                              <  2:3531943/93‑1 (MQ=255)

                                                            cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                             <  1:5540019/93‑1 (MQ=255)

                                                            cAGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACa                             <  1:2949975/93‑1 (MQ=255)

                                                             aGCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGCCGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAg                            <  2:2999164/93‑1 (MQ=255)

                                                              gCGATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGt                           <  2:5243001/93‑1 (MQ=255)

                                                               cGATAATGCTCGGGTTGGTTGNCNCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGtt                          <  1:1573905/93‑1 (MQ=255)

                                                                gATAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTc                         <  2:6098218/93‑1 (MQ=255)

                                                                  tAATGCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCa                       >  2:6644710/1‑93 (MQ=255)

                                                                      gCTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgat                   <  2:2863807/93‑1 (MQ=255)

                                                                       cTCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCatgatg                  <  2:2800751/93‑1 (MQ=255)

                                                                       cTCGGGTTGGTTGACATCCTGGCAATGGGGAATATGCGTGCCCGACGTTCGACTTCTGCGACGTTAGCGGCGTCCAGACACAGTTTCatgatg                  >  2:1066900/1‑93 (MQ=255)

                                                                        tCGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGt                 >  1:901731/1‑93 (MQ=255)

                                                                         cGGGTTGGTTGTCACCCTGGGAATGGGGAATATGCGTGCCAGACGTTCGATTTCTGTGACGTTAGAGGTGTCGCGATACAGTTACATGATGtt                >  2:1632509/1‑93 (MQ=255)

                                                                         cGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt                <  2:6940007/93‑1 (MQ=255)

                                                                         cGGGTTGGTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGtt                <  2:1164595/93‑1 (MQ=255)

                                                                          gggTTAGTGATCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGttt               <  1:6558153/93‑1 (MQ=255)

                                                                           ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              >  2:7086081/1‑93 (MQ=255)

                                                                           ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              <  1:1140962/93‑1 (MQ=255)

                                                                           ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              >  1:1650044/1‑93 (MQ=255)

                                                                           ggttggttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTc              <  2:2187027/93‑1 (MQ=255)

                                                                            gttggttgTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTcc             <  2:5810108/93‑1 (MQ=255)

                                                                             ttggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt            <  2:1739321/90‑1 (MQ=255)

                                                                             ttggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCt            <  2:4022680/90‑1 (MQ=255)

                                                                               ggTTGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCttt          >  2:4626998/2‑93 (MQ=255)

                                                                                 ttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTa        >  2:1456942/1‑93 (MQ=255)

                                                                                 ttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTa        <  1:6711263/93‑1 (MQ=255)

                                                                                 ttGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTAGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTa        <  1:3225169/93‑1 (MQ=255)

                                                                                  tGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAc       <  1:7811488/93‑1 (MQ=255)

                                                                                  tGTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAc       <  1:5086938/93‑1 (MQ=255)

                                                                                   ggcACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  1:1050754/91‑1 (MQ=255)

                                                                                   gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  2:5993173/93‑1 (MQ=255)

                                                                                   gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  2:5425932/93‑1 (MQ=255)

                                                                                   gTCACCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAca      <  1:1009055/93‑1 (MQ=255)

                                                                                      aCCCTGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagtt   >  2:7930046/1‑93 (MQ=255)

                                                                                       ccctGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttc  <  2:1472111/93‑1 (MQ=255)

                                                                                       ccctGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTAcagttc  >  1:3521963/1‑93 (MQ=255)

CCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTC  >  NC_000913/4137129‑4137308

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