BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,885,561	(C)_5→4	intergenic (‑319/+40)	ygcB ← / ← cysH	hypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,885,557	0	C	.	100.0%	834.3 / NA	180	intergenic (‑315/+44)	ygcB/cysH	hypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand): ref base C (0/0); new base . (103/77); total (103/77)

Only 100 of 204 total aligned reads displayed.

GCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAC  >  NC_000913/2885466‑2885644

gCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcc                                                                                        >  1:4428240/1‑93 (MQ=255)

gCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcc                                                                                        >  2:207480/1‑93 (MQ=255)

   tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                     >  2:5285526/1‑92 (MQ=255)

   tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                     >  2:5191312/1‑92 (MQ=255)

   tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                     >  2:1582730/1‑92 (MQ=255)

   tCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGcccca                                                                                     >  2:4214685/1‑92 (MQ=255)

    cAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccat                                                                                    >  1:1640773/1‑91 (MQ=255)

     aaGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatc                                                                                   >  2:600447/1‑90 (MQ=255)

     aaGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatc                                                                                   <  2:4452412/93‑4 (MQ=255)

      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                  <  2:2282954/93‑5 (MQ=255)

      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                  >  2:5903662/1‑89 (MQ=255)

      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                  >  1:5469558/1‑89 (MQ=255)

      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                  <  2:1097414/93‑5 (MQ=255)

      aGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGccccatca                                                                                  >  2:1519306/1‑89 (MQ=255)

       gCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑TCCcatcat                                                                                >  2:5389391/1‑93 (MQ=255)

       gCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCcatcat                                                                                <  1:4611315/93‑1 (MQ=255)

        cACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCAtt                                                                               >  2:4676375/1‑93 (MQ=255)

          cccGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTc                                                                             >  2:5079246/1‑93 (MQ=255)

             gACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGa                                                                          >  1:1097957/1‑93 (MQ=255)

              aCGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGAc                                                                         <  2:4579584/93‑1 (MQ=255)

                  cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACagag                                                                     >  1:4201142/1‑93 (MQ=255)

                  cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACagag                                                                     <  1:3642927/93‑1 (MQ=255)

                  cATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACagag                                                                     <  1:2683636/93‑1 (MQ=255)

                  cATTCTATACCTGATAATTCTTCCCGCGGCGTGCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACagag                                                                     <  1:5112257/93‑1 (MQ=255)

                    ttCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGc                                                                   >  1:2226814/1‑93 (MQ=255)

                      cTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGt                                                                 >  2:2978882/1‑93 (MQ=255)

                       tataCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGtt                                                                <  2:199732/93‑1 (MQ=255)

                          aCCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTaa                                                             >  2:1666665/1‑93 (MQ=255)

                          aCCTGATAATTCTGCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTaa                                                             <  2:2580209/93‑1 (MQ=255)

                             tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAAttt                                                          >  1:6669674/1‑93 (MQ=255)

                             tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAAttt                                                          >  1:4116532/1‑93 (MQ=255)

                             tGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAAttt                                                          >  1:4157466/1‑93 (MQ=255)

                              gATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTg                                                         >  1:2161405/1‑93 (MQ=255)

                               aTAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGt                                                        >  2:6033190/1‑93 (MQ=255)

                                tAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTc                                                       <  1:306444/93‑1 (MQ=255)

                                 aaTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                      <  2:3084464/93‑1 (MQ=255)

                                 aaTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTcc                                                      <  2:6309790/93‑1 (MQ=255)

                                  aTTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCg                                                     <  2:3370894/93‑1 (MQ=255)

                                    tcttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGc                                                   >  1:5797154/1‑93 (MQ=255)

                                     cttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCa                                                  >  1:4768967/1‑93 (MQ=255)

                                      ttcCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCaa                                                 >  1:4088907/1‑93 (MQ=255)

                                        cccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAata                                               >  2:1351040/1‑93 (MQ=255)

                                         ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                              >  1:2524010/1‑93 (MQ=255)

                                         ccGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAatat                                              >  2:1940705/1‑93 (MQ=255)

                                          cgcgGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAtt                                             >  2:3746013/1‑93 (MQ=255)

                                           gcggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAttt                                            <  2:5854869/93‑1 (MQ=255)

                                           gcggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATAttt                                            <  2:4124420/93‑1 (MQ=255)

                                            cggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTa                                           <  1:1948345/93‑1 (MQ=255)

                                            cggcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTa                                           <  1:3231346/93‑1 (MQ=255)

                                              gcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                         >  2:2518956/1‑93 (MQ=255)

                                              gcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                         <  1:1892699/93‑1 (MQ=255)

                                              gcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                         >  1:949013/1‑93 (MQ=255)

                                              gcgTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcc                                         >  1:483350/1‑93 (MQ=255)

                                                gTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCt                                       <  1:2033232/93‑1 (MQ=255)

                                                 ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTAcccat                                      <  1:6411780/93‑3 (MQ=255)

                                                 ttCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCtt                                      >  1:4295757/1‑93 (MQ=255)

                                                  tCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                     >  2:5970849/1‑93 (MQ=255)

                                                  tCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                     >  1:1516208/1‑93 (MQ=255)

                                                  tCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTc                                     <  1:2208452/93‑1 (MQ=255)

                                                      gTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCgtgt                                 >  1:625444/1‑93 (MQ=255)

                                                       tttAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTa                                >  1:6688568/1‑93 (MQ=255)

                                                        ttAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTaa                               <  1:5341758/93‑1 (MQ=255)

                                                        ttAAAATCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTaa                               <  2:1811029/93‑1 (MQ=255)

                                                           aaaTCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAAccc                            >  2:2609623/1‑93 (MQ=255)

                                                            aaTCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                           <  1:6643678/93‑1 (MQ=255)

                                                            aaTCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                           >  2:3724353/1‑93 (MQ=255)

                                                            aaTCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                           >  2:611051/1‑93 (MQ=255)

                                                            aaTCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCca                           <  2:2844647/93‑1 (MQ=255)

                                                             aTCAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcac                          >  1:3928062/1‑93 (MQ=255)

                                                              tCAGGAATGCTGCTCGCCAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGCTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcaca                         >  1:4312186/1‑93 (MQ=255)

                                                              tCAGGAATGCTGCTCGCAAAATAAAGACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcaca                         >  1:1373357/1‑93 (MQ=255)

                                                              acAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCcaca                         <  1:814538/92‑1 (MQ=255)

                                                               cAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAt                        >  1:522124/1‑93 (MQ=255)

                                                               cAGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAt                        <  1:6039903/93‑1 (MQ=255)

                                                                aGGAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACAtt                       >  1:1840354/1‑93 (MQ=255)

                                                                 ggAATGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTTACCCACATTc                      >  1:275323/1‑93 (MQ=255)

                                                                   aaTGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTccc                    <  2:6267222/93‑1 (MQ=255)

                                                                    aTGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCt                   <  1:781910/93‑1 (MQ=255)

                                                                    aTGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCt                   >  1:3494864/1‑93 (MQ=255)

                                                                    aTGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCt                   >  1:3782348/1‑93 (MQ=255)

                                                                    aTGCTGCTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCt                   >  1:611072/1‑93 (MQ=255)

                                                                     tgctgcTCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCtt                  <  1:1281493/93‑1 (MQ=255)

                                                                      gctgctCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCttt                 <  1:3997429/93‑1 (MQ=255)

                                                                        tgctCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACCTTCACttttt               >  2:5103624/1‑93 (MQ=255)

                                                                        tgctCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCttttt               >  1:3591564/1‑93 (MQ=255)

                                                                           tCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg            >  1:5286563/1‑93 (MQ=255)

                                                                           tCGCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAg            >  2:4880568/1‑93 (MQ=255)

                                                                             gCAAAATAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGcc          >  1:2418884/1‑93 (MQ=255)

                                                                               aaaaTAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaa        >  2:5633025/1‑93 (MQ=255)

                                                                                aaaTAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCaaa       >  1:4674918/1‑93 (MQ=255)

                                                                                  aTAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGa     >  1:1051936/1‑93 (MQ=255)

                                                                                   tAAAAACG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaa    <  1:405354/93‑1 (MQ=255)

                                                                                    aaaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa   <  1:3168149/93‑1 (MQ=255)

                                                                                    aaaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGaaa   <  1:6109016/93‑1 (MQ=255)

                                                                                     aaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  >  2:5721541/1‑93 (MQ=255)

                                                                                     aaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  <  2:1300699/93‑1 (MQ=255)

                                                                                     aaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  <  1:1293447/93‑1 (MQ=255)

                                                                                     aaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  >  2:3932981/1‑93 (MQ=255)

                                                                                     aaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  <  2:609051/93‑1 (MQ=255)

                                                                                     aaaaCG‑CCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAc  >  1:3899954/1‑93 (MQ=255)

GCCTCAAGCACCCGACGCCATTCTATACCTGATAATTCTTCCCGCGGCGTTCACGTTTAAAATCAGGAATGCTGCTCGCAAAATAAAAACGCCCCCATCATTTCTGACAGAGGCGTTTAATTTGTCCGGCAATATTTACCCTTCGTGTAACCCACATTCCCTTTTTAAGCCAAAGAAAC  >  NC_000913/2885466‑2885644

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