BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,971,075	2 bp→GC	coding (942‑943/1068 nt)	rffG →	dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,971,075	0	C	G	99.6%	903.7 / NA	251	C314W (TGC→TGG)	rffG	dTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base G (105/145); minor base T (0/1); total (105/146)
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
*	NC_000913	3,971,076	0	G	C	100.0%	906.8 / NA	251	V315L (GTG→CTG)	rffG	dTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (105/146); total (105/146)

Only 100 of 265 total aligned reads displayed.

TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162

||

 gACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGTCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCg                                                                                   <  1:3170413/93‑1 (MQ=255)

 gACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCg                                                                                   >  1:1874804/1‑93 (MQ=255)

   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 <  1:3425352/93‑1 (MQ=255)

   ccTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTCATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCa                                                                                 >  2:2394430/1‑93 (MQ=255)

    cTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAg                                                                                >  2:5617225/1‑93 (MQ=255)

     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  2:46002/1‑93 (MQ=255)

     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  2:5786274/1‑93 (MQ=255)

     tGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAgg                                                                               >  1:979374/1‑93 (MQ=255)

      gATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGa                                                                              <  1:196252/93‑1 (MQ=255)

        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGaaa                                                                            >  1:5590421/1‑93 (MQ=255)

        tCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAAATTGGCTGGCTGCCGCAGGaac                                                                            >  1:2448453/1‑92 (MQ=255)

         cACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAc                                                                           <  2:4155798/93‑1 (MQ=255)

          aCCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAAcc                                                                          <  2:3731343/93‑1 (MQ=255)

           ccTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCt                                                                         <  1:2073705/93‑1 (MQ=255)

            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        >  2:2028905/1‑93 (MQ=255)

            cTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCtt                                                                        <  2:3513308/93‑1 (MQ=255)

             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  1:2363050/93‑1 (MQ=255)

             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       >  1:4053925/1‑93 (MQ=255)

             tttGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCttt                                                                       <  2:4898749/93‑1 (MQ=255)

                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   <  1:5880976/93‑1 (MQ=255)

                 tCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGaaa                                                                   >  1:3378922/1‑93 (MQ=255)

                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 <  1:3974608/93‑1 (MQ=255)

                   gCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGt                                                                 >  1:1106048/1‑93 (MQ=255)

                    cTGACCGTCCGGGGCATGAACTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTg                                                                <  2:4339805/93‑1 (MQ=255)

                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTTCCCGTGAACTTGGCTGGCTGCCGGCGGAAACCTTTGAAAGTgg                                                               >  1:2707488/1‑93 (MQ=255)

                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTgg                                                               <  1:828064/93‑1 (MQ=255)

                     tGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTgg                                                               >  2:3435418/1‑93 (MQ=255)

                      gACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGa                                                              >  1:2013472/1‑93 (MQ=255)

                       gccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  2:1354213/92‑1 (MQ=255)

                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:618863/93‑1 (MQ=255)

                       aCCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGaa                                                             <  1:5469066/93‑1 (MQ=255)

                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            <  1:4174286/93‑1 (MQ=255)

                        ccGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAAt                                                            >  2:3857110/1‑93 (MQ=255)

                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGGATGc                                                          >  1:3598725/1‑93 (MQ=255)

                          gTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGc                                                          <  1:502578/93‑1 (MQ=255)

                           tCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCTGCAGGAAACCTTTGAAAGTGGAATGCg                                                         <  2:3921780/93‑1 (MQ=255)

                            ccGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGt                                                        >  2:1889921/1‑93 (MQ=255)

                            ccGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGt                                                        <  1:1196336/93‑1 (MQ=255)

                             cGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTa                                                       <  2:515545/93‑1 (MQ=255)

                             cGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTa                                                       >  1:1561804/1‑93 (MQ=255)

                               gggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaa                                                     <  1:3248113/93‑1 (MQ=255)

                                ggCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTaaaa                                                    <  1:4629007/93‑1 (MQ=255)

                                 gCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAAc                                                   <  2:2778644/93‑1 (MQ=255)

                                  cATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACg                                                  <  1:5648533/93‑1 (MQ=255)

                                  cATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGACGGAAACCTTTGAAAGTGGAATGCGCaaaagg                                                  >  2:1698602/1‑91 (MQ=255)

                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                <  2:3743925/93‑1 (MQ=255)

                                    tGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGt                                                >  2:2763517/1‑93 (MQ=255)

                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGtt                                               >  2:74270/1‑93 (MQ=255)

                                     gATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCGTTGAAAGTGGAATGCGGAAAACGGtt                                               >  1:5500573/1‑93 (MQ=255)

                                       tCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCa                                             <  2:47413/93‑1 (MQ=255)

                                        cTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAg                                            >  1:2829308/1‑93 (MQ=255)

                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAACCCTTTGAAAGTGCACTGCGTAAACAGGTTCCGt                                           >  1:2310738/1‑93 (MQ=255)

                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGt                                           >  1:4807399/1‑93 (MQ=255)

                                         tGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGt                                           <  2:4771448/93‑1 (MQ=255)

                                          gCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTg                                          <  2:1301283/93‑1 (MQ=255)

                                            gTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGt                                        <  1:904990/93‑1 (MQ=255)

                                             ttATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTa                                       <  2:6466/93‑1 (MQ=255)

                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      <  1:2678969/93‑1 (MQ=255)

                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      >  1:1479354/1‑93 (MQ=255)

                                              tATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTAt                                      <  1:5105840/93‑1 (MQ=255)

                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  1:4415156/93‑1 (MQ=255)

                                               aTGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATc                                     <  1:4640861/93‑1 (MQ=255)

                                                 gCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTg                                   >  2:5511050/1‑93 (MQ=255)

                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTTGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  2:1180638/1‑93 (MQ=255)

                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  2:2660842/1‑93 (MQ=255)

                                                   cATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGc                                 >  1:5830430/1‑93 (MQ=255)

                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  2:2252962/1‑93 (MQ=255)

                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  1:21725/1‑93 (MQ=255)

                                                    aTTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCt                                >  2:5216575/1‑93 (MQ=255)

                                                      tGATGCTTCGAAAATTGCCCGTTAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTaa                              <  1:1704138/93‑1 (MQ=255)

                                                       gATGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAAt                             >  2:2936814/1‑93 (MQ=255)

                                                        aTGCTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTAGCTCATg                            >  2:5403150/1‑93 (MQ=255)

                                                           cTTCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGaaa                         <  2:3000601/93‑1 (MQ=255)

                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        >  2:5351733/1‑93 (MQ=255)

                                                            ttCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAg                        <  2:2278622/93‑1 (MQ=255)

                                                             tCGAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGc                       <  1:257105/93‑1 (MQ=255)

                                                               gAAAATTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtg                     >  1:5744577/1‑93 (MQ=255)

                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    <  2:3717246/93‑1 (MQ=255)

                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    >  1:5210434/1‑93 (MQ=255)

                                                                aaaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtgg                    <  2:245017/93‑1 (MQ=255)

                                                                 aaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggt                   <  2:3193517/93‑1 (MQ=255)

                                                                 aaaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggt                   <  2:3599983/93‑1 (MQ=255)

                                                                  aaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtg                  <  2:2274888/93‑1 (MQ=255)

                                                                  aaTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtg                  >  1:6028191/1‑93 (MQ=255)

                                                                   aTTGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCtggtgg                 >  1:114175/1‑93 (MQ=255)

                                                                    ttGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGa                <  2:599936/93‑1 (MQ=255)

                                                                     tGCCCGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGaa               <  2:4614730/93‑1 (MQ=255)

                                                                      gCCCGTGAACTTGGCTGTCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAg              <  1:3734448/93‑1 (MQ=255)

                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCa            <  1:4230413/93‑1 (MQ=255)

                                                                        ccGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGGATCTGGCTAATGAAAGCTGGTGGAAGCa            <  1:2990837/93‑1 (MQ=255)

                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAg           <  1:3819375/93‑1 (MQ=255)

                                                                         cGTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAg           >  2:952654/1‑93 (MQ=255)

                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  1:5932383/93‑1 (MQ=255)

                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          >  1:4237006/1‑93 (MQ=255)

                                                                          gTGAACTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAgg          <  2:69928/93‑1 (MQ=255)

                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTga       >  2:4295349/1‑92 (MQ=255)

                                                                             aaCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGc       >  2:5256613/1‑93 (MQ=255)

                                                                              aCTTGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCa      <  1:1350089/93‑1 (MQ=255)

                                                                                ttGGCTGGCTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAgg    >  2:2882194/1‑93 (MQ=255)

                                                                                 tggctggcTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGa   <  1:1235283/93‑1 (MQ=255)

||

TGACCTGATCACCTTTGTCGCTGACCGTCCGGGGCATGATCTGCGTTATGCCATTGATGCTTCGAAAATTGCCCGTGAACTTGGCTGCGTGCCGCAGGAAACCTTTGAAAGTGGAATGCGTAAAACGGTTCAGTGGTATCTGGCTAATGAAAGCTGGTGGAAGCAGGTGCAGGAC  >  NC_000913/3970988‑3971162

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