BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	866,644	A→G	N285S (AAT→AGT)	ybiK →	putative asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	866,644	0	A	G	99.5%	628.4 / NA	199	N285S (AAT→AGT)	ybiK	putative asparaginase
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (92/106); total (93/106)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 207 total aligned reads displayed.

GTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAG  >  NC_000913/866554‑866733

gTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTg                                                                                         <  2:1949423/93‑1 (MQ=255)

gTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTg                                                                                         <  2:4907330/93‑1 (MQ=255)

  ctcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTc                                                                                       <  2:4490330/93‑1 (MQ=255)

   tcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcg                                                                                      <  2:3984458/93‑1 (MQ=255)

   tcGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcg                                                                                      <  1:2387897/93‑1 (MQ=255)

    cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgc                                                                                     <  2:6168099/93‑1 (MQ=255)

    cgcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgc                                                                                     >  2:3401896/1‑93 (MQ=255)

     gcgGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTcgcg                                                                                    >  2:146691/1‑93 (MQ=255)

       ggAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCt                                                                                  >  1:3014010/1‑93 (MQ=255)

       ggAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCt                                                                                  <  2:3727307/93‑1 (MQ=255)

        gAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTa                                                                                 >  1:1447366/1‑93 (MQ=255)

         aaGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTAc                                                                                <  1:5641987/93‑1 (MQ=255)

         aaGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTAc                                                                                <  2:3389320/93‑1 (MQ=255)

          aGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTAcc                                                                               >  2:1318564/1‑93 (MQ=255)

           gCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCg                                                                              >  2:5301951/1‑93 (MQ=255)

             cTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGtt                                                                            >  2:2316756/1‑93 (MQ=255)

             cTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGtt                                                                            >  1:3078997/1‑93 (MQ=255)

              tGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGttt                                                                           <  1:599457/93‑1 (MQ=255)

                cGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTaa                                                                         <  1:1649501/93‑1 (MQ=255)

                  aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       <  1:1927495/93‑1 (MQ=255)

                  aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       >  2:5360757/1‑93 (MQ=255)

                  aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAaca                                                                       >  1:1279801/1‑93 (MQ=255)

                   gCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAacac                                                                      >  1:744262/1‑93 (MQ=255)

                     gggTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCg                                                                    >  2:2653223/1‑93 (MQ=255)

                     gggCAGTAATGGAAAAAATCCCTGCGCTTGGCGGTAGCGGTGGATTAATAGCTATCGACCATGAAGGGAGAGTAGCGCTACCGTTTAACACCg                                                                    >  2:5612001/1‑93 (MQ=255)

                      ggTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGa                                                                   <  1:1548498/93‑1 (MQ=255)

                       gtagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  >  2:5922623/1‑93 (MQ=255)

                       gtagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  <  2:2049558/93‑1 (MQ=255)

                       gtagtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  >  2:1781004/1‑93 (MQ=255)

                       gtagtaATGGAAAAACTCCCTGCGCTCGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGaa                                                                  <  1:1494503/93‑1 (MQ=255)

                         agtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGTCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAgg                                                                <  2:3564749/93‑1 (MQ=255)

                         agtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAgg                                                                >  2:1181795/1‑93 (MQ=255)

                          gtaATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGa                                                               >  2:2109485/1‑93 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              <  2:2813360/93‑1 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              <  2:3872335/93‑1 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              <  1:4629482/93‑1 (MQ=255)

                           taATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGaa                                                              >  2:1673344/1‑93 (MQ=255)

                            aaTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGTAGGAAt                                                             >  1:3468139/1‑93 (MQ=255)

                            aaTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAAt                                                             <  1:3955399/93‑1 (MQ=255)

                            aaTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAAt                                                             >  2:100126/1‑93 (MQ=255)

                             aTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATg                                                            <  2:1934495/93‑1 (MQ=255)

                             aTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATg                                                            >  2:4835741/1‑93 (MQ=255)

                             aTGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATg                                                            >  1:1128697/1‑93 (MQ=255)

                              tGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGt                                                           >  1:1970098/1‑93 (MQ=255)

                               ggAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTaa                                                                         >  1:634059/1‑78 (MQ=255)

                               ggAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTaa                                                                         <  2:634059/78‑1 (MQ=255)

                                gAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTAt                                                         >  2:5402499/1‑93 (MQ=255)

                                 aaaaaCTCCCTGCGCTTGGCGGTAGCGGTGGGTTAATCGCTATCGACCATGAAGGGAGTGGCGCGCTACCGTTTAACACCGAAGGAATGTATc                                                        <  1:3914023/93‑1 (MQ=255)

                                 aaaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATc                                                        <  2:4118013/93‑1 (MQ=255)

                                 aaaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATc                                                        >  1:5019027/1‑93 (MQ=255)

                                  aaaaCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATcg                                                       <  1:3276122/93‑1 (MQ=255)

                                     aCTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCgcgc                                                    >  1:5222387/1‑93 (MQ=255)

                                      cTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGcc                                                   <  1:2815989/93‑1 (MQ=255)

                                      cTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGcc                                                   <  2:3242457/93‑1 (MQ=255)

                                       tCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCt                                                  <  2:2726333/93‑1 (MQ=255)

                                       tCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCt                                                  <  2:55420/93‑1 (MQ=255)

                                       tCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCt                                                  >  2:1732953/1‑93 (MQ=255)

                                         ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                <  2:6248287/93‑1 (MQ=255)

                                         ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                >  1:2610747/1‑93 (MQ=255)

                                         ccTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTgg                                                <  2:3527779/93‑1 (MQ=255)

                                          cTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTggg                                               <  1:1334634/93‑1 (MQ=255)

                                            gcgcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGc                                             <  1:4651984/93‑1 (MQ=255)

                                            gcgcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGc                                             >  1:3668146/1‑93 (MQ=255)

                                              gcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTa                                           <  1:1081423/93‑1 (MQ=255)

                                               cTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTAc                                          <  1:4387279/93‑1 (MQ=255)

                                                ttGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACg                                         >  2:2336414/1‑93 (MQ=255)

                                                ttGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACg                                         <  1:3535038/93‑1 (MQ=255)

                                                  ggCGGTAGCGGTGGCTTAATCGCTATCGACCATGAACGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCa                                       <  2:1504084/93‑1 (MQ=255)

                                                   gCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAg                                      <  1:2517172/93‑1 (MQ=255)

                                                    cGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAgg                                     <  2:708485/93‑1 (MQ=255)

                                                     ggTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGc                                    <  2:744920/93‑1 (MQ=255)

                                                       tAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGa                                  <  1:2051483/93‑1 (MQ=255)

                                                         gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  1:4139192/93‑1 (MQ=255)

                                                         gCGGTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATa                                <  1:2000475/93‑1 (MQ=255)

                                                           ggTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACg                              <  1:5859528/93‑1 (MQ=255)

                                                           ggTGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACg                              >  1:4560484/1‑93 (MQ=255)

                                                             tGGCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGcc                            >  2:5676654/1‑93 (MQ=255)

                                                               gCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCaa                          >  2:1555599/1‑93 (MQ=255)

                                                               gCTTAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCaa                          >  2:1099346/1‑93 (MQ=255)

                                                                 ttAATCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAacc                        >  2:5062942/1‑93 (MQ=255)

                                                                    aTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCCGGCGATACGCCCaccacc                     >  2:3081914/1‑93 (MQ=255)

                                                                    aTCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAaccacc                     <  1:4409482/93‑1 (MQ=255)

                                                                     tCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCg                    <  1:876598/93‑1 (MQ=255)

                                                                     tCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCg                    >  2:4739029/1‑93 (MQ=255)

                                                                     tCGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCg                    >  2:465046/1‑93 (MQ=255)

                                                                      cGCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCgg                   <  1:3770871/93‑1 (MQ=255)

                                                                       gCTATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGt                  <  2:505096/93‑1 (MQ=255)

                                                                         tATCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTAt                >  1:2183336/1‑93 (MQ=255)

                                                                          aTCGACCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATc               <  1:1616900/93‑1 (MQ=255)

                                                                              aCCATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTAcc           <  2:580525/93‑1 (MQ=255)

                                                                                cATGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGCCTACGCAGGCGATACGCCAACCACCGGTATCTACCGt         <  1:6235376/93‑1 (MQ=255)

                                                                                  tGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGa       >  2:4253611/1‑93 (MQ=255)

                                                                                  tGAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGa       >  1:5336976/1‑93 (MQ=255)

                                                                                   gAAGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaa      >  2:2967622/1‑93 (MQ=255)

                                                                                     aGGGAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaa    >  1:40853/1‑93 (MQ=255)

                                                                                      gggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaaa   >  2:4163609/1‑93 (MQ=255)

                                                                                      gggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGaaaaa   <  1:3513662/93‑1 (MQ=255)

                                                                                       ggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAg  <  2:5416374/93‑1 (MQ=255)

                                                                                       ggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAg  <  2:4471158/93‑1 (MQ=255)

                                                                                       ggAGTGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAg  >  1:912780/1‑93 (MQ=255)

GTCTCGCGGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAG  >  NC_000913/866554‑866733

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