BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,442,757	C→T	D561D (GAC→GAT)	ydbH →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,442,757	0	C	T	83.6%	671.3 / 91.5	262	D561D (GAC→GAT)	ydbH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (22/21); new base T (129/90); total (151/111)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 263 total aligned reads displayed.

CGTCGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGAACTTACG  >  NC_000913/1442667‑1442849

 gTCGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTTCATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGt                                                                                           >  2:4751910/1‑93 (MQ=255)

 gTCGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGt                                                                                           >  2:1846370/1‑93 (MQ=255)

  tCGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTa                                                                                          >  1:672738/1‑93 (MQ=255)

   cGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTAt                                                                                         >  2:4271966/1‑93 (MQ=255)

   cGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTAt                                                                                         >  1:3609024/1‑93 (MQ=255)

   cGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTAt                                                                                         <  2:6781812/93‑1 (MQ=255)

    gATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTAtt                                                                                        >  1:2469806/1‑93 (MQ=255)

      tGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCg                                                                                      <  1:1618211/93‑1 (MQ=255)

        ggCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTc                                                                                    >  1:3124324/1‑93 (MQ=255)

        ggCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTc                                                                                    <  1:1023326/93‑1 (MQ=255)

           gcgATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTgc                                                                                 <  1:1849668/93‑1 (MQ=255)

           gcgATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTgc                                                                                 <  1:1303058/93‑1 (MQ=255)

           gcgATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTgc                                                                                 <  1:3584468/93‑1 (MQ=255)

             gATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTgcgc                                                                               <  2:2106733/93‑1 (MQ=255)

              aTCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGcgcg                                                                              >  2:5827527/1‑93 (MQ=255)

               tCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCgg                                                                             <  2:2824888/93‑1 (MQ=255)

               tCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCgg                                                                             >  2:3384872/1‑93 (MQ=255)

               tCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCgg                                                                             <  2:4186770/93‑1 (MQ=255)

                 ctacctaTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTa                                                                           >  1:979387/1‑93 (MQ=255)

                  tacctaTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTaa                                                                          <  1:6319714/93‑1 (MQ=255)

                   acctaTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAAc                                                                         <  1:1608524/93‑1 (MQ=255)

                     ctaTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGGGCGGTAACGc                                                                       <  1:6080051/93‑1 (MQ=255)

                      taTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGcc                                                                      >  2:6629261/1‑93 (MQ=255)

                      taTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGcc                                                                      <  2:447141/93‑1 (MQ=255)

                       aTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCt                                                                     >  1:423132/1‑93 (MQ=255)

                       aTTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCt                                                                     >  2:1140271/1‑93 (MQ=255)

                        tttCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCtg                                                                    >  1:1361258/1‑93 (MQ=255)

                          tCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCtggt                                                                  >  1:5025312/1‑93 (MQ=255)

                          tCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCtggt                                                                  <  2:5243800/93‑1 (MQ=255)

                           ccTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCtggtg                                                                 >  1:3368733/1‑93 (MQ=255)

                            ctctTTAAAGGCGATTTACATGCTGGTGTGTTTGGCCCGGTTCGGGAAAATGGTCGCTGGGATGGCATTCGTCTGCGCGGTAACGCCtggtgg                                                                >  1:3868644/1‑93 (MQ=255)

                             tctTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGc                                                               >  1:4473874/1‑93 (MQ=255)

                             tctTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGc                                                               >  2:2604128/1‑93 (MQ=255)

                              ctTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGcc                                                              >  2:736106/1‑93 (MQ=255)

                               tttAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCt                                                             >  1:5464265/1‑93 (MQ=255)

                               tttAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCt                                                             <  2:6322365/93‑1 (MQ=255)

                               tttAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCt                                                             >  2:5297949/1‑93 (MQ=255)

                                ttAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTa                                                            >  2:1762963/1‑93 (MQ=255)

                                  aaaGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTaaa                                                          >  2:2413809/1‑93 (MQ=255)

                                   aaGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAAc                                                         >  2:1539541/1‑93 (MQ=255)

                                   aaGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAAc                                                         >  2:3554272/1‑93 (MQ=255)

                                        gATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCa                                                    <  2:5752474/93‑1 (MQ=255)

                                         aTTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCAc                                                   <  2:2244140/93‑1 (MQ=255)

                                         aTTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCAc                                                   >  2:3327733/1‑93 (MQ=255)

                                         aTTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCAc                                                   >  2:4300082/1‑93 (MQ=255)

                                         aTTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCAc                                                   >  2:2014103/1‑93 (MQ=255)

                                             aCATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGAc                                               <  1:4051141/93‑1 (MQ=255)

                                              cATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGAcc                                              >  1:6962541/1‑93 (MQ=255)

                                              cATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGAcc                                              >  1:7212240/1‑93 (MQ=255)

                                               aTGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCg                                             >  2:4841308/1‑93 (MQ=255)

                                                tGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGt                                            >  1:4403901/1‑93 (MQ=255)

                                                 gCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTa                                           <  2:3670506/93‑1 (MQ=255)

                                                 gCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTa                                           <  1:1975921/93‑1 (MQ=255)

                                                  cTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTAt                                          >  1:1140363/1‑93 (MQ=255)

                                                   tGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTAtt                                         >  2:3034513/1‑93 (MQ=255)

                                                   tGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTAtt                                         >  1:409096/1‑93 (MQ=255)

                                                     gTGAGATTGGCCCGGTTAGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCATTGACCGTATTcc                                       <  1:538191/93‑1 (MQ=255)

                                                       gagaTTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAg                                     >  1:4699475/1‑93 (MQ=255)

                                                        agaTTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGc                                    <  2:878028/93‑1 (MQ=255)

                                                        agaTTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGc                                    <  1:526490/93‑1 (MQ=255)

                                                        agaTTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGc                                    <  1:7520727/93‑1 (MQ=255)

                                                        agaTTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTACTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGc                                    >  2:2391442/1‑93 (MQ=255)

                                                         gaTTGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGcc                                   >  1:3672957/1‑93 (MQ=255)

                                                           ttGGCCCGGTTCGGGTAAGTGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGc                                 <  2:5995639/93‑1 (MQ=255)

                                                           ttGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGc                                 >  2:3511550/1‑93 (MQ=255)

                                                           ttGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGc                                 >  2:3886231/1‑93 (MQ=255)

                                                           ttGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGc                                 >  2:7169263/1‑93 (MQ=255)

                                                            tGGCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCt                                <  2:6698332/93‑1 (MQ=255)

                                                              gCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTgg                              <  1:2138236/93‑1 (MQ=255)

                                                              gCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTgg                              >  2:1513/1‑93 (MQ=255)

                                                              gCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTgg                              <  2:3930232/93‑1 (MQ=255)

                                                              gCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTgg                              <  2:461819/93‑1 (MQ=255)

                                                              gCCCGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTgg                              >  1:5131415/1‑93 (MQ=255)

                                                              gCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTgg                              >  2:1234832/1‑93 (MQ=255)

                                                               cccGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGt                             <  1:6379187/93‑1 (MQ=255)

                                                                 cGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGc                           <  2:1876523/93‑1 (MQ=255)

                                                                 cGGTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGc                           <  2:6268414/93‑1 (MQ=255)

                                                                 cGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGc                           <  1:6054960/93‑1 (MQ=255)

                                                                   gTTCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCa                         <  2:3114697/93‑1 (MQ=255)

                                                                   gTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCa                         <  1:887500/93‑1 (MQ=255)

                                                                    ttCGGGTAAATGGTTGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCAc                        <  2:2937112/93‑1 (MQ=255)

                                                                    ttCGGGTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCAc                        >  2:5156294/1‑93 (MQ=255)

                                                                       gggTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCg                     >  2:7046556/1‑93 (MQ=255)

                                                                       gggTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCg                     >  1:269507/1‑93 (MQ=255)

                                                                        ggTAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGa                    >  1:3016562/1‑93 (MQ=255)

                                                                          tAAATGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACt                  <  2:3682488/93‑1 (MQ=255)

                                                                          tAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACt                  <  1:3935002/93‑1 (MQ=255)

                                                                            aaTGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTgg                >  1:3199544/1‑93 (MQ=255)

                                                                            aaTGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTgg                <  2:615996/93‑1 (MQ=255)

                                                                             aTGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGa               >  2:5095796/1‑93 (MQ=255)

                                                                              tGGTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGaa              <  1:4876217/93‑1 (MQ=255)

                                                                                gTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGa            <  2:7401186/93‑1 (MQ=255)

                                                                                gTCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACAGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGa            <  1:665242/93‑1 (MQ=255)

                                                                                 tCGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGAt           >  1:5860584/1‑93 (MQ=255)

                                                                                  cGCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATg          <  2:1306710/93‑1 (MQ=255)

                                                                                   gCTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGa         <  1:4303988/93‑1 (MQ=255)

                                                                                    cTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGaa        >  2:5573784/1‑93 (MQ=255)

                                                                                    cTGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGaa        >  1:3957782/1‑93 (MQ=255)

                                                                                     tGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGAAc       >  1:2800124/1‑93 (MQ=255)

                                                                                     tGGGATGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGAAc       >  1:1447477/1‑93 (MQ=255)

CGTCGATGGGCGCGATCCTACCTATTTCCTCTTTAAAGGCGATTTACATGCTGGTGAGATTGGCCCGGTTCGGGTAAATGGTCGCTGGGACGGTATTCGTCTGCGCGGTAACGCCTGGTGGCCTAAACAATCACTGACCGTATTCCAGCCGCTGGTGCCACCCGACTGGAAGATGAACTTACG  >  NC_000913/1442667‑1442849

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