BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,932	A→G	A52A (GCT→GCC)	b4161 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,932	0	A	G	99.7%	985.2 / ‑5.1	321	A52A (GCT→GCC)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (151/168); total (151/170)
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.68e-01

Only 100 of 347 total aligned reads displayed.

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                         <  1:5280094/93‑1 (MQ=255)

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                         <  1:1634698/93‑1 (MQ=255)

 cGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCg                                                                                        <  1:7433357/93‑1 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCCGCGt                                                                                       >  2:1173428/1‑93 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                    <  1:1407704/93‑1 (MQ=255)

     aGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCa                                                                                    >  1:4680256/1‑93 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   >  2:3682739/1‑93 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   <  2:822221/93‑1 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATa                                                                                  >  2:6262223/1‑93 (MQ=255)

        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAc                                                                                 <  1:1005913/93‑1 (MQ=255)

        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAc                                                                                 <  2:339937/93‑1 (MQ=255)

        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGGACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAc                                                                                 >  2:4553764/1‑93 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                >  1:2699574/1‑93 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                <  1:590775/93‑1 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                <  2:5783159/93‑1 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGGACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTTCACCTCCGCGTGCCTAcc                                                                                >  1:2788722/1‑93 (MQ=255)

          aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCa                                                                               >  1:4590645/1‑93 (MQ=255)

          aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCa                                                                               <  2:1623274/93‑1 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              >  2:15886/1‑93 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              >  2:605977/1‑93 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              <  2:6477554/93‑1 (MQ=255)

             cgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAg                                                                            <  2:672432/93‑1 (MQ=255)

             cgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAg                                                                            >  1:4227249/1‑93 (MQ=255)

             cgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAg                                                                            <  2:2760742/93‑1 (MQ=255)

              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                           <  2:4780799/93‑1 (MQ=255)

              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                           <  1:663219/93‑1 (MQ=255)

               gTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCg                                                                          >  2:5128661/1‑93 (MQ=255)

                tCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCgg                                                                         >  1:993032/1‑93 (MQ=255)

                 cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  2:332066/93‑1 (MQ=255)

                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtga                                                                     <  1:4078772/93‑1 (MQ=255)

                     ggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgat                                                                    <  2:906323/93‑1 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   >  2:3522454/1‑93 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   >  2:4891678/1‑93 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCAGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   >  2:1336068/1‑93 (MQ=255)

                      gTTACCAGCGAACGGATGGGACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCATATCTGCGTGCATACCAAAGCGGCtgatg                                                                   >  2:7251807/1‑93 (MQ=255)

                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  <  2:7638856/93‑1 (MQ=255)

                          ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATTCCAAAGCGGCTGATGACGa                                                               <  1:7660737/93‑1 (MQ=255)

                          ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGa                                                               <  2:2308331/93‑1 (MQ=255)

                          ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGa                                                               <  2:568840/93‑1 (MQ=255)

                          ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGa                                                               <  1:7552790/93‑1 (MQ=255)

                              cGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAcc                                                           <  2:1985959/93‑1 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                          >  1:2980277/1‑93 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCt                                                          >  1:5034338/1‑93 (MQ=255)

                               gAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAACt                                                          >  2:6443431/1‑93 (MQ=255)

                                 aCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtc                                                        <  1:3421999/93‑1 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                       <  1:1492113/93‑1 (MQ=255)

                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      >  1:4442275/1‑93 (MQ=255)

                                   ggATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcat                                                      <  1:5541734/93‑1 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                     >  1:3986179/1‑93 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                     >  1:2290353/1‑93 (MQ=255)

                                     aTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                    >  1:5618247/1‑93 (MQ=255)

                                     aTGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatca                                                    <  2:1737565/93‑1 (MQ=255)

                                      tGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAg                                                   >  1:2495617/1‑93 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:4091445/93‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:5227296/93‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  1:7105915/93‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  1:6742960/93‑1 (MQ=255)

                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGc                                             >  2:3858831/1‑93 (MQ=255)

                                            cgCGAATATTGCAGCGATGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGc                                             <  2:2961929/92‑1 (MQ=255)

                                               gAATATTGCAGCGGTGAACGTCGCCCTCGGCGGATTCCACATCGGCGTGCATACCAAAACGGCTGATGACGAAACCTTAAAAAGGCTCGCCaa                                          >  1:6449251/1‑93 (MQ=255)

                                                aaTATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaaa                                         >  1:6272754/1‑93 (MQ=255)

                                                aaTATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaaa                                         <  1:6043456/93‑1 (MQ=255)

                                                 atatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAAc                                        >  2:7848460/1‑93 (MQ=255)

                                                  tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       <  1:44283/93‑1 (MQ=255)

                                                  tatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACa                                       >  1:2461773/1‑93 (MQ=255)

                                                      tcAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAtt                                   <  1:1976625/92‑1 (MQ=255)

                                                       cAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTg                                  <  2:5419297/93‑1 (MQ=255)

                                                        cgCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgt                                 <  1:4886905/92‑1 (MQ=255)

                                                         gCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtc                                >  2:6313932/1‑93 (MQ=255)

                                                          cGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcg                               <  1:3877589/93‑1 (MQ=255)

                                                           ggTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgt                              >  2:2435900/1‑93 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                             <  1:3504644/93‑1 (MQ=255)

                                                              gAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTcgtcgt                           <  2:5693098/93‑1 (MQ=255)

                                                                aCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAg                         >  1:5244504/1‑93 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  1:2083537/93‑1 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCACATAGGGGTGCATACCAGCGCGGCTGATGACGCTACCTTCCTCAGGCTCGACAAACAGATTGTCGGCGGAggg                       >  2:3151316/1‑91 (MQ=255)

                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                      >  1:1097866/1‑93 (MQ=255)

                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                      >  1:3829059/1‑93 (MQ=255)

                                                                   tCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCg                      >  1:5134763/1‑93 (MQ=255)

                                                                     gCCATCGGCGGATTCCACATCGGCGTGCCTACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTTGCTGTCGTCggc                    >  1:2484757/1‑92 (MQ=255)

                                                                      ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGc                   >  1:6880237/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACCTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                  <  2:7044916/93‑1 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                 >  2:2572582/1‑93 (MQ=255)

                                                                        aTCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCtt                 >  1:3758875/1‑93 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCttt                >  1:3092788/1‑92 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc                <  2:4652125/93‑1 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc                >  1:7912723/1‑93 (MQ=255)

                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             <  1:7329094/93‑1 (MQ=255)

                                                                             cggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt            >  2:7399161/1‑93 (MQ=255)

                                                                             cggATTCCACATCGGCGTGCATAACAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCt            <  2:4332526/93‑1 (MQ=255)

                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt           <  1:260467/93‑1 (MQ=255)

                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt           >  1:828379/1‑93 (MQ=255)

                                                                               gATTCCACATCCGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          <  2:5902639/93‑1 (MQ=255)

                                                                                aTTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAg         <  2:1952931/93‑1 (MQ=255)

                                                                                aTTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAg         >  2:1830581/1‑93 (MQ=255)

                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        <  1:3219894/93‑1 (MQ=255)

                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        <  2:5549231/93‑1 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       <  1:3174474/93‑1 (MQ=255)

                                                                                     acaTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtt    >  2:6647753/1‑93 (MQ=255)

                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  >  2:6836938/1‑93 (MQ=255)

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

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