BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,853,551	T→C	Q16R (CAG→CGG)	yidG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,853,551	0	T	C	90.4%	786.8 / 61.8	292	Q16R (CAG→CGG)	yidG	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (12/16); new base C (122/142); total (134/158)
Fisher's exact test for biased strand distribution p-value = 8.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 296 total aligned reads displayed.

CCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAA  >  NC_000913/3853459‑3853643

ccGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCt                                                                                              >  1:2687432/1‑93 (MQ=255)

  gCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGc                                                                                            >  1:962750/1‑93 (MQ=255)

   ccTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCa                                                                                           <  1:4600565/93‑1 (MQ=255)

    cTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAg                                                                                          <  1:3995638/93‑1 (MQ=255)

    cTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAg                                                                                          <  2:4077573/93‑1 (MQ=255)

       gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTcc                                                                                       >  2:882948/1‑93 (MQ=255)

       gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTcc                                                                                       <  1:2031450/93‑1 (MQ=255)

       gTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTcc                                                                                       <  2:4730713/93‑1 (MQ=255)

         gCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTg                                                                                     <  1:770257/93‑1 (MQ=255)

         gCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTg                                                                                     <  1:3519757/93‑1 (MQ=255)

          ccAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTgg                                                                                    >  1:4818405/1‑93 (MQ=255)

          ccAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTgg                                                                                    <  2:731363/93‑1 (MQ=255)

           cAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTggg                                                                                   >  1:2532572/1‑93 (MQ=255)

            aGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGt                                                                                  <  1:1670614/93‑1 (MQ=255)

            aGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGt                                                                                  <  1:6048586/93‑1 (MQ=255)

            aGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGt                                                                                  <  1:4138936/93‑1 (MQ=255)

             gTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTc                                                                                 >  1:2899716/1‑93 (MQ=255)

              ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCg                                                                                >  2:5579499/1‑93 (MQ=255)

              ttGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCg                                                                                >  1:30948/1‑93 (MQ=255)

                gtgtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGc                                                                              <  1:6419611/93‑1 (MQ=255)

                 tgtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCa                                                                             >  1:3403483/1‑93 (MQ=255)

                  gtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCaa                                                                            <  2:1131576/93‑1 (MQ=255)

                  gtTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCaa                                                                            <  2:5238333/93‑1 (MQ=255)

                   tttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAAt                                                                           >  2:3255415/1‑93 (MQ=255)

                   tttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAAt                                                                           <  1:1472606/93‑1 (MQ=255)

                    ttGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATg                                                                          >  1:880396/1‑93 (MQ=255)

                      gATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCg                                                                        <  1:1624940/93‑1 (MQ=255)

                        tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTc                                                                      <  2:1315836/93‑1 (MQ=255)

                        tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTc                                                                      >  2:4241553/1‑93 (MQ=255)

                        tAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTc                                                                      <  1:6508132/93‑1 (MQ=255)

                          gCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTg                                                                    <  2:31972/93‑1 (MQ=255)

                           ccAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTgg                                                                   <  1:4263060/93‑1 (MQ=255)

                           ccAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTgg                                                                   <  1:5826322/93‑1 (MQ=255)

                            cAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGc                                                                  <  2:2701280/93‑1 (MQ=255)

                             aaCGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCt                                                                 >  1:2327704/1‑93 (MQ=255)

                                gCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTg                                                              >  1:5026889/1‑93 (MQ=255)

                                 ccATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGc                                                             <  2:750876/93‑1 (MQ=255)

                                 ccATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTga                                                             <  2:4637282/93‑2 (MQ=255)

                                  cATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCg                                                            >  2:477287/1‑93 (MQ=255)

                                  cATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCg                                                            >  2:7079841/1‑93 (MQ=255)

                                   aTCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCgg                                                           >  2:5794061/1‑93 (MQ=255)

                                    tCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGc                                                          >  2:3368671/1‑93 (MQ=255)

                                    tCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGc                                                          >  2:6163387/1‑93 (MQ=255)

                                     cAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCt                                                         <  2:3804836/93‑1 (MQ=255)

                                      aaCGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTa                                                        <  2:2259395/93‑1 (MQ=255)

                                        cgcgCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATc                                                      <  1:5201904/93‑1 (MQ=255)

                                         gcgcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATcc                                                     >  1:7252421/1‑93 (MQ=255)

                                         gcgcCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATcc                                                     <  1:4194345/93‑1 (MQ=255)

                                            ccGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGc                                                  <  2:7464997/93‑1 (MQ=255)

                                              gTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCAt                                                >  2:188643/1‑93 (MQ=255)

                                              gTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCAt                                                >  2:293767/1‑93 (MQ=255)

                                              gTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCAt                                                >  1:7417952/1‑93 (MQ=255)

                                               tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATa                                               >  2:5101891/1‑93 (MQ=255)

                                               tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATa                                               <  2:4294831/93‑1 (MQ=255)

                                               tAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATa                                               >  1:1180255/1‑93 (MQ=255)

                                                aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTTCGGCTATCCGGCATac                                              >  2:5356347/1‑93 (MQ=255)

                                                aGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATac                                              <  2:7280661/93‑1 (MQ=255)

                                                 gCCCAGCATGGTACGGAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCAGGCAGaca                                             >  1:5554182/1‑93 (MQ=255)

                                                   ccAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATacaac                                           <  1:3532036/93‑1 (MQ=255)

                                                    cAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATacaaca                                          >  2:6303048/1‑93 (MQ=255)

                                                     aGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAacac                                         >  2:1046983/1‑93 (MQ=255)

                                                     aGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAacac                                         >  2:3639772/1‑93 (MQ=255)

                                                       cATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCa                                       >  1:6950072/1‑93 (MQ=255)

                                                           gTACGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTcc                                   <  1:2164745/93‑1 (MQ=255)

                                                             aCGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCa                                 >  1:6597498/1‑93 (MQ=255)

                                                             aCGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCa                                 <  2:6846320/93‑1 (MQ=255)

                                                              cGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCAt                                >  2:954323/1‑93 (MQ=255)

                                                              cGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCAt                                <  2:2819524/93‑1 (MQ=255)

                                                              cGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCAt                                <  1:2207058/93‑1 (MQ=255)

                                                              cGAAACCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCAt                                <  2:1069669/93‑1 (MQ=255)

                                                                aaaCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                              <  1:1778880/93‑1 (MQ=255)

                                                                aaaCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATaa                              <  1:5112950/93‑1 (MQ=255)

                                                                  aCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                            >  2:5896410/1‑93 (MQ=255)

                                                                  aCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                            >  1:7104632/1‑93 (MQ=255)

                                                                  aCCACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGCCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACg                            >  2:576626/1‑93 (MQ=255)

                                                                   ccACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGa                           <  1:3706234/93‑1 (MQ=255)

                                                                    cACGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGAt                          <  2:2490801/93‑1 (MQ=255)

                                                                     aCGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATc                         >  1:1544189/1‑93 (MQ=255)

                                                                      cGCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCa                        <  1:6732017/93‑1 (MQ=255)

                                                                       gCCAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATcgc                       <  2:2623975/93‑3 (MQ=255)

                                                                        ccAGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGACCCATAACGATCAcc                      <  2:311772/93‑1 (MQ=255)

                                                                          aGTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGc                    <  2:6481261/93‑1 (MQ=255)

                                                                           gTGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGcc                   >  2:7218328/1‑93 (MQ=255)

                                                                            tGATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCa                  <  2:6813683/93‑1 (MQ=255)

                                                                             gATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc                 <  1:7388776/93‑1 (MQ=255)

                                                                             gATGTGCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCAc                 <  1:5365715/93‑1 (MQ=255)

                                                                               tgtgCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGa               >  2:5426560/1‑93 (MQ=255)

                                                                                gtgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGAc              <  2:2374885/93‑1 (MQ=255)

                                                                                gtgCGCTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGAc              <  2:2337043/93‑1 (MQ=255)

                                                                                  gcgcTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACCCCAGTCCCATAACGATCACCGCCACGACCa            <  1:6986720/93‑1 (MQ=255)

                                                                                  gcgcTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCa            >  2:457080/1‑93 (MQ=255)

                                                                                   cgcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt           <  2:3988952/93‑1 (MQ=255)

                                                                                   cgcTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCAt           <  2:7037838/93‑1 (MQ=255)

                                                                                    gcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc          >  1:4406643/1‑93 (MQ=255)

                                                                                    gcTCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATc          <  1:1113456/93‑1 (MQ=255)

                                                                                      tCCGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCaa        <  2:1254828/93‑1 (MQ=255)

                                                                                       ccGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAat       <  1:6799400/93‑1 (MQ=255)

                                                                                        cGGCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAata      <  1:3785584/93‑1 (MQ=255)

                                                                                         ggCCGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAatat     >  1:6246643/1‑93 (MQ=255)

                                                                                            tGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTaa  <  2:1968335/93‑1 (MQ=255)

CCGCCTGGTGCCAGTTGTGTTTGATAGCCAACGCCATCAACGCGCCGTAGCCCAGCATGGTACGAAACCACGCCAGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAA  >  NC_000913/3853459‑3853643

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