BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,469,173	T→C	D122G (GAC→GGC)	b2354 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,469,173	0	T	C	98.7%	495.9 / ‑2.6	155	D122G (GAC→GGC)	b2354	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (73/80); total (74/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 162 total aligned reads displayed.

CACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTC  >  NC_000913/2469083‑2469262

cACCGCGCCGGATGTCAACTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGccc                                                                                         <  2:795500/93‑1 (MQ=255)

 aCCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCa                                                                                        <  1:3645149/93‑1 (MQ=255)

 aCCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCa                                                                                        >  1:3429999/1‑93 (MQ=255)

     cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATa                                                                                    >  2:1935505/1‑93 (MQ=255)

     cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATa                                                                                    >  2:1850633/1‑93 (MQ=255)

       ccGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAAt                                                                                  <  1:1451491/93‑1 (MQ=255)

        cGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATc                                                                                 >  2:3504690/1‑93 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGATGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                <  2:1923940/93‑1 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATcc                                                                                >  2:4514455/1‑93 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                <  2:1064408/93‑1 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                >  1:1954582/1‑93 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                >  2:1131815/1‑93 (MQ=255)

         ggATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATcc                                                                                <  2:3138633/93‑1 (MQ=255)

          gATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCa                                                                               >  2:2962337/1‑93 (MQ=255)

          gATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCa                                                                               >  1:325740/1‑93 (MQ=255)

           aTGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAg                                                                              >  2:168634/1‑93 (MQ=255)

              tCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCca                                                                           <  1:2777642/93‑1 (MQ=255)

               cATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATCATCCAGCcac                                                                          >  1:1672382/1‑93 (MQ=255)

               cATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCcac                                                                          <  1:4573134/93‑1 (MQ=255)

                aTCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCcaca                                                                         <  1:3002013/93‑1 (MQ=255)

                 tCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACaa                                                                        >  2:4609804/1‑93 (MQ=255)

                  cTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACaaa                                                                       >  2:3401639/1‑93 (MQ=255)

                   tgtccTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAAt                                                                      <  1:2179721/90‑1 (MQ=255)

                   tGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAAt                                                                      >  2:3039280/1‑93 (MQ=255)

                    ggCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAAtt                                                                     >  2:2575523/1‑93 (MQ=255)

                    ggCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAAtt                                                                     >  1:2253125/1‑93 (MQ=255)

                     gCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATta                                                                    <  2:4190392/93‑1 (MQ=255)

                     gCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATta                                                                    <  1:887551/93‑1 (MQ=255)

                       cTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATtata                                                                  <  2:759009/93‑1 (MQ=255)

                       cTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATtata                                                                  >  1:788430/1‑93 (MQ=255)

                       cTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTAGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATtata                                                                  <  2:216673/93‑1 (MQ=255)

                       cTGAACTGCAGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGTCCAGCTCAAGTGCGCCCAGATAATCCAGCCACAAATtata                                                                  <  2:3517465/93‑1 (MQ=17)

                         taaCTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATAtt                                                                <  2:1489975/92‑1 (MQ=255)

                         gAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATAtt                                                                <  2:85986/93‑1 (MQ=255)

                           aCTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTgc                                                              >  2:787270/1‑93 (MQ=255)

                            cTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTgcg                                                             <  1:4318073/93‑1 (MQ=255)

                            cTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTgcg                                                             <  1:3968000/93‑1 (MQ=255)

                               ccGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAGGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCa                                                          <  2:3983422/93‑1 (MQ=255)

                                 ggaggaGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGt                                                        >  1:3175309/1‑93 (MQ=255)

                                 ggaggaGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGt                                                        >  2:4437642/1‑93 (MQ=255)

                                  gaggagTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGtt                                                       >  1:1877830/1‑93 (MQ=255)

                                   aggagTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTc                                                      <  1:1426013/93‑1 (MQ=255)

                                       gTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCa                                                  >  2:1182761/1‑93 (MQ=255)

                                         aGGCCATTCAATATCTGGCGCATTGAAAGTATCGACCAGTTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTAAcc                                                <  2:2831966/93‑1 (MQ=17)

                                         aGGCCATGCAAGATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCAcc                                                <  2:3294835/93‑1 (MQ=255)

                                          ggCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCt                                               <  2:2469121/93‑1 (MQ=255)

                                           gCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCtt                                              >  1:628954/1‑93 (MQ=255)

                                           gCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCtt                                              >  1:2966576/1‑93 (MQ=255)

                                           gCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCtt                                              >  2:3574501/1‑93 (MQ=255)

                                           gCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCtt                                              <  1:1997359/93‑1 (MQ=255)

                                           gCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCtt                                              >  2:3716350/1‑93 (MQ=255)

                                              aTTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCa                                           <  1:3293950/93‑1 (MQ=255)

                                              aTTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCa                                           >  2:472652/1‑93 (MQ=255)

                                               ttCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAg                                          <  1:4732902/93‑1 (MQ=255)

                                               ttCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAg                                          >  2:4603108/1‑93 (MQ=255)

                                                 cAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgac                                        <  1:3923199/93‑1 (MQ=255)

                                                  aaTATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacg                                       >  1:1706410/1‑93 (MQ=255)

                                                   atatCGGGCGCCCTGGAAGTATCGACCAGCTCCAGTGGGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                      <  1:2150144/93‑1 (MQ=25)

                                                    tatCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacgac                                     >  1:3706979/1‑93 (MQ=255)

                                                    tatCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacgac                                     >  2:3002831/1‑93 (MQ=255)

                                                       cTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCaa                                  >  2:2238579/1‑93 (MQ=255)

                                                       cTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCaa                                  <  2:3402411/93‑1 (MQ=255)

                                                        tGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAAt                                 <  2:3692191/93‑1 (MQ=255)

                                                          gcgcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAg                               <  2:2529635/93‑1 (MQ=255)

                                                           gcaaCTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGc                              <  2:3665230/90‑1 (MQ=255)

                                                           cgcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGc                              <  2:3083450/93‑1 (MQ=255)

                                                            gcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGAGCAATAGcg                             >  2:2956480/1‑92 (MQ=255)

                                                            gcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGcc                             >  2:487157/1‑93 (MQ=255)

                                                            gcACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGcc                             <  2:4067380/93‑1 (MQ=255)

                                                             cACTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCg                            <  1:2766621/93‑1 (MQ=255)

                                                              aCTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGc                           <  1:3425444/93‑1 (MQ=255)

                                                              aCTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGc                           <  2:4005198/93‑1 (MQ=255)

                                                               cTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCt                          <  1:2147124/93‑1 (MQ=255)

                                                               cTGGAAGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCt                          <  1:3900046/93‑1 (MQ=255)

                                                                  gaggTATCGACCAGCTCCATGGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTa                       <  2:1588847/90‑1 (MQ=38)

                                                                    aGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAcc                     <  1:3700976/93‑1 (MQ=255)

                                                                    aGTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTAcc                     >  2:3187488/1‑93 (MQ=255)

                                                                    aGTATCGACCACCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTTCGCCAGGTCCTCACCTTTCAGACGAGCAATAGCCGCTTTatc                     >  2:4160802/1‑91 (MQ=37)

                                                                     gTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCa                    <  1:4088427/93‑1 (MQ=255)

                                                                     gTATCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCa                    >  2:143350/1‑93 (MQ=255)

                                                                       aTCGACCAGCTCCAGCGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc                  >  1:948894/1‑93 (MQ=255)

                                                                        tCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGcc                 >  1:3833538/1‑93 (MQ=255)

                                                                        tCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGcc                 >  1:3124575/1‑93 (MQ=255)

                                                                        tCGACCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGcc                 <  1:3606425/93‑1 (MQ=255)

                                                                           aCCAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCAt              <  1:4182710/93‑1 (MQ=255)

                                                                            ccAGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCAtt             <  1:3323121/93‑1 (MQ=255)

                                                                              aGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGt           <  1:2239359/93‑1 (MQ=255)

                                                                              aGCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGt           <  1:2461024/93‑1 (MQ=255)

                                                                               gCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGtt          >  2:4614101/1‑93 (MQ=255)

                                                                               gCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGtt          <  1:3176600/93‑1 (MQ=255)

                                                                               gCTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGtt          >  2:634640/1‑93 (MQ=255)

                                                                                cTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGttt         <  2:1251205/93‑1 (MQ=255)

                                                                                cTCCAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGttt         >  1:4113748/1‑93 (MQ=255)

                                                                                 cccAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa        <  2:1667779/92‑1 (MQ=255)

                                                                                  ccAGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTAc       <  2:4400993/93‑1 (MQ=255)

                                                                                    aGTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTg     >  1:3484015/1‑93 (MQ=255)

                                                                                     gTGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGt    >  2:1025257/1‑93 (MQ=255)

                                                                                      tGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGtt   <  2:65782/93‑1 (MQ=255)

                                                                                      tGCGCCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGtt   >  2:2744602/1‑93 (MQ=255)

                                                                                       gcgcCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTc  <  2:2531500/93‑1 (MQ=255)

CACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACTGTTC  >  NC_000913/2469083‑2469262

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