BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,421	2 bp→GC	coding (381‑382/903 nt)	rarD ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,001,421	0	C	G	100.0%	440.1 / NA	137	G128R (GGA→CGA)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (0/0); new base G (63/74); total (63/74)
*	NC_000913	4,001,422	0	G	C	100.0%	437.2 / NA	137	D127E (GAC→GAG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base C (63/74); total (63/74)

Only 100 of 144 total aligned reads displayed.

TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508

||

taataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAAtcgctcgc                                                                                    <  1:1574761/93‑1 (MQ=255)

 aataGGTAGCGAACAAAAAGTCCACAGCTGAACTAACAGGCGACATATCGCCATAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGcc                                                                                   <  2:2401612/93‑1 (MQ=255)

  ataGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCg                                                                                  >  2:3448897/1‑93 (MQ=255)

     ggTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAgg                                                                               <  2:1517197/93‑1 (MQ=255)

      gTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGa                                                                              <  1:3999832/93‑1 (MQ=255)

        aGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaa                                                                            <  1:2368185/93‑1 (MQ=255)

         gCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGaaaa                                                                           <  2:2962037/93‑1 (MQ=255)

          cGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAat                                                                          <  1:1131443/93‑1 (MQ=255)

           gAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatc                                                                         >  1:1253455/1‑93 (MQ=255)

            aaCGAAAAGTCCACAGCGCGACTAACACGCCACATATCGCCAGAATCACAGCCAGCCATTGCATCCGGCTGAATCGCTCGCCGAGGAAAatca                                                                        <  1:1897523/93‑1 (MQ=255)

            aaCCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatca                                                                        <  1:3662401/93‑1 (MQ=255)

             aCCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcat                                                                       >  1:1866523/1‑93 (MQ=255)

             aCCAAAAGGCCACAGCTGGACTAACACGCCACAAATCGCCAGAATCACCGCCAGCCAGTGCATCCGGCGGAACCGCGCGCCGAGGAAAaacac                                                                       >  1:3608452/1‑92 (MQ=255)

              ccAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAatcatc                                                                      <  2:3475897/93‑1 (MQ=255)

               cAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcc                                                                     <  1:17502/93‑1 (MQ=255)

                aaaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATccc                                                                    <  1:3840933/93‑1 (MQ=255)

                 aaaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATcccc                                                                   >  1:3405782/1‑93 (MQ=255)

                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  1:3276801/93‑1 (MQ=255)

                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  <  1:3276648/93‑1 (MQ=255)

                  aaGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCa                                                                  >  1:3975373/1‑93 (MQ=255)

                   aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                                                 <  2:1457517/93‑1 (MQ=255)

                   aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                                                 <  2:772335/93‑1 (MQ=255)

                   aGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAg                                                                 <  1:913/93‑1 (MQ=255)

                    gTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGc                                                                >  2:3114273/1‑93 (MQ=255)

                      ccACAGCTGTACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcac                                                              <  1:1441985/93‑1 (MQ=255)

                      ccACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGcac                                                              >  1:4270487/1‑93 (MQ=255)

                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            <  2:2429813/93‑1 (MQ=255)

                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            >  2:262589/1‑93 (MQ=255)

                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            >  1:202466/1‑93 (MQ=255)

                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            >  1:1607450/1‑93 (MQ=255)

                        acaGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACaa                                                            >  1:112298/1‑93 (MQ=255)

                           gCTGGACTAACACGCCACATATCGCCAGAATCACCGCTAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGt                                                         >  1:1187715/1‑93 (MQ=255)

                             tGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTc                                                       >  1:3657780/1‑93 (MQ=255)

                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                                     >  2:3849723/1‑93 (MQ=255)

                               gACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAc                                                     >  2:3187193/1‑93 (MQ=255)

                                aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                                    >  1:3080783/1‑93 (MQ=255)

                                aCTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCAcc                                                    <  1:386390/93‑1 (MQ=255)

                                 cTAACACGCCACATATCTCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   <  2:1970411/93‑1 (MQ=255)

                                 cTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   >  1:2511539/1‑93 (MQ=255)

                                 cTAACACGCCACATATCGCCAGAATCACCCCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCa                                                   <  1:920867/93‑1 (MQ=255)

                                 cTAACACACCACATATCACCAGAATCACAGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCAACCCCAGCACAATGTTCACCa                                                   <  2:2371457/93‑1 (MQ=255)

                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCGTCCCCAGCACAATGTTCACCAg                                                  >  2:1634735/1‑93 (MQ=255)

                                  tAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAg                                                  >  1:4592486/1‑93 (MQ=255)

                                   aaCACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAgg                                                 >  1:2941899/1‑92 (MQ=255)

                                      acGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                                              >  2:233020/1‑93 (MQ=255)

                                      acGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                                              <  2:2518440/93‑1 (MQ=255)

                                      acGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCggg                                              >  2:3842003/1‑93 (MQ=255)

                                        gCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGtt                                            >  2:4705599/1‑93 (MQ=255)

                                         ccACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTta                                           <  1:1372540/93‑1 (MQ=255)

                                          cacaTATCGACAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaa                                          <  1:377795/93‑1 (MQ=255)

                                           acaTATCGCCAGAATCACCGCCGGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaat                                         >  1:812467/1‑93 (MQ=255)

                                            caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaata                                        <  1:4276280/93‑1 (MQ=255)

                                            caTATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTtaata                                        <  2:3108247/93‑1 (MQ=255)

                                              tatCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATaaa                                      <  1:1698887/93‑1 (MQ=255)

                                               atCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAg                                     >  1:698801/1‑93 (MQ=255)

                                                tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                                    <  1:1021962/93‑1 (MQ=255)

                                                tCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                                    <  2:4563942/93‑1 (MQ=255)

                                                ccaccAGAATCACCGCTAGCCATTGCATCCGGCGGCATCGCCCGCCGAGGAAAATCATCCCCAGCAGAATGTTCACCAGCGGGTTAATAAAGt                                    <  2:4091179/90‑1 (MQ=255)

                                                ccGCCAGAATCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGt                                    <  2:3614724/92‑1 (MQ=255)

                                                 cGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTa                                   <  2:1813083/93‑1 (MQ=255)

                                                    cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAAcc                                >  2:1757142/1‑93 (MQ=255)

                                                    cAGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGGAAca                                >  1:3478645/1‑92 (MQ=255)

                                                     aGAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCa                               <  2:4137989/93‑1 (MQ=255)

                                                      gAATCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCaa                              >  1:3900833/1‑93 (MQ=255)

                                                       aaTCCCCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             <  2:1248887/93‑1 (MQ=255)

                                                       aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             <  2:2953381/93‑1 (MQ=255)

                                                       aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             <  2:2720692/93‑1 (MQ=255)

                                                       aaTCACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAg                             >  2:2839204/1‑93 (MQ=255)

                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCt                          >  2:2797986/1‑93 (MQ=255)

                                                          cACCGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGAAACCAAGGCt                          >  1:4098809/1‑93 (MQ=255)

                                                           cccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTc                         <  2:252584/92‑1 (MQ=255)

                                                            ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATCAAGTAACCAAGGCTcg                        >  2:1645667/1‑92 (MQ=255)

                                                            ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg                        >  2:563726/1‑92 (MQ=255)

                                                            ccGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcg                        >  2:4156106/1‑92 (MQ=255)

                                                             cGCCAGCCATTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgc                       <  2:2628979/93‑3 (MQ=255)

                                                              gccagccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTcgct                      <  1:1062458/93‑4 (MQ=255)

                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   <  2:428401/93‑1 (MQ=255)

                                                                 agccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTcc                   <  2:1330816/93‑1 (MQ=255)

                                                                   ccaTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAg                 <  1:1371727/93‑1 (MQ=255)

                                                                    caTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGc                <  1:357146/93‑1 (MQ=255)

                                                                     aTTGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCa               <  1:3764844/93‑1 (MQ=255)

                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  1:2101339/1‑93 (MQ=255)

                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  1:3262185/1‑93 (MQ=255)

                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  1:3381054/1‑93 (MQ=255)

                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  2:3235281/1‑93 (MQ=255)

                                                                      ttGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              >  2:1466374/1‑93 (MQ=255)

                                                                      atGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCat              <  2:3632647/92‑1 (MQ=255)

                                                                       tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata             <  1:3095957/93‑1 (MQ=255)

                                                                       tGCATCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCata             <  1:3145117/93‑1 (MQ=255)

                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          >  1:4622384/1‑93 (MQ=255)

                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          <  1:1081199/93‑1 (MQ=255)

                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          <  2:3390975/93‑1 (MQ=255)

                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          >  2:2895138/1‑93 (MQ=255)

                                                                          aTCCGGCGGAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATgg          <  1:4087845/93‑1 (MQ=255)

                                                                             cggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGa       <  1:32496/93‑1 (MQ=255)

                                                                              ggcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAt      >  2:2864321/1‑93 (MQ=255)

                                                                               gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt     <  2:857886/93‑1 (MQ=255)

                                                                               gcggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGAtt     <  1:2980114/93‑1 (MQ=255)

                                                                                cggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTg    <  1:1838256/93‑1 (MQ=255)

                                                                                 ggAATCGCTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCGCTTCCAGCATATGGTGATTGt   <  1:530839/93‑1 (MQ=255)

||

TAATAGGTAGCGAACCAAAAGTCCACAGCTGGACTAACACGCCACATATCGCCAGAATCACCGCCAGCCATTGCATCCGGCGGAATCCGTCGCCGAGGAAAATCATCCCCAGCACAATGTTCACCAGCGGGTTAATAAAGTAACCAAGGCTCCGTTCCAGCATATGGTGATTGTT  >  NC_000913/4001334‑4001508

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