BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,808,348	G→A	pseudogene (33/1178 nt)	ygaY →	putative transporter YgaY, N‑terminal fragment

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,808,348	0	G	A	100.0%	1390.8 / NA	326	pseudogene (33/1178 nt)	ygaY	putative transporter YgaY, N‑terminal fragment
Reads supporting (aligned to +/- strand): ref base G (0/0); new base A (191/135); total (191/135)

Only 100 of 326 total aligned reads displayed.

TCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTT  >  NC_000913/2808250‑2808445

tCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACt                                                                                                 <  2:4258997/101‑1 (MQ=255)

tCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACt                                                                                                 <  2:4259202/101‑1 (MQ=255)

tCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACt                                                                                                 <  1:4259122/101‑1 (MQ=255)

tCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACt                                                                                                 <  1:4259029/101‑1 (MQ=255)

tCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACt                                                                                                 >  1:4259412/1‑101 (MQ=255)

 cccATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTg                                                                                                >  2:4259414/1‑101 (MQ=255)

   cATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGAt                                                                                              >  1:4259416/1‑101 (MQ=255)

    atatTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATc                                                                                             >  2:4259422/1‑101 (MQ=255)

    atatTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATc                                                                                             >  1:4259419/1‑101 (MQ=255)

      atTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGt                                                                                           >  1:4259428/1‑101 (MQ=255)

       ttCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTg                                                                                          >  1:4259435/1‑101 (MQ=255)

       ttCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTg                                                                                          >  2:4259430/1‑101 (MQ=255)

        tcatcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGc                                                                                         <  1:4258923/101‑1 (MQ=255)

          atcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTg                                                                                       <  1:4258992/101‑1 (MQ=255)

          atcaACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTg                                                                                       >  1:4259438/1‑101 (MQ=255)

             aaCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATg                                                                                    >  2:4259443/1‑101 (MQ=255)

             aaCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATg                                                                                    >  1:4259442/1‑101 (MQ=255)

              aCATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGt                                                                                   <  2:4258748/101‑1 (MQ=255)

               cATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTc                                                                                  <  1:4259146/101‑1 (MQ=255)

               cATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTc                                                                                  >  1:4259449/1‑101 (MQ=255)

                aTCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCt                                                                                 >  2:4259453/1‑101 (MQ=255)

                  ccGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTAt                                                                               >  1:4259454/1‑101 (MQ=255)

                   cGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATc                                                                              >  2:4259457/1‑101 (MQ=255)

                   cGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATc                                                                              <  2:4258938/101‑1 (MQ=255)

                     cTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGc                                                                            >  1:4259461/1‑101 (MQ=255)

                     cTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGc                                                                            <  2:4258651/101‑1 (MQ=255)

                     cTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGc                                                                            >  2:4259463/1‑101 (MQ=255)

                      tattatTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGcc                                                                           <  1:4258635/101‑1 (MQ=255)

                       attattGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCa                                                                          >  2:4259468/1‑101 (MQ=255)

                          attGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCg                                                                       >  1:4259473/1‑101 (MQ=255)

                          attGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCg                                                                       <  1:4258972/101‑1 (MQ=255)

                           ttGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCgg                                                                      >  2:4259476/1‑101 (MQ=255)

                           ttGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCgg                                                                      <  2:4258978/101‑1 (MQ=255)

                            tGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGt                                                                     >  1:4259478/1‑101 (MQ=255)

                              aTTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCt                                                                   >  2:4259481/1‑101 (MQ=255)

                              aTTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCt                                                                   >  2:4259484/1‑101 (MQ=255)

                               tttCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTg                                                                  >  1:4259486/1‑101 (MQ=255)

                               tttCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTg                                                                  <  2:4259019/101‑1 (MQ=255)

                                ttCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTgg                                                                 >  1:4259489/1‑101 (MQ=255)

                                  ccAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCg                                                               <  1:4258835/101‑1 (MQ=255)

                                   cAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCgg                                                              >  2:4259495/1‑101 (MQ=255)

                                       ttAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCgg                                                              >  2:4259500/1‑97 (MQ=255)

                                       ttAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGc                                                          <  2:4259283/101‑1 (MQ=255)

                                       ttAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGc                                                          >  1:4259498/1‑101 (MQ=255)

                                         aaTCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCa                                                        <  1:4259096/101‑1 (MQ=255)

                                          atcatcACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAg                                                       >  2:4259505/1‑101 (MQ=255)

                                          atcatcACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAg                                                       >  2:4259508/1‑101 (MQ=255)

                                           tcatcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGt                                                      >  1:4259509/1‑101 (MQ=255)

                                            catcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTa                                                     <  2:4258857/101‑1 (MQ=255)

                                            catcaCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTa                                                     <  1:4258774/101‑1 (MQ=255)

                                                aCCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTa                                                 >  1:4259515/1‑101 (MQ=255)

                                                 ccTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTAt                                                <  2:4258796/101‑1 (MQ=255)

                                                   tgatgaACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTa                                              >  2:4259516/1‑101 (MQ=255)

                                                       gaACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGcc                                          >  2:4259520/1‑101 (MQ=255)

                                                        aaCAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTATATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGccc                                         >  1:4259522/1‑101 (MQ=255)

                                                         aCAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCa                                        <  2:4259174/101‑1 (MQ=255)

                                                          cAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAg                                       >  2:4259525/1‑101 (MQ=255)

                                                            aaaaTAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCGGTAACTATTACGCCCAGcc                                     >  2:4259532/1‑101 (MQ=255)

                                                            aaaaTAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGcc                                     <  1:4259007/101‑1 (MQ=255)

                                                            aaaaTAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGcc                                     >  1:4259531/1‑101 (MQ=255)

                                                             aaaTAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                    <  1:4258903/101‑1 (MQ=255)

                                                             aaaTAATGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCa                                    >  1:4259537/1‑101 (MQ=255)

                                                              aataatGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCAt                                   >  1:4259540/1‑101 (MQ=255)

                                                                taatGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTg                                 >  2:4259545/1‑101 (MQ=255)

                                                                taatGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTg                                 >  1:4259546/1‑101 (MQ=255)

                                                                 aatGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCgg                                                              >  2:4259550/1‑71 (MQ=255)

                                                                  atGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCt                               <  2:4258924/101‑1 (MQ=255)

                                                                   tGACTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTc                              >  1:4259553/1‑101 (MQ=255)

                                                                      cTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGac                           >  1:4259556/1‑101 (MQ=255)

                                                                      cTAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGac                           <  1:4259243/101‑1 (MQ=255)

                                                                       tAAACCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGaca                          <  2:4259261/101‑1 (MQ=255)

                                                                        aaaCCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGacac                         >  1:4259559/1‑101 (MQ=255)

                                                                          aCCTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCa                       >  1:4259562/1‑101 (MQ=255)

                                                                           ccTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACGCCAt                      <  2:4259057/101‑1 (MQ=255)

                                                                           ccTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCAt                      <  2:4259121/101‑1 (MQ=255)

                                                                            cTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATc                     <  1:4259063/101‑1 (MQ=255)

                                                                            cTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATc                     <  1:4259156/101‑1 (MQ=255)

                                                                            cTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATc                     >  2:4259572/1‑101 (MQ=255)

                                                                            cTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATc                     >  2:4259566/1‑101 (MQ=255)

                                                                            cTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATc                     >  1:4259570/1‑101 (MQ=255)

                                                                            cTAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATc                     <  1:4258815/101‑1 (MQ=255)

                                                                             tAATCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcg                    <  2:4258842/101‑1 (MQ=255)

                                                                              aaTCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcgc                   >  1:4259580/1‑101 (MQ=255)

                                                                              aaTCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcgc                   <  2:4259149/101‑1 (MQ=255)

                                                                              aaTCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcgc                   >  1:4259577/1‑101 (MQ=255)

                                                                              aaTCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcgc                   <  2:4259088/101‑1 (MQ=255)

                                                                                tCATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcgcgc                 <  1:4259113/101‑1 (MQ=255)

                                                                                 cATGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATcgcgcg                <  1:4258982/101‑1 (MQ=255)

                                                                                  aTGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGt               >  2:4259585/1‑101 (MQ=255)

                                                                                  aTGAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGt               >  1:4259586/1‑101 (MQ=255)

                                                                                    gAGCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTaa             >  2:4259591/1‑101 (MQ=255)

                                                                                      gCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACt           <  2:4259170/101‑1 (MQ=255)

                                                                                      gCTTAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACt           >  1:4259594/1‑101 (MQ=255)

                                                                                        ttAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACttt         <  1:4259208/101‑1 (MQ=255)

                                                                                        ttAGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACttt         >  2:4259595/1‑101 (MQ=255)

                                                                                          aGCCCGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTc       <  1:4258936/101‑1 (MQ=255)

                                                                                            cccGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTccc     >  1:4259598/1‑101 (MQ=255)

                                                                                             ccGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCt    <  2:4259222/101‑1 (MQ=255)

                                                                                             ccGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCt    <  2:4259078/101‑1 (MQ=255)

                                                                                              cGGCACTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCtt   <  1:4259191/101‑1 (MQ=255)

TCCCATATTCATCAACATCCGCTATTATTGATTTCCAGCTTAATCATCACCTGATGAACAAAAATAATGACTAAACCTAATCATGAGCTTAGCCCGGCGCTGATCGTGCTGATGTCTATCGCCACCGGTCTGGCGGTAGCCAGTAACTATTACGCCCAGCCATTGCTCGACACCATCGCGCGTAACTTTTCCCTTT  >  NC_000913/2808250‑2808445

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