BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	500,613	A→G	E422G (GAG→GGG)	gsk →	inosine‑guanosine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	500,613	0	A	G	97.9%	448.3 / ‑1.3	144	E422G (GAG→GGG)	gsk	inosine‑guanosine kinase
Reads supporting (aligned to +/- strand): ref base A (3/0); new base G (61/80); total (64/80)
Fisher's exact test for biased strand distribution p-value = 8.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 147 total aligned reads displayed.

CATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTT  >  NC_000913/500523‑500705

cATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGc                                                                                            <  1:1086937/93‑1 (MQ=255)

cATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGc                                                                                            >  1:4671201/1‑93 (MQ=255)

  tCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGt                                                                                          >  1:3747716/1‑93 (MQ=255)

  tCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGt                                                                                          <  2:3061901/93‑1 (MQ=255)

  tCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGt                                                                                          >  1:2894708/1‑93 (MQ=255)

   cGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTg                                                                                         >  2:2081566/1‑93 (MQ=255)

    gTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGa                                                                                        >  2:1280886/1‑93 (MQ=255)

    gTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGa                                                                                        >  1:4603196/1‑93 (MQ=255)

    gTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGa                                                                                        >  2:3927517/1‑93 (MQ=255)

       aGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGa                                                                                     <  1:3400042/93‑1 (MQ=255)

         cgcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTGCTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACa                                                                                   <  1:2317677/93‑1 (MQ=255)

          gcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAg                                                                                  <  2:1842866/93‑1 (MQ=255)

          gcAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAg                                                                                  >  1:1031904/1‑93 (MQ=255)

            aGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGcc                                                                                >  2:3914085/1‑93 (MQ=255)

            aGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGcc                                                                                <  1:3756137/93‑1 (MQ=255)

            aGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGcc                                                                                <  1:1419387/93‑1 (MQ=255)

            aGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGcc                                                                                >  1:203651/1‑93 (MQ=255)

             ggTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCt                                                                               >  1:2528407/1‑93 (MQ=255)

             ggTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCt                                                                               <  1:394456/93‑1 (MQ=255)

              gtgtgtAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTg                                                                              <  1:1741344/93‑1 (MQ=255)

                 tgtAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGaa                                                                           <  2:2143822/93‑1 (MQ=255)

                 tgtAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGaa                                                                           >  2:55387/1‑93 (MQ=255)

                 tgtAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGaa                                                                           <  2:1669576/93‑1 (MQ=255)

                 tgtAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGaa                                                                           >  2:735333/1‑93 (MQ=255)

                   tAAATATGCTAACCGTGTGGGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAaga                                                                         <  2:2612457/93‑1 (MQ=255)

                   tAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAaga                                                                         <  1:2728432/93‑1 (MQ=255)

                    aaaTATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAagag                                                                        <  2:1777356/93‑1 (MQ=255)

                    aaaTATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAagag                                                                        >  1:4675925/1‑93 (MQ=255)

                      atatGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTc                                                                      <  1:880928/93‑1 (MQ=255)

                       tatGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCt                                                                     >  2:423401/1‑93 (MQ=255)

                        atGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCtt                                                                    <  1:2646890/93‑1 (MQ=255)

                        atGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCtt                                                                    <  1:2353045/93‑1 (MQ=255)

                          gCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTAc                                                                  >  1:316591/1‑93 (MQ=255)

                           cTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACt                                                                 <  1:3666505/93‑1 (MQ=255)

                           cTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACt                                                                 >  1:1651423/1‑93 (MQ=255)

                            tAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTg                                                                >  2:524258/1‑93 (MQ=255)

                             aaCCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTgg                                                               <  1:4621810/93‑1 (MQ=255)

                             aaCCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTgg                                                               >  2:3133648/1‑93 (MQ=255)

                             aaCCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTgg                                                               <  1:4016944/93‑1 (MQ=255)

                                 gtgtGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATc                                                           <  1:2476200/93‑1 (MQ=255)

                                  tgtgAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCg                                                          <  1:3308760/93‑1 (MQ=255)

                                   gtgAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGt                                                         >  2:4231181/1‑93 (MQ=255)

                                   gtgAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGt                                                         <  2:4720868/93‑1 (MQ=255)

                                    tgAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGtt                                                        <  2:1574686/93‑1 (MQ=255)

                                     gAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTa                                                       <  2:1907090/93‑1 (MQ=255)

                                     gAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTa                                                       >  1:537590/1‑93 (MQ=255)

                                      aGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTaa                                                      <  2:3442616/93‑1 (MQ=255)

                                       gCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAg                                                     >  1:3823065/1‑93 (MQ=255)

                                       gCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAg                                                     <  1:1131030/93‑1 (MQ=255)

                                       gCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAg                                                     >  2:1126624/1‑93 (MQ=255)

                                        cTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGt                                                    >  2:2025652/1‑93 (MQ=255)

                                        cTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGt                                                    <  1:20097/93‑1 (MQ=255)

                                          aTCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTa                                                  >  1:4293658/1‑93 (MQ=255)

                                          aTCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTa                                                  >  2:3006559/1‑93 (MQ=255)

                                          aTCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTa                                                  >  2:282637/1‑93 (MQ=255)

                                          aTCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTa                                                  <  2:804678/93‑1 (MQ=255)

                                           tCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAt                                                 >  1:554952/1‑93 (MQ=255)

                                            cAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtc                                                <  1:4720254/93‑1 (MQ=255)

                                             aGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcg                                               <  1:1913577/93‑1 (MQ=255)

                                             aGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcg                                               <  1:313437/93‑1 (MQ=255)

                                             aGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcg                                               <  1:313438/93‑1 (MQ=255)

                                              ggTGCTGAGCCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcgt                                              <  2:515813/93‑1 (MQ=255)

                                               gTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcgtc                                             <  1:3027514/93‑1 (MQ=255)

                                                tACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTAtcgtcg                                            <  1:911424/93‑1 (MQ=255)

                                                tACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGGCTTACTGGGATCGTTACGTTAtcgtcg                                            >  1:1902236/1‑93 (MQ=255)

                                                 aCTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCgg                                           <  2:3687206/93‑1 (MQ=255)

                                                  cTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGt                                          <  1:4279634/93‑1 (MQ=255)

                                                    gAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTc                                        <  1:976907/93‑1 (MQ=255)

                                                        cAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAg                                    >  1:3007735/1‑93 (MQ=255)

                                                         aGCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAgg                                   >  2:9286/1‑93 (MQ=255)

                                                          gCATTCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGAGATCGTTAAGTTATCGTCGGTTCGTAGGc                                  >  1:1458002/1‑93 (MQ=255)

                                                              tCACCTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGa                              <  1:3692842/93‑1 (MQ=255)

                                                                 ccTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATaa                           <  1:116920/93‑1 (MQ=255)

                                                                 ccTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATaa                           <  2:2875509/93‑1 (MQ=255)

                                                                  cTCGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAg                          <  1:3932142/93‑1 (MQ=255)

                                                                    cGTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGc                        <  1:828574/93‑1 (MQ=255)

                                                                     gTTTAACGCGCGGCTTGCCGGGGTGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCg                       <  1:1659928/93‑1 (MQ=255)

                                                                     gTTTAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCg                       >  2:4630256/1‑93 (MQ=255)

                                                                       ttAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGtt                     >  1:3169222/1‑93 (MQ=255)

                                                                       ttAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGtt                     <  1:954217/93‑1 (MQ=255)

                                                                       ttAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGtt                     <  2:4649616/93‑1 (MQ=255)

                                                                        tAACGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTc                    <  1:674567/93‑1 (MQ=255)

                                                                         aaCGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTTGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCa                   >  1:907250/1‑93 (MQ=255)

                                                                         aaCGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCa                   <  2:1932368/93‑1 (MQ=255)

                                                                          aCGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCAc                  <  2:1125251/93‑1 (MQ=255)

                                                                          aCGCGCGGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCAAGATAAGGCGTTCAc                  <  1:3999956/93‑1 (MQ=255)

                                                                            gcgcgGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCAcgc                >  2:4387703/1‑93 (MQ=255)

                                                                               cgGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAATGCGTTCAcgccgc             >  1:2820112/1‑93 (MQ=255)

                                                                               cgGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCAcgccgc             <  2:3259495/93‑1 (MQ=255)

                                                                               cgGCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCAcgccgc             >  2:1238233/1‑93 (MQ=255)

                                                                                ggCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCa            >  2:2363296/1‑93 (MQ=255)

                                                                                 gCTTGCCGGGGCGTGAAGATAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCAt           >  2:3900589/1‑93 (MQ=255)

                                                                                 gCTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCAt           >  1:2258362/1‑93 (MQ=255)

                                                                                  cTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATc          <  1:2036917/93‑1 (MQ=255)

                                                                                  cTTGCCGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATc          <  2:2672132/93‑1 (MQ=255)

                                                                                    tGCGGGGGCGTTAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTg        <  2:560803/93‑1 (MQ=255)

                                                                                      ccGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGc      <  2:2080912/93‑1 (MQ=255)

                                                                                       cGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCa     <  1:1308081/93‑1 (MQ=255)

                                                                                       cGGGGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCa     <  2:884096/93‑1 (MQ=255)

                                                                                          aGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCAttt  <  1:4400073/93‑1 (MQ=255)

CATCGTTAGCGCAGGTGTGTAAATATGCTAACCGTGTGAGCTATCAGGTACTGAACCAGCATTCACCTCGTTTAACGCGCGGCTTGCCGGAGCGTGAAGACAGCCTGGAAGAGTCTTACTGGGATCGTTAAGTTATCGTCGGTTCGTAGGCCAGATAAGGCGTTCACGCCGCATCTGGCATTT  >  NC_000913/500523‑500705

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: `‑`