BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	991,040	A→G	N399S (AAC→AGC)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	991,040	0	A	G	95.4%	385.6 / 5.6	131	N399S (AAC→AGC)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand): ref base A (1/5); new base G (48/77); total (49/82)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 135 total aligned reads displayed.

CATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGT  >  NC_000913/990951‑991129

cATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCtt                                                                                        <  2:1928564/93‑1 (MQ=255)

  ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTcc                                                                                      <  1:1202804/93‑1 (MQ=255)

  ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTcc                                                                                      >  1:171197/1‑93 (MQ=255)

  ttGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTcc                                                                                      <  1:1571973/93‑1 (MQ=255)

   tGAGATGAACAACTTCTACCCCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCa                                                                                     <  2:3483697/93‑1 (MQ=255)

   tGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACGCCCTGCTTGGCGAAGAAAGCTTCCa                                                                                     <  2:2000103/93‑1 (MQ=255)

      gaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGaa                                                                                  >  1:1288510/1‑93 (MQ=255)

      gaTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGaa                                                                                  >  2:1018573/1‑93 (MQ=255)

       aTGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGaaa                                                                                 <  2:3466934/93‑1 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                >  1:334278/1‑93 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                >  1:347401/1‑93 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                <  1:489029/93‑1 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                >  1:972766/1‑93 (MQ=255)

        tGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAg                                                                                <  2:1457135/93‑1 (MQ=255)

         gAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAgg                                                                               <  1:3009157/93‑1 (MQ=255)

          aacaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAggg                                                                              >  2:4066488/1‑93 (MQ=255)

           acaacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGa                                                                             >  2:4603427/1‑93 (MQ=255)

            caacTTCTACACCCTGACCGTTTACGAGAAGGGTGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGAt                                                                            >  2:4067425/1‑93 (MQ=255)

            caacTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGAt                                                                            <  1:884678/93‑1 (MQ=255)

              acTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGc                                                                          >  2:3407476/1‑93 (MQ=255)

               cTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCa                                                                         <  2:2418263/93‑1 (MQ=255)

               cTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCa                                                                         <  1:4414008/93‑1 (MQ=255)

                 tCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGc                                                                       <  1:4537632/93‑1 (MQ=255)

                  cTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCt                                                                      <  1:4005379/93‑1 (MQ=255)

                   tACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCtt                                                                     <  2:784916/93‑1 (MQ=255)

                   tACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCtt                                                                     >  1:2801771/1‑93 (MQ=255)

                     cacCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTa                                                                   >  1:2817124/1‑93 (MQ=255)

                      acCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTAt                                                                  <  1:4525111/93‑1 (MQ=255)

                       ccgtGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTGCCAGAAAGGCATGCCGCTTTAtt                                                                 >  2:1316443/4‑93 (MQ=255)

                         cTGACGGTGTACGAGAAGGGCGCGGAAGTGATTCGCATTATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTc                                                               <  1:748023/93‑1 (MQ=255)

                         cTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTc                                                               <  2:551708/93‑1 (MQ=255)

                         cTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTc                                                               >  2:4194501/1‑93 (MQ=255)

                           gACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGTTGCAGCTTTATTTCGa                                                             <  2:4507311/93‑1 (MQ=255)

                             ccGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGc                                                           >  1:559025/1‑93 (MQ=255)

                              cGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCg                                                          <  1:1067747/93‑1 (MQ=255)

                               gTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGt                                                         >  2:773605/1‑93 (MQ=255)

                               gTTTACGAGAAGGGCGCAGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGt                                                         <  2:3687990/93‑1 (MQ=255)

                                tttACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTc                                                        <  2:1175306/93‑1 (MQ=255)

                                 ttACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCa                                                       <  1:1067234/93‑1 (MQ=255)

                                  tACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCat                                                      <  2:1632864/93‑1 (MQ=255)

                                   aCGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatg                                                     <  1:941982/93‑1 (MQ=255)

                                    cGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatga                                                    <  2:1814003/93‑1 (MQ=255)

                                    cGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatga                                                    >  1:3193212/1‑93 (MQ=255)

                                      ggaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatgatg                                                  <  1:1545853/92‑1 (MQ=255)

                                      agaAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCatgatg                                                  >  2:2425874/1‑93 (MQ=255)

                                         aGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTa                                               >  2:3801964/1‑93 (MQ=255)

                                          gggCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAg                                              <  2:69768/93‑1 (MQ=255)

                                           ggCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGt                                             <  2:3302510/93‑1 (MQ=255)

                                             cgcgGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGc                                           <  2:3858246/93‑1 (MQ=255)

                                              gcgGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCa                                          <  1:1616480/93‑1 (MQ=255)

                                              gcgGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCa                                          >  2:315565/1‑93 (MQ=255)

                                               cgGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAg                                         >  1:2322810/1‑93 (MQ=255)

                                               cgGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAg                                         <  2:3546484/93‑1 (MQ=255)

                                                 gAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCg                                       >  1:1692651/1‑93 (MQ=255)

                                                 gAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCg                                       >  1:899323/1‑93 (MQ=255)

                                                    gTGATTCGCATGCTAGACACCCTGCTTGGCAAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGAcc                                    <  1:3332925/93‑1 (MQ=255)

                                                    gTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGAcc                                    >  2:1082704/1‑93 (MQ=255)

                                                    gTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGAcc                                    <  2:2311951/93‑1 (MQ=255)

                                                      gATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGCAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtg                                  >  2:2558405/1‑93 (MQ=255)

                                                      gATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtg                                  <  1:2433273/93‑1 (MQ=255)

                                                      gATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtg                                  >  2:3760682/1‑93 (MQ=255)

                                                      gATTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtg                                  >  2:394191/1‑93 (MQ=255)

                                                       aTTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtgt                                 <  2:2610900/93‑1 (MQ=255)

                                                       aTTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtgt                                 <  1:2505348/93‑1 (MQ=255)

                                                       aTTCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCtgt                                 <  1:802507/93‑1 (MQ=255)

                                                         tCGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                               <  1:3353583/93‑1 (MQ=255)

                                                         tCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTga                               <  2:3319044/93‑1 (MQ=255)

                                                          cGCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgac                              <  1:1188231/93‑1 (MQ=255)

                                                           gCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacg                             >  1:818822/1‑93 (MQ=255)

                                                           gCATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacg                             <  1:3520757/93‑1 (MQ=255)

                                                            cATGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacga                            <  1:2229289/93‑1 (MQ=255)

                                                             aTGATCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacgac                           <  1:4710122/93‑1 (MQ=255)

                                                             aTGATCCACACCCTGCTTGGAGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTgacgac                           <  1:515631/93‑1 (MQ=255)

                                                                aTCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATTATGGTAGTGCAGCGTCCTGTGACGACttt                        <  2:1108497/93‑1 (MQ=255)

                                                                 tCCACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTtg                       >  2:3961874/1‑93 (MQ=255)

                                                                  ccACACCCTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTCACGACTTtgt                      >  2:3125839/1‑93 (MQ=255)

                                                                       cccTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAgg                 >  2:3957931/1‑93 (MQ=255)

                                                                       cccTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAgg                 <  1:3196251/93‑1 (MQ=255)

                                                                       cccTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAgg                 >  1:2211945/1‑93 (MQ=255)

                                                                        ccTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                <  1:2844065/93‑1 (MQ=255)

                                                                        ccTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                <  1:4586881/93‑1 (MQ=255)

                                                                        ccTGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGc                <  1:2549650/93‑1 (MQ=255)

                                                                          tGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGa              >  1:2870675/1‑93 (MQ=255)

                                                                          tGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGa              >  2:4737640/1‑93 (MQ=255)

                                                                          tGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGa              >  1:4736181/1‑93 (MQ=255)

                                                                          tGCTTGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGa              <  2:1663691/93‑1 (MQ=255)

                                                                             ttGGCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATgg           >  1:2056616/1‑93 (MQ=255)

                                                                               tgCGAAAAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCTTCATTATGGTAGTGCATCGACCTGTGACTACTTTGTGCAGGCGATGGaa         <  1:2384939/92‑1 (MQ=255)

                                                                                gCGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAg        <  1:2958969/93‑1 (MQ=255)

                                                                                 cGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGa       <  2:1611099/93‑1 (MQ=255)

                                                                                 cGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGa       <  1:4021826/93‑1 (MQ=255)

                                                                                 cGAAGAAAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGa       >  2:2841262/1‑93 (MQ=255)

                                                                                  gaagaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGAt      <  1:1113975/93‑1 (MQ=255)

                                                                                   aagaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGCCGACTTTGTGCAGGCGATGGTAGATg     >  2:2545078/1‑93 (MQ=255)

                                                                                   aagaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATg     >  2:2066611/1‑93 (MQ=255)

                                                                                   aagaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATg     <  2:37314/93‑1 (MQ=255)

                                                                                    agaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATgt    >  2:574130/1‑92 (MQ=255)

                                                                                    agaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGc    <  2:3744189/93‑1 (MQ=255)

                                                                                    agaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGc    <  2:3308554/93‑1 (MQ=255)

                                                                                     gaaAGCTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCg   >  2:3805121/1‑93 (MQ=255)

CATTGAGATGAACAACTTCTACACCCTGACCGTTTACGAGAAGGGCGCGGAAGTGATTCGCATGATCCACACCCTGCTTGGCGAAGAAAACTTCCAGAAAGGGATGCAGCTTTATTTCGAGCGTCATGATGGTAGTGCAGCGACCTGTGACGACTTTGTGCAGGCGATGGAAGATGCGT  >  NC_000913/990951‑991129

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