BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,101,362	T→A	V111D (GTT→GAT)	mutY →	adenine glycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,101,362	0	T	A	99.6%	939.0 / ‑3.9	240	V111D (GTT→GAT)	mutY	adenine glycosylase
Reads supporting (aligned to +/- strand): ref base T (1/0); new base A (106/133); total (107/133)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

Only 100 of 241 total aligned reads displayed.

GCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACG  >  NC_000913/3101272‑3101451

 gtATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATgc                                                                                        <  2:4918644/92‑1 (MQ=255)

 cTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATgc                                                                                        <  2:4237820/93‑1 (MQ=255)

  tATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATgca                                                                                       <  2:402634/93‑1 (MQ=255)

    ttACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATgcagc                                                                                     <  1:6020985/93‑1 (MQ=255)

    ttACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATgcagc                                                                                     <  2:3389401/93‑1 (MQ=255)

       cGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACt                                                                                  <  1:1471055/93‑1 (MQ=255)

       cGCCCGCGCGCGAAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACt                                                                                  <  2:1780701/93‑1 (MQ=255)

        gCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTg                                                                                 <  1:5552983/93‑1 (MQ=255)

        gCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTg                                                                                 >  1:4418261/1‑93 (MQ=255)

         cccGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGc                                                                                >  2:2916755/1‑93 (MQ=255)

          ccGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGcc                                                                               <  1:3992093/93‑1 (MQ=255)

          ccGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGcc                                                                               <  1:6106251/93‑1 (MQ=255)

           cgcgcgcgCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcg                                                                              >  1:858309/1‑93 (MQ=255)

           cgcgcgcgCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcg                                                                              >  2:63152/1‑93 (MQ=255)

           cgcgcgcgCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcg                                                                              <  2:1714850/93‑1 (MQ=255)

             cgcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcggg                                                                            >  1:652547/1‑93 (MQ=255)

             cgcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcggg                                                                            <  1:3739588/93‑1 (MQ=255)

              gcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggc                                                                           >  2:1400452/1‑93 (MQ=255)

              gcgcgcAATCGGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAATCCTTGGAGGAAGATGCAGCACTGCcgggc                                                                           <  1:1593390/93‑1 (MQ=255)

               cgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcg                                                                          >  2:4282055/1‑93 (MQ=255)

                gcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgt                                                                         >  1:3879673/1‑93 (MQ=255)

                gcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgt                                                                         >  1:4193384/1‑93 (MQ=255)

                  gcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcg                                                                       <  2:4104040/93‑1 (MQ=255)

                  gcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcg                                                                       <  1:3893237/93‑1 (MQ=255)

                   cAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgg                                                                      >  1:3280632/1‑93 (MQ=255)

                   cAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgg                                                                      >  1:2310043/1‑93 (MQ=255)

                    aaTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtgggg                                                                     >  1:1737801/1‑93 (MQ=255)

                    aaTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcggg                                                                     >  1:5354921/1‑93 (MQ=255)

                    aaTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcggg                                                                     >  2:88808/1‑93 (MQ=255)

                     aTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgggc                                                                    >  1:4329608/1‑93 (MQ=255)

                     aTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgggc                                                                    <  2:5416173/93‑1 (MQ=255)

                     aTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgggc                                                                    >  1:4993646/1‑93 (MQ=255)

                      tCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgggcg                                                                   >  2:4997911/1‑93 (MQ=255)

                       cTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCcgggcgtcgggcgt                                                                  <  2:5866782/93‑1 (MQ=255)

                           aTAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCa                                                              >  1:1036054/1‑93 (MQ=255)

                            tAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCAc                                                             <  1:2717660/93‑1 (MQ=255)

                             aaaGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCAcc                                                            >  1:1756785/1‑93 (MQ=255)

                             aaaGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCAcc                                                            >  1:3854654/1‑93 (MQ=255)

                                gCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCa                                                         >  1:3157963/1‑93 (MQ=255)

                                 cGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAg                                                        <  1:1855326/93‑1 (MQ=255)

                                  ggCACAACAAGTGGCGCCCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCacg                                                       >  1:509412/1‑91 (MQ=255)

                                  ggCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAgg                                                       <  1:1232820/93‑1 (MQ=255)

                                   gCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGgc                                                      <  2:5081560/93‑1 (MQ=255)

                                   gCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGgc                                                      <  2:3600574/93‑1 (MQ=255)

                                    cacaACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGgcg                                                     <  1:988670/93‑1 (MQ=255)

                                    cacaACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGgcg                                                     <  1:1516477/93‑1 (MQ=255)

                                      caacaaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGcgcg                                                   <  1:4513812/93‑1 (MQ=255)

                                      caacaaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGcgcg                                                   >  2:4425207/1‑93 (MQ=255)

                                         caaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGAtt                                                >  2:1993136/1‑93 (MQ=255)

                                          aaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtc                                               >  1:1089824/1‑93 (MQ=255)

                                           aGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtct                                              >  2:488076/1‑93 (MQ=255)

                                            gTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtctc                                             <  1:5197526/93‑1 (MQ=255)

                                             tGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtctct                                            >  2:1255957/1‑93 (MQ=255)

                                             tGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtctct                                            <  2:5265207/93‑1 (MQ=255)

                                             tGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtctct                                            <  2:5213768/93‑1 (MQ=255)

                                              ggCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtctctc                                           <  2:2188401/93‑1 (MQ=255)

                                              ggCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATtctctc                                           <  2:3864502/93‑1 (MQ=255)

                                                cGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGc                                         <  2:5347897/93‑1 (MQ=255)

                                                  aCCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCtt                                       >  2:147566/1‑93 (MQ=255)

                                                   ccTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCttt                                      <  2:6077947/93‑1 (MQ=255)

                                                   ccTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCttt                                      >  1:2884134/1‑93 (MQ=255)

                                                   ccTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCttt                                      <  2:4165382/93‑1 (MQ=255)

                                                     ttACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTtct                                    >  1:2074357/1‑93 (MQ=255)

                                                     ttACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTtct                                    >  1:2240763/1‑93 (MQ=255)

                                                     ttACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTtct                                    >  1:5783013/1‑93 (MQ=255)

                                                       acacGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTctct                                  >  2:4607402/1‑93 (MQ=255)

                                                        cacGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTg                                 <  1:5885972/93‑1 (MQ=255)

                                                        cacGGCGGTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTg                                 <  2:3691941/93‑1 (MQ=255)

                                                          cggcggTAAATTCCCTGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTggg                               <  2:481247/93‑1 (MQ=255)

                                                               gTAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGc                          >  2:5583862/1‑93 (MQ=255)

                                                                tAAATTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGAGATTCTCTCGCTTTCTCTGGGTAAGCa                         <  2:5583013/93‑1 (MQ=255)

                                                                   aTTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACtt                      <  2:2984210/93‑1 (MQ=255)

                                                                   aTTCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACtt                      <  2:75616/93‑1 (MQ=255)

                                                                    ttCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACttt                     <  1:4857541/93‑1 (MQ=255)

                                                                    ttCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACttt                     <  1:1418323/93‑1 (MQ=255)

                                                                     tCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTc                    <  2:5700199/93‑1 (MQ=255)

                                                                     tCCCGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTc                    <  2:107360/93‑1 (MQ=255)

                                                                      cccGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTcc                   <  2:4353313/93‑1 (MQ=255)

                                                                      cccGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTcc                   <  2:2364225/93‑1 (MQ=255)

                                                                       ccGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGGGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCg                  >  2:3433389/1‑93 (MQ=255)

                                                                        cGGAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGa                 <  2:6379566/93‑1 (MQ=255)

                                                                         ggAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGAt                <  1:3552387/93‑1 (MQ=255)

                                                                          gAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGAtt               <  2:2705057/93‑1 (MQ=255)

                                                                          gAAACCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGAtt               <  2:5760894/93‑1 (MQ=255)

                                                                           aaaCCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtc              <  1:1007813/93‑1 (MQ=255)

                                                                            aaCCTTTGAGCAAGATGCAGCACTGCCGGGCGTCGGGCGTCCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtct             <  1:5748550/93‑1 (MQ=255)

                                                                             aCCTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtctc            <  2:2376784/93‑1 (MQ=255)

                                                                             aCCTTTGAAGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtctc            >  1:6308260/1‑93 (MQ=255)

                                                                              ccTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCg           <  2:165365/93‑1 (MQ=255)

                                                                              ccTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCg           >  2:1430272/1‑93 (MQ=255)

                                                                               cTTTGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGa          <  1:6332980/93‑1 (MQ=255)

                                                                                tttGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGAc         <  2:3219305/93‑1 (MQ=255)

                                                                                  tGAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACgg       <  1:1995738/93‑1 (MQ=255)

                                                                                   gAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGt      <  1:2444491/93‑1 (MQ=255)

                                                                                   gAGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGt      >  1:4687985/1‑93 (MQ=255)

                                                                                    aGGAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTa     <  2:2127549/93‑1 (MQ=255)

                                                                                    aGGAAGAGGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTa     <  2:5985810/93‑1 (MQ=255)

                                                                                     ggAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTaa    >  2:259358/1‑93 (MQ=255)

                                                                                      gAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAAc   <  2:251604/93‑1 (MQ=255)

                                                                                      gAAGATGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAAc   >  2:2242413/1‑93 (MQ=255)

GCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACG  >  NC_000913/3101272‑3101451

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