BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,001,717	+C	coding (86/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,717	1	.	C	100.0%	993.3 / NA	319	T29R (ACG→AGG)	rarD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (145/174); total (145/174)

Only 100 of 358 total aligned reads displayed.

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACG‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                       <  1:1019880/93‑3 (MQ=255)

  cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgct                                                                                       <  2:4574289/93‑3 (MQ=255)

   tcatcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACgctg                                                                                      >  1:2337846/1‑90 (MQ=255)

      tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                                   >  1:449898/1‑93 (MQ=255)

      tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                                   >  1:1530112/1‑93 (MQ=255)

      tcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                                   >  1:2928623/1‑93 (MQ=255)

       caGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAATTAATCAACTTGTAGTACGCTGgcgc                                                                                  >  2:4974041/1‑93 (MQ=255)

        aGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCt                                                                                 <  2:5620820/93‑1 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCta                                                                                <  1:2735975/93‑1 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCta                                                                                <  2:7566549/93‑1 (MQ=255)

         gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCta                                                                                <  1:5601178/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                              <  1:1242471/93‑1 (MQ=255)

           accaccATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCtata                                                                              <  1:6883393/93‑1 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCGCGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                             <  2:3106471/93‑1 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGGGCTATAc                                                                             >  2:6709757/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                             >  2:7811738/1‑93 (MQ=255)

            ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAc                                                                             <  1:5938816/93‑1 (MQ=255)

             caccaTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATAcc                                                                            <  1:5546778/93‑1 (MQ=255)

                 aTAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGGAGTAAATCAACTTTAAGTACGCTGGCGCTATACCCCaa                                                                        >  1:5002972/1‑93 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGGCCGTAGTAAATAAACTTTAAGTAAGCTGGCGCTATACCCCaaa                                                                       >  1:1755246/1‑93 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaa                                                                       >  2:1919281/1‑93 (MQ=255)

                  tAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGCAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaa                                                                       >  2:4813049/1‑93 (MQ=255)

                    aaGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaata                                                                     >  2:5084193/1‑93 (MQ=255)

                     aGAAAAACGACCAGAACACGCGATGCCGCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataa                                                                    <  1:2507138/93‑1 (MQ=255)

                      gAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaa                                                                   >  2:2469053/1‑93 (MQ=255)

                        aaaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaaa                                                                 >  2:2559283/1‑92 (MQ=255)

                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                <  1:6095738/93‑1 (MQ=255)

                         aaaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCaaataaaata                                                                <  1:2575969/93‑1 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAaataaaataa                                                               <  2:6426281/93‑1 (MQ=255)

                          aaCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAaataaaataa                                                               <  2:4061496/93‑1 (MQ=255)

                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                                              <  2:7924726/93‑1 (MQ=255)

                           aCGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAg                                                              >  1:2028617/1‑93 (MQ=255)

                            cGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGc                                                             <  2:2690381/93‑1 (MQ=255)

                             gACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAACg                                                            >  1:6728611/1‑93 (MQ=255)

                              atcAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                           <  2:6342437/91‑1 (MQ=255)

                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCgg                                                           <  2:2847443/93‑1 (MQ=255)

                              aCCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTATACCCCAAATAAAATAAGcag                                                           >  1:3436621/1‑91 (MQ=255)

                               ccAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTGGAGGACGCTGGCGCTATACCCCAAATAAAATAAACGGc                                                          >  2:4447889/1‑93 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGcca                                                        >  1:5223653/1‑91 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        >  1:3002676/1‑93 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        >  2:3244969/1‑93 (MQ=255)

                                 aGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCaa                                                        >  1:3473366/1‑93 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGGACGCTGGGGCTATACCCCAAATAAAATAAGCGGGAag                                                       >  1:1971481/1‑93 (MQ=255)

                                  gATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGGAAATCAACTTGAAGTAAGCTGGAGCTATACCCCAAATAAAATAAGCGGCAag                                                       >  1:3647049/1‑93 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  2:6045005/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  2:3606634/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      <  1:4854740/93‑1 (MQ=255)

                                   aTCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAaga                                                      >  1:6243416/1‑93 (MQ=255)

                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                     <  2:7409147/93‑1 (MQ=255)

                                    tCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAagag                                                     >  2:7657490/1‑93 (MQ=255)

                                        gcgATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAg                                                 <  1:6307092/93‑1 (MQ=255)

                                           aTGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTaa                                              >  1:6341825/1‑93 (MQ=255)

                                            tGCGTCAGGATTTGATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             >  1:7158855/1‑93 (MQ=255)

                                            tGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGGACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAAt                                             >  2:617396/1‑93 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAATTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:7849262/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:5692905/93‑1 (MQ=255)

                                             gCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  2:1212317/93‑1 (MQ=255)

                                             gCGGCAGAATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATa                                            <  1:7651348/93‑1 (MQ=255)

                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATCCCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                           <  1:1310575/93‑1 (MQ=255)

                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                           <  2:7263641/93‑1 (MQ=255)

                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                           <  1:4058565/93‑1 (MQ=255)

                                              cGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATAc                                           <  1:2181835/93‑1 (MQ=255)

                                               gTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACg                                          <  1:5290982/93‑1 (MQ=255)

                                                tCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCCAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                         >  2:4730920/1‑93 (MQ=255)

                                                  aGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGccc                                       <  1:4657604/93‑1 (MQ=255)

                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCCAATAAAATAAGCGGCAAGGGCGAGTAATACGCCCt                                      >  2:1725179/1‑93 (MQ=255)

                                                   ggATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                      >  1:3204460/1‑93 (MQ=255)

                                                       ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                  <  2:3401883/93‑1 (MQ=255)

                                                       ttCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcg                                  <  2:3801904/93‑1 (MQ=255)

                                                         cATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCcgcg                                <  2:2102724/93‑1 (MQ=255)

                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                               <  2:2960775/93‑1 (MQ=255)

                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                               <  2:6936807/93‑1 (MQ=255)

                                                          aTCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgt                               >  2:5306508/1‑93 (MQ=255)

                                                           tCGGCGTGTACGTAGTAAATCAACGTGAAGTACGCTGGCGCTATACCCCAAAAAAAATAAGCGGCAAGAGCGAGTAATACGTCCTGCCGCgtt                              <  2:3995965/93‑1 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              <  2:3721587/93‑1 (MQ=255)

                                                           tCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                              >  1:899991/1‑93 (MQ=255)

                                                            cggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                             <  1:2454679/93‑1 (MQ=255)

                                                            cggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                             <  2:6789853/93‑1 (MQ=255)

                                                             ggcggGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                            >  1:3763655/1‑93 (MQ=255)

                                                                 ggCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGtttt                        <  1:3871462/93‑1 (MQ=255)

                                                                   cACGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc                      <  2:2363513/93‑1 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:1568756/1‑93 (MQ=255)

                                                                    aCGTAGTAAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGca                     >  2:5759112/1‑93 (MQ=255)

                                                                      gtagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                   >  1:4183147/1‑93 (MQ=255)

                                                                       tagtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcc                  >  2:6549716/1‑93 (MQ=255)

                                                                        agtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                 >  1:4912678/1‑93 (MQ=255)

                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatccat                >  1:7495436/1‑93 (MQ=255)

                                                                         gtaAATCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGGAATACGCCCTGCCGCGTGTGTGGTGcaaccac                >  1:7002612/1‑92 (MQ=255)

                                                                           aaaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg              >  1:1203325/1‑93 (MQ=255)

                                                                            aaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             >  1:7632789/1‑93 (MQ=255)

                                                                            aaTCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCgg             <  1:6498090/93‑1 (MQ=255)

                                                                              tCAACTTGAAGTACGCTGGCGCTATACCCCAAATAAAAAAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGCTTGTTTTGCATCCATCGGGa           >  2:5434346/1‑93 (MQ=255)

                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATTAGCGGCAAGAGCGAGGAAAACGCCCCGCTGCGTTTGTTTTGCAGCCTTCGGGTAt         >  2:7041654/1‑93 (MQ=255)

                                                                                aaCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAAt         <  1:1643251/93‑1 (MQ=255)

                                                                                 aCTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa        >  1:7078743/1‑93 (MQ=255)

                                                                                  cTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       >  2:3470725/1‑93 (MQ=255)

                                                                                  cTTGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATAc       <  1:1369330/93‑1 (MQ=255)

                                                                                   ttGAAGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt      >  1:3286567/1‑93 (MQ=255)

                                                                                       aGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt  >  1:774283/1‑93 (MQ=255)

                                                                                       aGTACGCTGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCAtt  >  1:414597/1‑93 (MQ=255)

TGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACG‑TGGCGCTATACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATT  >  NC_000913/4001625‑4001803

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