BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	350,318	A→G	Q361Q (CAA→CAG)	prpC →	putative citrate synthase; propionate metabolism?

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	350,318	0	A	G	100.0%	492.3 / NA	157	Q361Q (CAA→CAG)	prpC	putative citrate synthase; propionate metabolism?
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (57/100); total (57/100)

Only 100 of 164 total aligned reads displayed.

CTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAA  >  NC_000913/350228‑350406

cTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCg                                                                                        <  1:114218/93‑1 (MQ=255)

  acaacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTc                                                                                      <  1:4116340/93‑1 (MQ=255)

  acaacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTc                                                                                      <  1:258462/93‑1 (MQ=255)

  acaacaTGACGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTc                                                                                      >  1:492531/1‑93 (MQ=255)

    aacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAg                                                                                    <  2:3229403/93‑1 (MQ=255)

    aacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAg                                                                                    <  1:2924334/93‑1 (MQ=255)

    aacaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAg                                                                                    <  1:4655283/93‑1 (MQ=255)

     acaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGGGCACATTATCGAACCGCGTCAgg                                                                                   >  2:269002/1‑93 (MQ=255)

     acaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAgg                                                                                   >  1:1522688/1‑93 (MQ=255)

     acaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAgg                                                                                   <  2:2032367/93‑1 (MQ=255)

     acaTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAgg                                                                                   <  2:2561644/93‑1 (MQ=255)

      caTGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGa                                                                                  <  1:2353017/93‑1 (MQ=255)

       atgatgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGac                                                                                 >  1:1439428/1‑93 (MQ=255)

        tgatgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGaca                                                                                >  2:1961847/1‑93 (MQ=255)

         gatgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGacaa                                                                               <  2:1563854/93‑1 (MQ=255)

          atgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGacaac                                                                              <  2:4236652/93‑1 (MQ=255)

          atgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGacaac                                                                              >  1:327552/1‑93 (MQ=255)

           tgGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGacaaca                                                                             <  2:3320198/93‑1 (MQ=255)

               cGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAAt                                                                         <  1:1068412/93‑1 (MQ=255)

               cGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAAt                                                                         <  1:569512/93‑1 (MQ=255)

                gTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAAtt                                                                        <  2:838038/93‑1 (MQ=255)

                 ttCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTa                                                                       <  1:1398673/93‑1 (MQ=255)

                 ttCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTa                                                                       <  1:120911/93‑1 (MQ=255)

                   cccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATc                                                                     <  2:3485889/93‑1 (MQ=255)

                   cccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATc                                                                     >  2:2099869/1‑93 (MQ=255)

                    ccACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGGGTCAAGACCACAAAATTATcc                                                                    >  2:295953/1‑93 (MQ=255)

                     cACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCg                                                                   <  1:2459532/93‑1 (MQ=255)

                      aCCGAGATGTTCACACCACTGGTTGTTATTGCCCGCGTCACCGGCTGGGCGGCGCACATTACCGAACACCGTCCGGTACACAAAATTATTCGt                                                                  >  2:1496097/1‑93 (MQ=255)

                       ccGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTc                                                                 <  1:2741184/93‑1 (MQ=255)

                       ccGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTc                                                                 >  2:3601200/1‑93 (MQ=255)

                        cGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTcc                                                                <  1:4659292/93‑1 (MQ=255)

                         gagaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCt                                                               >  1:1801585/1‑93 (MQ=255)

                         gagaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCt                                                               >  1:3526414/1‑93 (MQ=255)

                          agaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCtt                                                              >  2:4363868/1‑93 (MQ=255)

                           gaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTc                                                             <  2:2554832/93‑1 (MQ=255)

                           gaTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTc                                                             >  2:2609561/1‑93 (MQ=255)

                            aTGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTcc                                                            <  2:1115843/93‑1 (MQ=255)

                             tGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCg                                                           >  1:3737486/1‑93 (MQ=255)

                               ttCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGcc                                                         >  1:2074751/1‑93 (MQ=255)

                                 cacaCCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCaa                                                       >  2:583986/1‑93 (MQ=255)

                                  acacCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAAt                                                      <  1:2848347/93‑1 (MQ=255)

                                    accacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTa                                                    <  2:1567370/93‑1 (MQ=255)

                                     ccacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCATTCCGACAATTAt                                                   >  2:2143289/1‑93 (MQ=255)

                                      ctcTGTTTGTTATCGCCCACGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATg                                                  <  1:2359186/91‑1 (MQ=255)

                                      cacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATg                                                  <  1:3223216/93‑1 (MQ=255)

                                      cacTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATg                                                  >  2:953227/1‑93 (MQ=255)

                                       acTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGt                                                 <  1:2004891/93‑1 (MQ=255)

                                        cTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGtt                                                <  1:4396565/93‑1 (MQ=255)

                                            ttgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGAc                                            >  2:3333680/1‑93 (MQ=255)

                                             tgttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGAcc                                           >  1:2303469/1‑93 (MQ=255)

                                              gttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCg                                          >  2:4204567/1‑93 (MQ=255)

                                               ttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCgg                                         >  1:29674/1‑93 (MQ=255)

                                               ttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCgg                                         >  2:1309662/1‑93 (MQ=255)

                                               ttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCgg                                         <  1:2198837/93‑1 (MQ=255)

                                               ttATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCgg                                         >  2:628299/1‑93 (MQ=255)

                                                 aTCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGaa                                       <  2:2849960/93‑1 (MQ=255)

                                                 aTCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGaa                                       >  2:1498998/1‑93 (MQ=255)

                                                   cGCCCGCGTCCCCGGCTGGGCGGCGCCCATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGa                                     <  1:2820294/93‑1 (MQ=255)

                                                   cGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACGACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGa                                     >  1:393605/1‑93 (MQ=255)

                                                   cGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGa                                     <  1:1834020/93‑1 (MQ=255)

                                                   cGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGa                                     <  1:2374815/93‑1 (MQ=255)

                                                    gCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAc                                    <  1:2155095/93‑1 (MQ=255)

                                                     cccGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAcc                                   >  2:4648677/1‑93 (MQ=255)

                                                     cccGCGTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGAcc                                   >  1:3013129/1‑93 (MQ=255)

                                                        gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                <  1:1833404/93‑1 (MQ=255)

                                                        gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                >  2:236599/1‑93 (MQ=255)

                                                        gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                <  1:2934770/93‑1 (MQ=255)

                                                        gcgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGcc                                <  1:1554934/93‑1 (MQ=255)

                                                         cgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGccc                               >  2:413244/1‑93 (MQ=255)

                                                         cgTCACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGccc                               >  2:843371/1‑93 (MQ=255)

                                                            cACCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGtt                            <  2:2031927/93‑1 (MQ=255)

                                                             aCCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCTTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGttt                           >  2:4487141/1‑93 (MQ=255)

                                                             aCCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGttt                           >  2:21922/1‑93 (MQ=255)

                                                             aCCGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGttt                           >  1:536174/1‑93 (MQ=255)

                                                              ccGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTg                          <  1:4472957/93‑1 (MQ=255)

                                                              ccGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTg                          <  1:3160335/93‑1 (MQ=255)

                                                               cGGCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGt                         <  1:2291029/93‑1 (MQ=255)

                                                                ggCTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTc                        <  1:970398/93‑1 (MQ=255)

                                                                  cTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcgc                      <  1:2982665/93‑1 (MQ=255)

                                                                  cTGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcgc                      <  2:2936435/93‑1 (MQ=255)

                                                                   tGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcgcg                     <  1:309561/93‑1 (MQ=255)

                                                                   tGGGCGGCGCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTcgcg                     >  2:4544204/1‑93 (MQ=255)

                                                                       cggcgCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTgg                 <  1:3936924/93‑1 (MQ=255)

                                                                        ggcgCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGa                <  1:896194/93‑1 (MQ=255)

                                                                        ggcgCACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTGCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGa                <  2:369361/93‑1 (MQ=255)

                                                                           gcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATaa             <  1:3287983/93‑1 (MQ=255)

                                                                           gcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATaa             >  2:3926316/1‑93 (MQ=255)

                                                                           gcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATaa             <  1:1125080/93‑1 (MQ=255)

                                                                           gcACATTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATaa             <  1:2412594/93‑1 (MQ=255)

                                                                            cacaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAg            <  1:124193/93‑1 (MQ=255)

                                                                             acaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgc           <  1:1080087/93‑1 (MQ=255)

                                                                             acaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgc           <  1:1586807/93‑1 (MQ=255)

                                                                             acaTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgc           >  1:428323/1‑93 (MQ=255)

                                                                              caTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgcg          >  1:3992846/1‑93 (MQ=255)

                                                                               aTTATCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAgcgc         <  1:3209039/93‑1 (MQ=255)

                                                                               aTTATCGAACAGCGTCAGGACAACAAAATGATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGGCCGCCCGTTTGTCGCGCTGGATAAgcgc         <  2:4433095/93‑1 (MQ=255)

                                                                                 tagcGAACAGCGTCAGGAGAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCa       <  1:4179016/90‑1 (MQ=255)

                                                                                  aTCGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAg      <  2:4061235/93‑1 (MQ=255)

                                                                                    cGAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTa    <  1:1688047/93‑1 (MQ=255)

                                                                                     gAACAGCGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTaa   <  1:360021/93‑1 (MQ=255)

CTACAACATGATGGGCGTTCCCACCGAGATGTTCACACCACTGTTTGTTATCGCCCGCGTCACCGGCTGGGCGGCGCACATTATCGAACAACGTCAGGACAACAAAATTATCCGTCCTTCCGCCAATTATGTTGGACCGGAAGACCGCCCGTTTGTCGCGCTGGATAAGCGCCAGTAAA  >  NC_000913/350228‑350406

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