BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,108,133	A→G	G208G (GGA→GGG)	cdh →	CDP‑diacylglycerol phosphotidylhydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,108,133	0	A	G	100.0%	717.3 / NA	227	G208G (GGA→GGG)	cdh	CDP‑diacylglycerol phosphotidylhydrolase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (103/124); total (103/124)

Only 100 of 238 total aligned reads displayed.

GCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTC  >  NC_000913/4108044‑4108222

  gcgTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTa                                                                                      >  2:1012466/1‑93 (MQ=255)

     tCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACgg                                                                                   >  1:188228/1‑93 (MQ=255)

     gcGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACgg                                                                                   <  2:1578698/92‑1 (MQ=255)

      cGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACggg                                                                                  >  1:4575508/1‑93 (MQ=255)

       gggTAACGGAAGGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  1:3913032/93‑1 (MQ=255)

       gggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  1:1857121/93‑1 (MQ=255)

       gggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGc                                                                                 <  1:2841963/93‑1 (MQ=255)

        ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                <  1:3220501/93‑1 (MQ=255)

        ggTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCt                                                                                <  2:4674579/93‑1 (MQ=255)

         gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               <  1:251179/93‑1 (MQ=255)

         gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               <  2:1637547/93‑1 (MQ=255)

         gTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTg                                                                               <  1:2324544/93‑1 (MQ=255)

          tAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGCACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTgg                                                                              <  1:868851/93‑1 (MQ=255)

           aaCGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGc                                                                             <  1:3523699/93‑1 (MQ=255)

             cGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGa                                                                           >  2:3323504/1‑93 (MQ=255)

             cGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGa                                                                           <  1:2566778/93‑1 (MQ=255)

              ggAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGAt                                                                          <  2:3740740/93‑1 (MQ=255)

               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         <  2:1160156/93‑1 (MQ=255)

               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         <  1:2028049/93‑1 (MQ=255)

               gAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATg                                                                         >  2:450321/1‑93 (MQ=255)

                aaaGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGTAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATgg                                                                        <  2:558701/93‑1 (MQ=255)

                   gCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTAc                                                                     <  2:3639322/93‑1 (MQ=255)

                    cGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGAGCTACGGGCTGGCGATGGTACg                                                                    >  1:1460571/1‑93 (MQ=255)

                     gAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACgg                                                                   <  2:2070246/93‑1 (MQ=255)

                     gAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACgg                                                                   >  1:485618/1‑93 (MQ=255)

                      aaCTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGc                                                                  <  2:652025/93‑1 (MQ=255)

                       aCTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCa                                                                 <  1:2997130/93‑1 (MQ=255)

                        cTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCag                                                                <  2:3001641/93‑1 (MQ=255)

                        cTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCag                                                                >  1:1075250/1‑93 (MQ=255)

                         tGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCaga                                                               <  1:552935/93‑1 (MQ=255)

                         tGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCaga                                                               >  2:4133797/1‑93 (MQ=255)

                           gTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCagagg                                                             >  2:1203354/1‑92 (MQ=255)

                            tACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTg                                                            >  1:4254202/1‑93 (MQ=255)

                              cAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGAt                                                          <  2:3077109/93‑1 (MQ=255)

                              cAACGCACCCCATTTACGATGCTGGCAGAAGAAGTCCCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGAt                                                          <  1:1201512/93‑1 (MQ=255)

                               agcGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATa                                                         <  1:1685358/91‑1 (MQ=255)

                               aaCGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATa                                                         >  1:3181772/1‑93 (MQ=255)

                                aCGCAGCCCATTTATGATGCTCGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATaa                                                        <  2:4435784/93‑1 (MQ=255)

                                    aGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCa                                                    >  1:1192247/1‑93 (MQ=255)

                                     gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   >  2:4010146/1‑93 (MQ=255)

                                     gCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAt                                                   <  1:4099707/93‑1 (MQ=255)

                                      cccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAtt                                                  <  1:2160627/93‑1 (MQ=255)

                                      cccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAtt                                                  >  1:1692380/1‑93 (MQ=255)

                                       ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 <  1:1990243/93‑1 (MQ=255)

                                       ccATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCAttt                                                 <  1:2390005/93‑1 (MQ=255)

                                        cATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTg                                                >  1:1774048/1‑93 (MQ=255)

                                        cATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTg                                                <  1:2809290/93‑1 (MQ=255)

                                        cATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGGGATAACTCATTTg                                                >  2:2702952/1‑93 (MQ=255)

                                         aTTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGt                                               >  1:2147171/1‑93 (MQ=255)

                                         aTTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGt                                               >  1:510993/1‑93 (MQ=255)

                                         aTTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGt                                               <  1:2690784/93‑1 (MQ=255)

                                          tttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTa                                              <  1:274934/93‑1 (MQ=255)

                                           ttATGATGCTGTCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAt                                             <  2:2723492/93‑1 (MQ=255)

                                           ttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAt                                             >  2:909586/1‑93 (MQ=255)

                                           ttATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTAt                                             >  2:1692852/1‑93 (MQ=255)

                                             atgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTg                                           >  1:87965/1‑93 (MQ=255)

                                             atgatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTg                                           >  1:3658515/1‑93 (MQ=255)

                                               gatgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCt                                         <  2:2211502/93‑1 (MQ=255)

                                                atgCTGGAAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTg                                        <  1:458003/93‑1 (MQ=255)

                                                 tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGGGTGACAACTCATTTGTATTTCTgg                                       >  2:3891117/1‑93 (MQ=255)

                                                 tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTgg                                       <  1:587046/93‑1 (MQ=255)

                                                 tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTgg                                       >  1:4449892/1‑93 (MQ=255)

                                                 tgCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTgg                                       >  2:1836906/1‑93 (MQ=255)

                                                  gCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGGGATAACTCATTTGTATTGCTGGc                                      >  1:942108/1‑93 (MQ=255)

                                                    tGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGa                                    >  1:417564/1‑93 (MQ=255)

                                                    tGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGa                                    >  1:662094/1‑93 (MQ=255)

                                                     ggCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGac                                   >  1:4017941/1‑93 (MQ=255)

                                                        agaagaAGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGAcaca                                >  1:1386280/1‑93 (MQ=255)

                                                         gaagaaGTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACaa                               >  2:4625033/1‑93 (MQ=255)

                                                          aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              <  2:4107194/93‑1 (MQ=255)

                                                          aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              >  2:1144543/1‑93 (MQ=255)

                                                          aagaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAAc                              <  1:290210/93‑1 (MQ=255)

                                                           agaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACg                             >  2:3447228/1‑93 (MQ=255)

                                                            gaagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGAAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGa                            <  1:3597396/93‑1 (MQ=255)

                                                            gaagTACCTGAGGCGCGGGAACAGATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGa                            <  2:238997/93‑1 (MQ=255)

                                                             aagTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaa                           <  2:313143/93‑1 (MQ=255)

                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          <  1:3375483/93‑1 (MQ=255)

                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGaaa                          <  2:696071/93‑1 (MQ=255)

                                                              agTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTACTGGCGACACAACaaaa                          >  2:4538067/1‑93 (MQ=255)

                                                               gTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAt                         <  1:1233507/93‑1 (MQ=255)

                                                               gTACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAAt                         <  1:1108603/93‑1 (MQ=255)

                                                                tACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATc                        <  2:2244392/93‑1 (MQ=255)

                                                                tACCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATc                        <  1:1006082/93‑1 (MQ=255)

                                                                 aCCTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATCTGTATTGCTGGCGACACAACGAAATCt                       >  2:3580562/1‑93 (MQ=255)

                                                                  ccTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGGACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTa                      >  2:4317690/1‑93 (MQ=255)

                                                                   cTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACTGCAGAGTGATAACTCATTTGTATTGCTGGCGGCACAACGAAATCTAc                     >  1:1638135/1‑93 (MQ=255)

                                                                   cTGAGGCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTAc                     <  1:4300707/93‑1 (MQ=255)

                                                                    tGAGGCGCGGGAACATGTGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACt                    >  2:2513624/1‑93 (MQ=255)

                                                                       ggCGCGGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGag                 >  1:3530247/1‑92 (MQ=255)

                                                                         cgcgGGAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGc               <  2:2946468/93‑1 (MQ=255)

                                                                          gcgGGAACATATGGGGCGCTACGGACTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCt              <  2:2874343/93‑1 (MQ=255)

                                                                            gggAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTa            <  1:3355925/93‑1 (MQ=255)

                                                                             ggAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTaa           <  1:3050702/93‑1 (MQ=255)

                                                                              gAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAc          >  2:2405699/1‑93 (MQ=255)

                                                                              gAACATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAAc          <  1:3105130/93‑1 (MQ=255)

                                                                                 cATATGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGt       <  2:2838726/93‑1 (MQ=255)

                                                                                  atatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTg      >  1:3423846/1‑93 (MQ=255)

                                                                                   tatGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGc     <  2:1385115/93‑1 (MQ=255)

                                                                                    atGGGGCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCt    <  2:1933890/93‑1 (MQ=255)

                                                                                      ggggCGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTc  <  2:534542/93‑1 (MQ=255)

GCGCGTCGGGTAACGGAAAGCGAACTGGTACAACGCAGCCCATTTATGATGCTGGCAGAAGAAGTACCTGAGGCGCGGGAACATATGGGACGCTACGGGCTGGCGATGGTACGGCAGAGTGATAACTCATTTGTATTGCTGGCGACACAACGAAATCTACTGACGCTTAACCGTGCTTC  >  NC_000913/4108044‑4108222

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