BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	500,396	(A)_6→5	coding (1048/1305 nt)	gsk →	inosine‑guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	500,391	0	A	.	99.6%	1135.6 / ‑3.3	242	coding (1043/1305 nt)	gsk	inosine‑guanosine kinase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (124/117); minor base G (0/1); total (124/118)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.09e-01

Only 100 of 265 total aligned reads displayed.

AGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAAC  >  NC_000913/500301‑500479

aGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaaa                                                                                        <  1:4639029/93‑1 (MQ=255)

 gTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaaaa                                                                                       >  1:3392708/1‑93 (MQ=255)

  tATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaaaaa                                                                                      >  2:2405791/1‑93 (MQ=255)

   aTGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaaaaat                                                                                     <  1:285309/93‑2 (MQ=255)

   aTGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaaaaat                                                                                     >  1:1171584/1‑92 (MQ=255)

    tGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGaaaaatc                                                                                    <  1:3852141/93‑3 (MQ=255)

       gTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATg                                                                                >  2:3195461/1‑93 (MQ=255)

        tttAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGa                                                                               >  1:6018348/1‑93 (MQ=255)

         ttAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGaa                                                                              <  2:4893508/93‑1 (MQ=255)

          tAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAac                                                                             >  2:4333102/1‑93 (MQ=255)

           aGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAaca                                                                            <  2:2800771/93‑1 (MQ=255)

            gCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAacac                                                                           >  1:1896864/1‑93 (MQ=255)

            gCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAacac                                                                           <  1:6193277/93‑1 (MQ=255)

            gCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAacac                                                                           >  1:705166/1‑93 (MQ=255)

            gCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAacac                                                                           <  1:3243537/93‑1 (MQ=255)

                cgcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAAt                                                                       <  1:5707052/93‑1 (MQ=255)

                cgcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAAt                                                                       <  2:5985534/93‑1 (MQ=255)

                cgcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAAt                                                                       <  1:1582640/93‑1 (MQ=255)

                cgcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAAt                                                                       <  1:783751/93‑1 (MQ=255)

                 gcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATg                                                                      <  1:4505736/93‑1 (MQ=255)

                 gcCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATACGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATg                                                                      >  2:4530113/1‑93 (MQ=255)

                    aTGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAg                                                                   <  1:2723498/93‑1 (MQ=255)

                       cgcCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCgg                                                                >  2:2092538/1‑93 (MQ=255)

                        gcCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCggg                                                               <  2:5963484/93‑1 (MQ=255)

                        gcCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCggg                                                               >  2:6271606/1‑93 (MQ=255)

                         ccACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCgggg                                                              >  2:175320/1‑93 (MQ=255)

                           acaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGa                                                            >  1:936705/1‑93 (MQ=255)

                           acaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGa                                                            >  2:2160206/1‑93 (MQ=255)

                           acaAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCACGCACACTAATGGAGCggagga                                                            >  1:915252/1‑93 (MQ=255)

                            caAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGag                                                           >  2:6467511/1‑92 (MQ=255)

                            caAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGAt                                                           <  2:5675970/93‑1 (MQ=255)

                              aGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATgg                                                         >  1:1309005/1‑93 (MQ=255)

                              aGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATgg                                                         >  2:1958312/1‑93 (MQ=255)

                               ggATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGgc                                                        >  1:4185387/1‑93 (MQ=255)

                                  ttGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCa                                                     <  2:2573736/93‑1 (MQ=255)

                                  ttGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCa                                                     <  2:640313/93‑1 (MQ=255)

                                    gCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCAtt                                                   <  1:3550837/93‑1 (MQ=255)

                                      cAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTgg                                                 >  1:734570/1‑93 (MQ=255)

                                      cAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTgg                                                 >  1:2204701/1‑93 (MQ=255)

                                      cAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTgg                                                 <  2:1882695/93‑1 (MQ=255)

                                       aGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGc                                                <  2:5231403/93‑1 (MQ=255)

                                         aaTCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAg                                              <  1:5425901/93‑1 (MQ=255)

                                            ccGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGt                                           >  2:6198592/1‑93 (MQ=255)

                                            ccGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGt                                           <  1:195570/93‑1 (MQ=255)

                                             cGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGtt                                          >  1:1033419/1‑93 (MQ=255)

                                               cTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgc                                        <  2:6340069/93‑1 (MQ=255)

                                               cTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgc                                        <  2:4188555/93‑1 (MQ=255)

                                               cTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgc                                        <  2:2956034/93‑1 (MQ=255)

                                               cTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgc                                        <  2:2313478/93‑1 (MQ=255)

                                               cTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgc                                        <  2:1072225/93‑1 (MQ=255)

                                               cTGCGTGGATAATAGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGAtgc                                        >  1:4409358/1‑93 (MQ=255)

                                                 gCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTtgctg                                      >  2:5535783/1‑93 (MQ=255)

                                                   gtgtATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCa                                    <  1:5953180/93‑1 (MQ=255)

                                                    tgtATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCAt                                   <  2:3879071/93‑1 (MQ=255)

                                                    tgtATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCAt                                   <  2:5486936/93‑1 (MQ=255)

                                                     gtATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATg                                  >  2:5726400/1‑93 (MQ=255)

                                                     gtATATTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATg                                  >  1:5066917/1‑93 (MQ=255)

                                                      tataTTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGa                                 >  1:3823307/1‑93 (MQ=255)

                                                        tatTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACa                               <  2:3645755/93‑1 (MQ=255)

                                                        tatTCGCACATTGCGCAGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACa                               <  2:5162744/93‑1 (MQ=255)

                                                         atTCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACAt                              <  1:3484628/93‑1 (MQ=255)

                                                           tCGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTa                            <  2:4106170/93‑1 (MQ=255)

                                                            cGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTAc                           <  2:1344376/93‑1 (MQ=255)

                                                            cGCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTAc                           >  1:5598440/1‑93 (MQ=255)

                                                             gCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTAcc                          <  1:6373219/93‑1 (MQ=255)

                                                             gCACATTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTAcc                          >  1:3811611/1‑93 (MQ=255)

                                                              cacaTTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCg                         <  2:782878/93‑1 (MQ=255)

                                                              cacaTTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCg                         <  1:4783192/93‑1 (MQ=255)

                                                              caaaTTGCGCCGCACATGGGCGGGCCGG‑AAAAATCATGAACACCAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCg                         <  2:4675173/90‑1 (MQ=255)

                                                               acaTTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGc                        >  2:286898/1‑93 (MQ=255)

                                                                 aTTGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCa                      >  1:1783412/1‑93 (MQ=255)

                                                                  ttGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCaa                     <  1:2679531/93‑1 (MQ=255)

                                                                   tGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAAc                    <  1:1058856/93‑1 (MQ=255)

                                                                   tGCGCCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAAc                    <  1:5267160/93‑1 (MQ=255)

                                                                    gcgcCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACa                   <  2:1588857/93‑1 (MQ=255)

                                                                    gcgcCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACa                   >  1:1557049/1‑93 (MQ=255)

                                                                     cgcCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAg                  <  1:1390341/93‑1 (MQ=255)

                                                                      gcCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGc                 >  2:1633930/1‑93 (MQ=255)

                                                                      gcCGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGc                 <  1:615168/93‑1 (MQ=255)

                                                                       ccGTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCt                <  1:4150510/93‑1 (MQ=255)

                                                                       ccGTACATGGGCGGCCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCt                <  2:3635564/93‑1 (MQ=255)

                                                                         gTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTAc              >  2:2212645/1‑93 (MQ=255)

                                                                         gTACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTAc              >  1:3119675/1‑93 (MQ=255)

                                                                          tACATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTAcc             >  2:3224692/1‑93 (MQ=255)

                                                                           aCATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCa            <  1:3095999/93‑1 (MQ=255)

                                                                           aCATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCa            <  2:3465845/93‑1 (MQ=255)

                                                                           aCATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCa            <  1:1936451/93‑1 (MQ=255)

                                                                            cATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCAt           <  1:5860355/93‑1 (MQ=255)

                                                                            cATGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCAt           <  2:1145774/93‑1 (MQ=255)

                                                                             aTGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATc          <  1:3338882/93‑1 (MQ=255)

                                                                             aTGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATc          >  1:3716181/1‑93 (MQ=255)

                                                                             aTGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATc          <  1:2328110/93‑1 (MQ=255)

                                                                             aTGGGCGGGCCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATc          <  1:4319481/93‑1 (MQ=255)

                                                                               gggcgggcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGt        >  2:5775297/1‑93 (MQ=255)

                                                                               gggcgggcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGt        >  1:2433110/1‑93 (MQ=255)

                                                                                ggcgggcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTa       >  1:4854398/1‑93 (MQ=255)

                                                                                  cgggcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGc     <  1:5185491/93‑1 (MQ=255)

                                                                                    ggcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCaa   <  2:6136607/93‑1 (MQ=255)

                                                                                    ggcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCaa   >  2:1474513/1‑93 (MQ=255)

                                                                                     gcCGG‑AAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAAc  >  2:5275788/1‑93 (MQ=255)

AGTATGAGTTTAGCCGCGCCATGCGCCACAAGGATTGCCAGAATCCGCTGCGTGTATATTCGCACATTGCGCCGTACATGGGCGGGCCGGAAAAAATCATGAACACTAATGGAGCGGGGGATGGCGCATTGGCAGCGTTGCTGCATGACATTACCGCCAACAGCTACCATCGTAGCAAC  >  NC_000913/500301‑500479

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