BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,695	T→C	K581R (AAG→AGG)	frvR ←	putative frv operon regulatory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,085,695	0	T	C	100.0%	906.7 / NA	277	K581R (AAG→AGG)	frvR	putative frv operon regulatory protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (147/130); total (147/130)

Only 100 of 285 total aligned reads displayed.

GAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTA  >  NC_000913/4085605‑4085784

gAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCtgt                                                                                        <  2:1846558/93‑1 (MQ=255)

gAAATCGTTGATGACTATAATAGCCTGCGTTACACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCtgt                                                                                        <  2:2670077/93‑1 (MQ=255)

  aaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGa                                                                                      <  1:2772000/93‑1 (MQ=255)

   aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGaa                                                                                     >  2:1574290/1‑93 (MQ=255)

   aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGaa                                                                                     <  1:4749744/93‑1 (MQ=255)

   aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGaa                                                                                     <  1:3038175/93‑1 (MQ=255)

   aTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGaa                                                                                     <  1:2850135/93‑1 (MQ=255)

     ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGt                                                                                   <  1:962590/93‑1 (MQ=255)

         gatgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACt                                                                               >  1:3759361/1‑93 (MQ=255)

          atgaTTATAATAGCCTGTTTAAAGTGTTTTTGATAATAAATGCAGTGTCGGAGGCGGCGC‑ACCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                                              <  2:4288579/93‑1 (MQ=255)

          atgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTc                                                                              <  2:814677/93‑1 (MQ=255)

           tgaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTcc                                                                             >  1:3637532/1‑93 (MQ=255)

            gaCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCa                                                                            >  2:4867976/1‑93 (MQ=255)

             attataATAGCCTGCTTTAACCGTTTTTGATAAAAAATGCAGTGTCGGATGCGGCGC‑ACCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                                           <  2:2214/91‑1 (MQ=255)

             aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                                           <  1:2397729/93‑1 (MQ=255)

              cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                                          >  1:4042530/1‑93 (MQ=255)

               tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                                         >  2:2062896/1‑93 (MQ=255)

               tataATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                                         >  2:1372792/1‑93 (MQ=255)

                ataataGCGTGCGTTAAACTTTTTGGATGATAAATGCAGTGTCGGAGGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                                        <  2:3345016/93‑1 (MQ=255)

                 taataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                                       >  1:4775905/1‑93 (MQ=255)

                  aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATg                                                                      >  2:5220116/1‑93 (MQ=255)

                    taGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGGAACTCCATAAATGcc                                                                    >  2:369171/1‑93 (MQ=255)

                     aGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCt                                                                   <  2:1666545/93‑1 (MQ=255)

                       cgtGCGTCAACCTTTTTTGATGATAAGTGCTGGGTGGGATGGGGCAC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                                                 <  1:1089285/91‑1 (MQ=255)

                       ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                                                 <  1:4283589/93‑1 (MQ=255)

                           cGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                                             <  1:2659113/93‑1 (MQ=255)

                             ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                           <  1:187117/93‑1 (MQ=255)

                             ttAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                                           >  1:1843006/1‑93 (MQ=255)

                              tAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                                          >  1:1641972/1‑93 (MQ=255)

                               aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                                         >  2:1387785/1‑93 (MQ=255)

                               aaaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                                         <  2:1983849/93‑1 (MQ=255)

                               aaaCGTTTTTGATGATAAATGCAGCGTCGGATGCGGCGCG‑TCCGACACAGTTCACCCCCTGTGAAGCAACTCCATAAATACCTCATATCCTGt                                                         >  2:3010936/1‑93 (MQ=255)

                                aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                                        >  1:5282092/1‑93 (MQ=255)

                                aaCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                                        >  1:4787607/1‑93 (MQ=255)

                                   gTTTTTGATGATTAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGCAGTAACTCCATAAATGCCTCATATCCTGTTAAc                                                     >  1:3408924/1‑93 (MQ=255)

                                    tttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                                                    >  2:6175217/1‑93 (MQ=255)

                                     ttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                                                   <  1:5558129/93‑1 (MQ=255)

                                       ttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                                 >  1:1598059/1‑93 (MQ=255)

                                       ttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCgg                                                 <  2:5799549/93‑1 (MQ=255)

                                         tatgaaAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCCCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                               <  2:2301346/92‑1 (MQ=255)

                                         gatgatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCa                                               >  2:3462158/1‑93 (MQ=255)

                                           tgatAAATGCAGTTTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                             <  2:5611253/93‑1 (MQ=255)

                                           tgatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                                             <  2:4557492/93‑1 (MQ=255)

                                            gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                                            >  2:450894/1‑93 (MQ=255)

                                            gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                                            <  2:497301/93‑1 (MQ=255)

                                            gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                                            <  2:3610379/93‑1 (MQ=255)

                                            gatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                                            >  1:2092178/1‑93 (MQ=255)

                                            gatAAATGCAGTGTCGGATGCGACGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaata                                            <  1:5078287/93‑1 (MQ=255)

                                              tAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                                          >  1:5049899/1‑93 (MQ=255)

                                               aaaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtc                                         >  1:3583364/1‑93 (MQ=255)

                                                aaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCCGTGAAGTAACTCCATAAATGCCTCATATCCTGTTCACCCGGCCATAAtct                                        >  2:3254971/1‑93 (MQ=255)

                                                 aTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                                       >  2:2278844/1‑93 (MQ=255)

                                                  tGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                      >  2:4234804/1‑93 (MQ=255)

                                                  tGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                                      <  2:3867960/93‑1 (MQ=255)

                                                    cAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc                                    <  2:4699669/93‑1 (MQ=255)

                                                    cAGTGCGGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgc                                    <  1:568304/93‑1 (MQ=255)

                                                     aGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc                                   <  2:5277289/93‑1 (MQ=255)

                                                      gtgtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgccg                                  >  1:2229845/1‑93 (MQ=255)

                                                        gtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa                                <  2:5932409/93‑1 (MQ=255)

                                                         tCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                               >  1:6244605/1‑93 (MQ=255)

                                                         tCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                               >  2:4775188/1‑93 (MQ=255)

                                                         tCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                               >  1:5950511/1‑93 (MQ=255)

                                                         tCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGAt                               <  1:3249069/93‑1 (MQ=255)

                                                          cGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATg                              <  1:5148091/93‑1 (MQ=255)

                                                           ggATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc                             >  2:2360508/1‑93 (MQ=255)

                                                             aTGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTg                           >  1:3324494/1‑93 (MQ=255)

                                                              tGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa                          <  1:189273/93‑1 (MQ=255)

                                                                cggcgC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGaca                        <  1:952608/93‑1 (MQ=255)

                                                                   gcgcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca                      <  2:3752626/89‑1 (MQ=255)

                                                                   gcgcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGAcaca                      >  1:1061407/5‑93 (MQ=255)

                                                                     gcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                    >  1:3071512/3‑93 (MQ=255)

                                                                     gcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                    >  2:6215893/3‑93 (MQ=255)

                                                                     gcaTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAat                    >  1:1736564/3‑93 (MQ=255)

                                                                      caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  1:4543353/92‑1 (MQ=255)

                                                                      caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   >  2:3582872/2‑93 (MQ=255)

                                                                      caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   >  2:3046922/2‑93 (MQ=255)

                                                                      caTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAata                   <  1:2949501/92‑1 (MQ=255)

                                                                       aTCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAatat                  <  1:1240505/93‑1 (MQ=255)

                                                                        tCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                 >  2:4065894/1‑93 (MQ=255)

                                                                        tCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATg                 >  1:2289098/1‑93 (MQ=255)

                                                                         ccGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc                <  2:6133606/93‑1 (MQ=255)

                                                                         ccGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGc                <  2:3956272/93‑1 (MQ=255)

                                                                          cGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt               >  1:5736159/1‑93 (MQ=255)

                                                                          cGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCt               <  2:6291850/93‑1 (MQ=255)

                                                                           gACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt              <  1:2635908/93‑1 (MQ=255)

                                                                           gACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCtt              >  2:1950021/1‑93 (MQ=255)

                                                                            acacAGTTCACCCCCTGTGAAGTAACTCGATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTg             <  1:646832/93‑1 (MQ=255)

                                                                             cccaGTGCACCCCCTGTGAAGTAACTCCATAAAGGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGc            <  1:4689850/91‑1 (MQ=255)

                                                                             cacaGTTCACCCCCTGTGAAGTAACTCCATAAATTCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGc            >  1:5446561/1‑93 (MQ=255)

                                                                              acaGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           >  1:4493519/1‑93 (MQ=255)

                                                                              acaGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGcc           <  1:824589/93‑1 (MQ=255)

                                                                                aGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAg         <  1:3866352/93‑1 (MQ=255)

                                                                                  ttCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt       >  2:4340703/1‑93 (MQ=255)

                                                                                  ttCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt       <  2:2904437/93‑1 (MQ=255)

                                                                                  ttCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTATTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGAt       <  1:6044779/93‑1 (MQ=255)

                                                                                     aCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc    <  2:733483/93‑1 (MQ=255)

                                                                                     aCCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGc    >  2:775406/1‑93 (MQ=255)

                                                                                      ccccCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt   <  1:4084259/93‑1 (MQ=255)

                                                                                      ccccCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCt   >  2:5048933/1‑93 (MQ=255)

                                                                                       ccccTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTa  >  2:1259135/1‑93 (MQ=255)

GAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACACAATATGCTTGCCAGATAGCTA  >  NC_000913/4085605‑4085784

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