BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,566,875	A→G	R29R (CGT→CGC)	glgC ←	glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,875	0	A	G	100.0%	704.0 / NA	220	R29R (CGT→CGC)	glgC	glucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (152/68); total (152/68)

Only 100 of 284 total aligned reads displayed.

AGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566785‑3566964

 gACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGt                                                                                        >  2:2713388/1‑93 (MQ=255)

 gACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGt                                                                                        >  1:1331808/1‑93 (MQ=255)

 gACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGt                                                                                        >  2:3962648/1‑93 (MQ=255)

  aCAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtc                                                                                       >  1:4023423/1‑93 (MQ=255)

  aCAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtc                                                                                       <  1:4349547/93‑1 (MQ=255)

   cAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtcc                                                                                      <  2:2803396/93‑1 (MQ=255)

   cAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtcc                                                                                      <  1:5135464/93‑1 (MQ=255)

    aGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtcct                                                                                     <  2:2238027/93‑1 (MQ=255)

     gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtcctc                                                                                    >  2:1526054/1‑93 (MQ=255)

     gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtcctc                                                                                    <  2:4758572/93‑1 (MQ=255)

     gcgcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGtcctc                                                                                    >  1:1038594/1‑93 (MQ=255)

       gcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTccc                                                                                  <  2:6074475/93‑1 (MQ=255)

       gcAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTccc                                                                                  <  2:6345214/93‑1 (MQ=255)

          aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               >  1:2002147/1‑93 (MQ=255)

          aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               >  2:1017744/1‑93 (MQ=255)

          aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               <  2:4279998/93‑1 (MQ=255)

          aaGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGcc                                                                               <  2:3414341/93‑1 (MQ=255)

           aGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCa                                                                              >  1:2843107/1‑93 (MQ=255)

             gcGATAATGCGGAACTTCCCGCGGAAGTGTACGGCCGGTTTTGCTCGCTTATTGTTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGt                                                                            <  2:4000014/92‑1 (MQ=255)

              cGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTa                                                                           >  1:5200713/1‑93 (MQ=255)

              cGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTa                                                                           <  2:2489208/93‑1 (MQ=255)

                aTAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATc                                                                         >  1:3274261/1‑93 (MQ=255)

                  aaTGCGGAACTTCCCGCCGAAGTGTACGGCCGTTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGCCCTCCCGCCAGTATCAg                                                                       <  1:960215/93‑1 (MQ=255)

                  aaTGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAg                                                                       >  1:3127382/1‑93 (MQ=255)

                     gCGGAACTTCCCGCCGAAGTGTGCGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGc                                                                    <  1:5404270/93‑1 (MQ=255)

                       ggAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCaa                                                                  >  1:2456949/1‑93 (MQ=255)

                         aaCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACa                                                                >  2:6272513/1‑93 (MQ=255)

                          aCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               >  2:3841613/1‑93 (MQ=255)

                          aCTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAg                                                               >  1:6266969/1‑93 (MQ=255)

                            ttcccgccgAATTGTACGGCGGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAt                                                             <  1:383035/90‑1 (MQ=255)

                              tccgccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAttt                                                           <  2:4503581/92‑1 (MQ=255)

                              aCCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGAttt                                                           >  1:2358867/1‑93 (MQ=255)

                               ccgcggAAGTGTACGGCGGTTTTTGCTCTTTTATTGGTAAAC‑TCCTCCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTc                                                          <  2:5068931/93‑1 (MQ=255)

                               ccgccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTc                                                          >  1:1400473/1‑93 (MQ=255)

                               ccgccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTc                                                          >  1:6649543/1‑93 (MQ=255)

                                  ccgAAGTGTACGGCCGTTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       <  1:3698260/93‑1 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCTCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       <  2:3213756/93‑1 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       >  2:5547970/1‑93 (MQ=255)

                                  ccgAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       >  1:5762259/1‑93 (MQ=255)

                                  ccgAAGTGTACGCCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAAt                                                       <  1:5361544/93‑1 (MQ=255)

                                    gAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATgg                                                     >  1:1074986/1‑93 (MQ=255)

                                     aaGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    >  1:801837/1‑93 (MQ=255)

                                     aaGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGc                                                    >  2:2850778/1‑93 (MQ=255)

                                        tgtACGGCGGGTTTTGCTCGCTTATGGGTTAAA‑TCCTTCAGGCGGGTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGc                                                 <  1:172156/93‑1 (MQ=255)

                                        tgtACGGCCTGTTTTGCTCGCTTATGGGTTAAA‑TCCTGCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGc                                                 <  2:1440795/93‑1 (MQ=255)

                                        tgtACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGc                                                 <  2:5395690/93‑1 (MQ=255)

                                        tgtACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGc                                                 <  2:6688467/93‑1 (MQ=255)

                                         gtgcGGCCGTTTTTGCTCGCTTTTTGTTTAAA‑CCCTGCGGGGGGTTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCt                                                <  1:1760169/90‑1 (MQ=255)

                                          tACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑CCCTTCAGGCGGGTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTg                                               <  1:2720274/93‑1 (MQ=255)

                                             ggCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcg                                            >  2:2442759/1‑93 (MQ=255)

                                                cGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcgcgc                                         >  2:2027707/1‑93 (MQ=255)

                                                cGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGgcgcgc                                         >  2:1224803/1‑93 (MQ=255)

                                                 gtttttGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGcc                                        <  1:6430650/91‑1 (MQ=255)

                                                  gTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCa                                       >  1:4582704/1‑93 (MQ=255)

                                                   ttttGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCaa                                      >  1:4296443/1‑93 (MQ=255)

                                                     ttGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACa                                    <  2:817472/93‑1 (MQ=255)

                                                       gCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTTGCGCGCCAACAtt                                  >  2:3480049/1‑93 (MQ=255)

                                                        cTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTa                                 <  2:596658/93‑1 (MQ=255)

                                                        cTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTa                                 >  1:2982498/1‑93 (MQ=255)

                                                         tCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                >  2:5741452/1‑93 (MQ=255)

                                                         tCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                >  1:2538188/1‑93 (MQ=255)

                                                         tCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTaa                                >  2:3065545/1‑93 (MQ=255)

                                                          cGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAg                               >  2:6626356/1‑93 (MQ=255)

                                                          cGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAg                               >  2:2923545/1‑93 (MQ=255)

                                                           gCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              >  1:3802928/1‑93 (MQ=255)

                                                           gCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              >  1:899914/1‑93 (MQ=255)

                                                           gCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              <  2:2202683/93‑1 (MQ=255)

                                                           gCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              >  2:4616447/1‑93 (MQ=255)

                                                           gCTTATTGGTTAAA‑TCCTTCAGGCGGGGACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGt                              >  2:420708/1‑93 (MQ=255)

                                                            ttttttGGTTAAA‑TCCTTCAGGGGGGTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTg                             <  1:4822417/89‑1 (MQ=255)

                                                              tATTGTTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           <  2:679539/93‑1 (MQ=255)

                                                              tATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           >  1:1844744/1‑93 (MQ=255)

                                                              tATTGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           >  2:10791/1‑93 (MQ=255)

                                                              tATTGGTAAAA‑TCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGAt                           <  1:5892022/93‑1 (MQ=255)

                                                                ttGGTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCg                         >  2:1831065/1‑93 (MQ=255)

                                                                ttGGTAAAA‑TCTTTCGGGGGGATACGGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCg                         <  1:1558224/93‑1 (MQ=255)

                                                                  ggTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGGGCGCCAACATTAAGTGATCGtt                       >  1:3516145/1‑93 (MQ=255)

                                                                   gtaaaa‑TCCTTCAGGCGGGTCCCGCGCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTCATCGttc                      <  1:3641585/90‑1 (MQ=255)

                                                                   gTTAAATTCCTTCGGGCGGGTCCCGCCCCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGtt                       <  1:269828/93‑1 (MQ=255)

                                                                   gTTAAA‑TCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttc                      <  2:5861755/93‑1 (MQ=255)

                                                                   gTTAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGCATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttc                      >  1:2045471/1‑93 (MQ=255)

                                                                    ttAAA‑TCCTTCAGGCGGTTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttct                     <  2:4956181/93‑1 (MQ=255)

                                                                    ttAAA‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttct                     <  1:1128898/93‑1 (MQ=255)

                                                                     tAAA‑GCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttctt                    <  2:126509/93‑1 (MQ=255)

                                                                      aaa‑TCTTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttc                   <  1:4649685/93‑1 (MQ=255)

                                                                      aaa‑GCCTTCAGGCGGGTCCCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGttcttc                   <  2:939450/93‑1 (MQ=255)

                                                                        a‑TCCTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTtctc                 >  2:5187565/1‑93 (MQ=255)

                                                                          tctttCAGGCGGGTCCCGCGTCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTctct                <  1:414922/90‑1 (MQ=255)

                                                                            cTTCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaa              >  1:5844376/1‑93 (MQ=255)

                                                                             ttCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTaaa             >  2:4044699/1‑93 (MQ=255)

                                                                             ttCAGACGGGCAACGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTAGAGCGCCAACATTAATTGGTCGTTCTTATCTcaa             >  1:3942385/1‑93 (MQ=255)

                                                                              tCAGGCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAAc            >  1:2253537/1‑93 (MQ=255)

                                                                                 ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  2:1684650/1‑93 (MQ=255)

                                                                                 ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  2:1319546/1‑93 (MQ=255)

                                                                                 ggCGGGTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTaa         >  1:4994253/1‑93 (MQ=255)

                                                                                 ggCGGGTACCGCCTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGGTCTTCTCTAAACTaa         >  1:5684077/1‑93 (MQ=255)

                                                                                     ggTACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     >  1:2585206/1‑93 (MQ=255)

                                                                                     ggTACCGCGCCCCCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCAt     <  2:533581/93‑1 (MQ=255)

                                                                                      gTACCGCGTCGTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATg    >  2:6384727/1‑93 (MQ=255)

                                                                                       tACCGCGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGa   <  2:3185181/93‑1 (MQ=255)

AGACAGCGCAAAGTCGATAATGCGGAACTTACCGCCGAAGTGTACGGCCGGTTTTGCTCGCTTATTGGTTAAA‑TCCTTCAGGCGGGTACCACGTCCTCCCGCCAGTATCAGGGCAACAGATTTCAATGGCAGCTGGCGCGCCAACATTAAGTGATCGTTCTTCTCTAAACTAACCATGAC  >  NC_000913/3566785‑3566964

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