BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,628,759	G→A	R488C (CGT→TGT)	yhiJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,628,759	0	G	A	99.8%	2945.7 / ‑5.0	790	R488C (CGT→TGT)	yhiJ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/1); new base A (426/362); total (427/363)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

Only 100 of 820 total aligned reads displayed.

AGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACGGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTAACA  >  NC_000913/3628667‑3628848

  aCTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACagg                                                                                         <  1:6983539/93‑1 (MQ=255)

  aCTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACagg                                                                                         <  2:3840502/93‑1 (MQ=255)

  aCTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACagg                                                                                         <  1:3644505/93‑1 (MQ=255)

   cTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGc                                                                                        <  2:1473153/93‑1 (MQ=255)

    ttCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCg                                                                                       >  1:1837683/1‑93 (MQ=255)

    ttCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCg                                                                                       >  1:5096661/1‑93 (MQ=255)

    ttCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCg                                                                                       <  2:6764839/93‑1 (MQ=255)

     tCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGt                                                                                      >  2:1407727/1‑93 (MQ=255)

       cAGGGAGTGGTTATCTCGGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTa                                                                                    <  2:5036048/93‑1 (MQ=255)

       cAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTa                                                                                    <  1:7617224/93‑1 (MQ=255)

       cAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTa                                                                                    >  2:8675841/1‑93 (MQ=255)

        aGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTaa                                                                                   <  2:9627662/93‑1 (MQ=255)

             gTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCt                                                                              >  1:205845/1‑93 (MQ=255)

              tGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTa                                                                             >  2:9045940/1‑93 (MQ=255)

              tGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTa                                                                             >  1:6103319/1‑93 (MQ=255)

              tGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTa                                                                             >  1:2517330/1‑93 (MQ=255)

               ggTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACGGGCGGTAAGACCTaa                                                                            >  2:6977608/1‑93 (MQ=255)

                 ttATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTaaaa                                                                          <  1:997578/93‑1 (MQ=255)

                  tATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAAc                                                                         <  2:9501450/93‑1 (MQ=255)

                    tctcAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTg                                                                       >  2:3653803/1‑93 (MQ=255)

                      tcAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGATCTAAAACTGTc                                                                     >  1:7266138/1‑93 (MQ=255)

                       cAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTcc                                                                    <  2:805281/93‑1 (MQ=255)

                         ggCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAGCATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTcccc                                                                  >  2:33695/1‑93 (MQ=255)

                          gCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCg                                                                 >  1:676067/1‑93 (MQ=255)

                           gccTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATCCTTCCAGCAGGTTTCCTGTCGCAGGCGTTAATACCTATAACTGGCCCCGt                                                                <  2:8779702/92‑1 (MQ=255)

                           cccTGCTCCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGt                                                                <  1:10026105/93‑1 (MQ=255)

                            ccTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGtt                                                               <  1:7151064/93‑1 (MQ=255)

                              cgCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTc                                                             <  2:2772671/92‑1 (MQ=255)

                                  ccAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCt                                                         >  1:9556763/1‑93 (MQ=255)

                                  ccAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCt                                                         <  1:8885203/93‑1 (MQ=255)

                                   cAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTa                                                        >  2:3961762/1‑93 (MQ=255)

                                   cAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTa                                                        >  2:6096339/1‑93 (MQ=255)

                                   cAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTa                                                        <  2:8543216/93‑1 (MQ=255)

                                    aGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTAt                                                       >  2:9314884/1‑93 (MQ=255)

                                    aGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTAt                                                       <  1:4877625/93‑1 (MQ=255)

                                     gggTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATg                                                      >  1:2737541/1‑93 (MQ=255)

                                      ggTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGc                                                     <  2:4834117/93‑1 (MQ=255)

                                      ggTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGc                                                     >  1:3530027/1‑93 (MQ=255)

                                       gTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCg                                                    <  1:6337775/93‑1 (MQ=255)

                                       gTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCg                                                    >  2:2041042/1‑93 (MQ=255)

                                         tAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGtt                                                  >  1:6227045/1‑93 (MQ=255)

                                         tAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGtt                                                  <  1:4341191/93‑1 (MQ=255)

                                          aaTTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTtc                                                 >  2:2456320/1‑93 (MQ=255)

                                            tttGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTtctc                                               >  1:8097715/1‑93 (MQ=255)

                                            tttGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTtctc                                               <  2:7795979/93‑1 (MQ=255)

                                            tttGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTtctc                                               <  1:2359080/93‑1 (MQ=255)

                                              tagCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCtt                                             <  1:5521978/91‑1 (MQ=255)

                                              tGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCtt                                             >  1:1193745/1‑93 (MQ=255)

                                              tGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCtt                                             <  1:921031/93‑1 (MQ=255)

                                              tGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACGGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCtt                                             <  1:6888450/93‑1 (MQ=255)

                                                gCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTcc                                           <  1:8510228/93‑1 (MQ=255)

                                                    ttATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTccacca                                       >  1:2232206/1‑93 (MQ=255)

                                                      atgcatgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCAtt                                     <  1:9368991/93‑1 (MQ=255)

                                                       tgcatgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCAttt                                    >  2:8275930/1‑93 (MQ=255)

                                                        tcatgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTCTCCCCGTTGCCCCTATGCGTTCTCTTCCTCCATTTc                                   <  1:5649834/92‑1 (MQ=255)

                                                        gcatgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTc                                   <  1:7747846/93‑1 (MQ=255)

                                                        gcatgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTc                                   <  2:185134/93‑1 (MQ=255)

                                                          atgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCg                                 <  2:5899675/93‑1 (MQ=255)

                                                           tgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGt                                >  1:6506626/1‑93 (MQ=255)

                                                           tgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGt                                >  1:583006/1‑93 (MQ=255)

                                                           tgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGt                                >  1:5305816/1‑93 (MQ=255)

                                                           tgcTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGt                                >  1:2925703/1‑93 (MQ=255)

                                                             cTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTa                              <  1:9034344/93‑1 (MQ=255)

                                                             cTGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTa                              <  2:3495338/93‑1 (MQ=255)

                                                              tGAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTAt                             >  1:1491498/1‑93 (MQ=255)

                                                               gAACATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTAtg                            >  1:9819518/1‑93 (MQ=255)

                                                                 aCATTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTAtgtg                          <  1:3900711/93‑1 (MQ=255)

                                                                   aTTCTTCCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGca                        <  2:1120436/93‑1 (MQ=255)

                                                                    ttcttcCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGcaa                       >  2:9816096/1‑93 (MQ=255)

                                                                    ttcttcCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGcaa                       <  2:2611165/93‑1 (MQ=255)

                                                                    ttcttcCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGcaa                       >  1:4115510/1‑93 (MQ=255)

                                                                    ttcttcCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGcaa                       >  2:568132/1‑93 (MQ=255)

                                                                       ttcCAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGcaacaa                    <  2:3088156/93‑1 (MQ=255)

                                                                         ccAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACt                  >  1:379443/1‑93 (MQ=255)

                                                                          cAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACtt                 <  2:4717385/93‑1 (MQ=255)

                                                                          cAGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACtt                 <  2:6346062/93‑1 (MQ=255)

                                                                           aGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTg                >  2:4153949/1‑93 (MQ=255)

                                                                           aGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTg                >  2:7233547/1‑93 (MQ=255)

                                                                           aGCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTg                >  1:7948803/1‑93 (MQ=255)

                                                                            gCATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGa               <  2:1691170/93‑1 (MQ=255)

                                                                             cATGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGAt              >  1:7430708/1‑93 (MQ=255)

                                                                              aTGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGAtt             >  2:2248370/1‑93 (MQ=255)

                                                                              aTGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGAtt             >  1:4546514/1‑93 (MQ=255)

                                                                              aTGTTTACTGTCACAGGCGTTAAGACCCAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGAtt             >  2:3267593/1‑93 (MQ=255)

                                                                               tGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGAttt            <  1:8355379/93‑1 (MQ=255)

                                                                               tGTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGAttt            <  2:661222/93‑1 (MQ=255)

                                                                                gtgtACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTta           <  2:5425355/90‑1 (MQ=255)

                                                                                gTTTACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTta           >  2:5277684/1‑93 (MQ=255)

                                                                                 tttACTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTtat          >  2:2153085/1‑93 (MQ=255)

                                                                                     cTGTCACAGGCGTTAAGTCCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCt      <  1:1678721/93‑1 (MQ=255)

                                                                                     cTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCt      >  2:7836079/1‑93 (MQ=255)

                                                                                     cTGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCt      >  1:6005115/1‑93 (MQ=255)

                                                                                      tGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTa     >  2:5196410/1‑93 (MQ=255)

                                                                                      tGTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTa     >  1:3331673/1‑93 (MQ=255)

                                                                                       gTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTaa    <  2:7591115/93‑1 (MQ=255)

                                                                                       gTCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTaa    <  1:8778121/93‑1 (MQ=255)

                                                                                        tCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTAAc   >  2:7391354/1‑93 (MQ=255)

                                                                                        tCACAGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTAAc   <  1:1390614/93‑1 (MQ=255)

                                                                                         cacaGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTAACa  >  1:6592526/1‑93 (MQ=255)

                                                                                         cacaGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTAACa  <  2:4551348/93‑1 (MQ=255)

AGACTTCCAGGGAGTGGTTATCTCAGGCCCTGCACCAGGGTTAATTTGGCAATTATGCATGCTGAACATTCTTCCAGCATGTTTACTGTCACGGGCGTTAAGACCTAAAACTGTCCCCGTTTCGCCTATGCGTTCTCTTCCACCATTTCCGTTATGTGCAACAACTTGATTTATACCTAACA  >  NC_000913/3628667‑3628848

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