BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,876,096	A→G	R78R (CGT→CGC)	cysD ←	sulfate adenylyltransferase subunit 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,876,096	0	A	G	100.0%	1323.0 / NA	317	R78R (CGT→CGC)	cysD	sulfate adenylyltransferase subunit 2
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (177/140); total (177/140)

Only 100 of 317 total aligned reads displayed.

GTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCA  >  NC_000913/2875998‑2876193

gTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCg                                                                                                 <  1:4679864/101‑1 (MQ=255)

 tATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGa                                                                                                >  1:4680428/1‑101 (MQ=255)

  aTGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGAt                                                                                               >  2:4680429/1‑101 (MQ=255)

   tGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATc                                                                                              <  1:4679939/101‑1 (MQ=255)

   tGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATc                                                                                              >  1:4680435/1‑101 (MQ=255)

    gTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATcg                                                                                             <  1:4680199/101‑1 (MQ=255)

      ttCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATcgcg                                                                                           >  2:4680438/1‑101 (MQ=255)

      ttCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATcgcg                                                                                           <  2:4680060/101‑1 (MQ=255)

       tCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCgg                                                                                          >  1:4680439/1‑101 (MQ=255)

       tCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCgg                                                                                          <  2:4680099/101‑1 (MQ=255)

          cgcTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAAc                                                                                       >  2:4680443/1‑101 (MQ=255)

          cgcTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAAc                                                                                       <  1:4680146/101‑1 (MQ=255)

           gcTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACt                                                                                      <  2:4679944/101‑1 (MQ=255)

            cTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTc                                                                                     >  1:4680446/1‑101 (MQ=255)

              gCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCAt                                                                                   >  2:4680450/1‑101 (MQ=255)

                cGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATAc                                                                                 >  1:4680453/1‑101 (MQ=255)

                 gTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACa                                                                                >  2:4680456/1‑101 (MQ=255)

                  tGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAt                                                                               >  1:4680461/1‑101 (MQ=255)

                  tGCACGAATGGATTAATCCCCATCGCCACGCCTCCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAt                                                                               >  2:4680460/1‑101 (MQ=255)

                   gCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAtc                                                                              >  1:4680467/1‑101 (MQ=255)

                    cACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAtct                                                                             >  2:4680470/1‑101 (MQ=255)

                    cACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAtct                                                                             <  1:4679650/101‑1 (MQ=255)

                     aCGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAtctc                                                                            >  1:4680473/1‑101 (MQ=255)

                      cGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTcg                                                                           <  1:4680071/101‑1 (MQ=255)

                      cGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTcg                                                                           <  2:4679924/101‑1 (MQ=255)

                       gAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTcgc                                                                          <  1:4680053/101‑1 (MQ=255)

                        aaTGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTcgcg                                                                         >  2:4680477/1‑101 (MQ=255)

                           ggATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGaa                                                                      <  2:4679963/101‑1 (MQ=255)

                           ggATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGaa                                                                      >  1:4680479/1‑101 (MQ=255)

                              ttAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAAttt                                                                   >  1:4680483/1‑101 (MQ=255)

                              ttAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAAttt                                                                   <  2:4679526/101‑1 (MQ=255)

                              ttAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAAttt                                                                   <  1:4680021/101‑1 (MQ=255)

                                   ccccATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGc                                                              >  1:4680490/1‑101 (MQ=255)

                                   ccccATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGc                                                              >  2:4680487/1‑101 (MQ=255)

                                    cccATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGcc                                                             <  1:4679954/101‑1 (MQ=255)

                                    cccATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGcc                                                             <  2:4680154/101‑1 (MQ=255)

                                    cccATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGcc                                                             >  2:4680493/1‑101 (MQ=255)

                                      cATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCgg                                                           >  1:4680497/1‑101 (MQ=255)

                                      cATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCgg                                                           <  2:4679908/101‑1 (MQ=255)

                                        ttgCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTa                                                         >  2:4680502/3‑101 (MQ=255)

                                        tCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTa                                                         >  1:4680503/1‑101 (MQ=255)

                                         cGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTAt                                                        <  2:4679773/101‑1 (MQ=255)

                                          gCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATc                                                       >  2:4680507/1‑101 (MQ=255)

                                          gCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATc                                                       <  1:4679982/101‑1 (MQ=255)

                                          gCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATc                                                       <  2:4680070/101‑1 (MQ=255)

                                           ccACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACAtct                                                                             <  1:4679819/78‑1 (MQ=255)

                                           ccACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCg                                                      <  2:4679852/101‑1 (MQ=255)

                                           ccACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCg                                                      >  1:4680512/1‑101 (MQ=255)

                                            cACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGa                                                     >  2:4680516/1‑101 (MQ=255)

                                            cACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGa                                                     >  1:4680515/1‑101 (MQ=255)

                                             aCGCCTTCCGGGTTTTTATTCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGAc                                                    <  1:4679972/101‑1 (MQ=255)

                                               gCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACAt                                                  >  2:4680520/1‑101 (MQ=255)

                                               gCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACAt                                                  <  2:4680006/101‑1 (MQ=255)

                                                ccTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATg                                                 >  1:4680523/1‑101 (MQ=255)

                                                 cTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATgc                                                <  2:4679813/101‑1 (MQ=255)

                                                  ttCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATgca                                               >  1:4680526/1‑101 (MQ=255)

                                                  ttCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATgca                                               <  1:4680208/101‑1 (MQ=255)

                                                    ccGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATgcagc                                             <  1:4680090/101‑1 (MQ=255)

                                                    ccGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATgcagc                                             <  1:4680389/101‑1 (MQ=255)

                                                     cGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATgcagca                                            >  2:4680529/1‑101 (MQ=255)

                                                       ggTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAAc                                          <  2:4679897/101‑1 (MQ=255)

                                                       ggTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAAc                                          >  1:4680535/1‑101 (MQ=255)

                                                          ttttATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACggg                                       >  1:4680538/1‑101 (MQ=255)

                                                           tttATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGa                                      >  2:4680541/1‑101 (MQ=255)

                                                            ttATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGaa                                     <  2:4680025/101‑1 (MQ=255)

                                                            ttATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGaa                                     >  1:4680543/1‑101 (MQ=255)

                                                                gCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGt                                                          <  2:4679984/76‑1 (MQ=255)

                                                                 cACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCa                                >  2:4680549/1‑101 (MQ=255)

                                                                 cACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCa                                >  1:4680548/1‑101 (MQ=255)

                                                                    cagcagTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCg                             <  1:4680172/101‑1 (MQ=255)

                                                                     agcagTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGt                            <  1:4680155/101‑1 (MQ=255)

                                                                      gcagTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTa                           >  1:4680553/1‑101 (MQ=255)

                                                                        agTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAg                               >  2:4680558/1‑95 (MQ=255)

                                                                        agTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTAcc                         >  1:4680557/1‑101 (MQ=255)

                                                                         gTTCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCt                        <  1:4679890/101‑1 (MQ=255)

                                                                          ttCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTg                       >  2:4680563/1‑101 (MQ=255)

                                                                          ttCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTg                       >  2:4680566/1‑101 (MQ=255)

                                                                          ttCGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTg                       >  1:4680567/1‑101 (MQ=255)

                                                                            cGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGa                     <  2:4680210/101‑1 (MQ=255)

                                                                            cGCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGa                     >  1:4680571/1‑101 (MQ=255)

                                                                             gCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTAc                          <  1:4680236/95‑1 (MQ=255)

                                                                             gCAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGAt                    <  2:4680046/101‑1 (MQ=255)

                                                                              cAGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATa                   >  2:4680576/1‑101 (MQ=255)

                                                                               aGCCGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATaa                  >  1:4680577/1‑101 (MQ=255)

                                                                                 ccGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATaaaa                >  2:4680582/1‑101 (MQ=255)

                                                                                 ccGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATaaaa                <  1:4680120/101‑1 (MQ=255)

                                                                                 ccGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATaaaa                <  1:4680015/101‑1 (MQ=255)

                                                                                  cGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAAc               <  1:4680034/101‑1 (MQ=255)

                                                                                  cGTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAAc               <  2:4680112/101‑1 (MQ=255)

                                                                                   gTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACg              >  1:4680587/1‑101 (MQ=255)

                                                                                   gTAGGCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACg              >  2:4680586/1‑101 (MQ=255)

                                                                                      ggCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCt           >  1:4680591/1‑101 (MQ=255)

                                                                                       gCTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCtt          <  2:4680083/101‑1 (MQ=255)

                                                                                        cTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTg         <  2:4680016/101‑1 (MQ=255)

                                                                                        cTTTAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTg         <  1:4679992/101‑1 (MQ=255)

                                                                                         tttAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTgc        >  2:4680594/1‑101 (MQ=255)

                                                                                         tttAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTgc        >  2:4680597/1‑101 (MQ=255)

                                                                                          ttAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTgcg       <  2:4680404/101‑1 (MQ=255)

                                                                                          ttAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTgcg       >  1:4680599/1‑101 (MQ=255)

                                                                                           tAGCAGTGCGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTgcgc      >  2:4680603/1‑101 (MQ=255)

GTATGTTTCGCGCTGCCGTGCACGAATGGATTAATCCCCATCGCCACGCCTTCCGGGTTTTTATGCACCAGCAGTTCGCAGCCGTAGGCTTTAGCAGTACGATCGCGGAACTCATACATCTCGCGGAATTTCCAGCCGGTATCGACATGCAGCAACGGGAAAGGCAGCGTACCTGGATAAAACGCCTTGCGCGCCA  >  NC_000913/2875998‑2876193

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