BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,810,674	T→C	L106P (CTG→CCG)	b3639 →	Dfp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,810,674	0	T	C	99.5%	668.9 / NA	213	L106P (CTG→CCG)	b3639	Dfp
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (116/95); minor base G (1/0); total (118/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

TTGAGCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAG  >  NC_000913/3810584‑3810762

  gAGCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTa                                                                                      >  1:5486216/1‑93 (MQ=255)

   aGCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTAt                                                                                     >  2:612606/1‑93 (MQ=255)

    gCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATc                                                                                    <  1:2730405/93‑1 (MQ=255)

    gCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATc                                                                                    <  1:4532731/93‑1 (MQ=255)

    gCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATc                                                                                    >  1:4123903/1‑93 (MQ=255)

         gTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATcgacga                                                                               <  1:3513022/93‑1 (MQ=255)

          tAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAt                                                                              <  1:1082203/93‑1 (MQ=255)

          tAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAt                                                                              >  2:2557504/1‑93 (MQ=255)

          tAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAt                                                                              >  1:2509984/1‑93 (MQ=255)

           aaaTGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAtt                                                                             >  1:5325204/1‑93 (MQ=255)

           aaaTGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAtt                                                                             <  2:3169361/93‑1 (MQ=255)

            aaTGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAttt                                                                            <  2:4338736/93‑1 (MQ=255)

            aaTGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAttt                                                                            >  2:5198418/1‑93 (MQ=255)

            aaTGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGAttt                                                                            >  1:559240/1‑93 (MQ=255)

             aTGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTg                                                                           >  1:22738/1‑93 (MQ=255)

              tGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGt                                                                          <  2:2879293/93‑1 (MQ=255)

              tGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGt                                                                          >  2:257994/1‑93 (MQ=255)

               gggCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTc                                                                         >  1:3176615/1‑93 (MQ=255)

                ggCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCt                                                                        <  1:5585007/93‑1 (MQ=255)

                 gCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCGGGTATCGACGATTTGTCTg                                                                       >  2:3937496/1‑93 (MQ=255)

                  cTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTgg                                                                      <  2:5399828/93‑1 (MQ=255)

                  cTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTgg                                                                      <  1:3299154/93‑1 (MQ=255)

                   tGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGc                                                                     <  1:5440870/93‑1 (MQ=255)

                   tGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGc                                                                     <  1:5496685/93‑1 (MQ=255)

                    gATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCt                                                                    >  2:4377905/1‑93 (MQ=255)

                    gATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCt                                                                    <  2:1521459/93‑1 (MQ=255)

                    gATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCt                                                                    <  1:4845452/93‑1 (MQ=255)

                     aTTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTa                                                                   >  2:1713688/1‑93 (MQ=255)

                     aTTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTa                                                                   <  2:3372492/93‑1 (MQ=255)

                      tttAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTac                                                                  >  2:88088/1‑93 (MQ=255)

                      tttAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTac                                                                  >  2:965384/1‑93 (MQ=255)

                      tttAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTac                                                                  >  2:4875353/1‑93 (MQ=255)

                      tttAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTac                                                                  <  1:2512433/93‑1 (MQ=255)

                      ttaaGTGATGCTCGCCGCGGCCACGGCAGATTTGATTGCCCGTGGTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTac                                                                  <  1:1252393/90‑1 (MQ=255)

                       ttAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTaca                                                                 <  1:221663/93‑1 (MQ=255)

                       ttAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTaca                                                                 <  2:3884270/93‑1 (MQ=255)

                          gTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCt                                                              >  2:5260401/1‑93 (MQ=255)

                            gATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTgc                                                            >  2:5617136/1‑93 (MQ=255)

                            gATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTgc                                                            <  2:4091412/93‑1 (MQ=255)

                             aTTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTgcg                                                           <  1:4274669/93‑1 (MQ=255)

                              ttCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTgcgc                                                          >  1:3008493/1‑93 (MQ=255)

                                 tcGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTg                                                       >  1:2992771/1‑93 (MQ=255)

                                   gCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTa                                                     >  1:801456/1‑93 (MQ=255)

                                   gCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTa                                                     >  2:4661899/1‑93 (MQ=255)

                                   gCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTa                                                     >  1:786845/1‑93 (MQ=255)

                                   gCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTa                                                     <  2:3240224/93‑1 (MQ=255)

                                    ccccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAg                                                    >  1:3390291/1‑93 (MQ=255)

                                     cccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGc                                                   <  1:1187651/93‑1 (MQ=255)

                                     cccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGc                                                   <  2:2096809/93‑1 (MQ=255)

                                     cccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGc                                                   <  2:4575828/93‑1 (MQ=255)

                                      ccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGcc                                                  <  1:4641225/93‑1 (MQ=255)

                                      ccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGcc                                                  >  2:4009818/1‑93 (MQ=255)

                                      ccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGcc                                                  >  2:5838394/1‑93 (MQ=255)

                                      ccTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGcc                                                  >  1:3351264/1‑93 (MQ=255)

                                         gCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTg                                               >  2:2213371/1‑93 (MQ=255)

                                          ccACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGc                                              <  2:1157933/93‑1 (MQ=255)

                                            aCGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTc                                            >  2:1434341/1‑93 (MQ=255)

                                              ggCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTccc                                          <  1:912479/93‑1 (MQ=255)

                                              ggCAGATTTGATTGCCCGTGTGGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTccc                                          <  1:1000737/93‑1 (MQ=255)

                                               gCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTcccc                                         <  2:1577155/93‑1 (MQ=255)

                                               gCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTcccc                                         >  1:5696530/1‑93 (MQ=255)

                                               gCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTcccc                                         >  2:2328165/1‑93 (MQ=255)

                                                cAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCg                                        >  2:197165/1‑93 (MQ=255)

                                                  gATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGcc                                      <  2:3031432/93‑1 (MQ=255)

                                                   aTTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCa                                     >  2:1497425/1‑93 (MQ=255)

                                                    tttGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCAt                                    <  1:2944410/93‑1 (MQ=255)

                                                     ttGATTGCCCGTATTGTTGCCGGAATGGCGATTGCCCCGTTCTCGAAGATTTGTCTGGCTACACCTGCGCCTGAAGCCGTGCTCCCCGCCATg                                   <  2:3615443/93‑1 (MQ=255)

                                                       gATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGaa                                 <  2:469745/93‑1 (MQ=255)

                                                         ttGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAAcc                               <  2:669092/93‑1 (MQ=255)

                                                          tGCCCGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACca                              <  2:3680731/93‑1 (MQ=255)

                                                            cccGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACcagc                            >  2:3510370/1‑93 (MQ=255)

                                                              cGTGTTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACcagcag                          >  1:1230381/1‑93 (MQ=255)

                                                               gtgtTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGa                         >  1:1581193/1‑93 (MQ=255)

                                                               gtgtGGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGa                         <  1:521107/93‑1 (MQ=255)

                                                                tgtTGCTGCCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGAt                        >  1:6254464/1‑93 (MQ=255)

                                                                    gctgcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTAc                    <  2:5108756/93‑1 (MQ=255)

                                                                    gctgcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTAc                    <  1:2374694/93‑1 (MQ=255)

                                                                      tgcCGGAATGGCGAATGACCGGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTAcct                  <  2:4185034/93‑2 (MQ=255)

                                                                      tgcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCg                  <  1:6170360/93‑1 (MQ=255)

                                                                      tgcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCg                  >  1:5969514/1‑93 (MQ=255)

                                                                       gcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGt                 <  2:611829/93‑1 (MQ=255)

                                                                       gcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGt                 <  1:1521637/93‑1 (MQ=255)

                                                                       gcCGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGt                 >  1:243925/1‑93 (MQ=255)

                                                                        ccGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTg                >  1:5004275/1‑93 (MQ=255)

                                                                        ccGCAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTg                <  2:2232525/93‑1 (MQ=255)

                                                                         cGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGc               <  1:4231741/93‑1 (MQ=255)

                                                                         cGGAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGc               >  1:498254/1‑93 (MQ=255)

                                                                          ggAATGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGcc              >  2:6036999/1‑93 (MQ=255)

                                                                            aaTGGCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGc            <  1:4458154/93‑1 (MQ=255)

                                                                               ggCGAATGACCCGGTATCGCCGATTTGTCTGGCTACCCCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGc         <  1:168574/93‑1 (MQ=255)

                                                                               ggCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGc         >  1:1786924/1‑93 (MQ=255)

                                                                               ggCGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGc         <  1:3891610/93‑1 (MQ=255)

                                                                                 ctaaTGACCCTGTATCGACGATTTGCCTGGCTACCCCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCa       <  1:4765276/91‑1 (MQ=255)

                                                                                 cGAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCa       <  2:4591589/93‑1 (MQ=255)

                                                                                  gAATGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGAAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCAc      >  1:720875/1‑93 (MQ=255)

                                                                                    aTGACCCGGTATCGACGATTTGTCTGGCTCCACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgc    <  1:5831194/93‑1 (MQ=255)

                                                                                    aTGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgc    <  1:257994/93‑1 (MQ=255)

                                                                                    aTGACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgc    >  2:837679/1‑93 (MQ=255)

                                                                                     ggACCCGGTATCGACGATTTGTCTGCCTACACCTGCGCCTGTAGCCGTCCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgca   <  1:3670829/92‑1 (MQ=255)

                                                                                      gACCCGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACgcag  >  1:5540279/1‑93 (MQ=255)

TTGAGCTGGGTAAATGGGCTGATTTAGTGATTCTCGCCCCTGCCACGGCAGATTTGATTGCCCGTGTTGCTGCCGGAATGGCGAATGACCTGGTATCGACGATTTGTCTGGCTACACCTGCGCCTGTAGCCGTGCTCCCCGCCATGAACCAGCAGATGTACCGTGCCGCTGCCACGCAG  >  NC_000913/3810584‑3810762

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