BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,085,675	Δ1 bp	intergenic (+38/+13)	yiiG → / ← frvR	hypothetical protein/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,675	0	G	.	100.0%	868.0 / NA	191	intergenic (+38/+13)	yiiG/frvR	hypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (123/68); total (123/68)

Only 100 of 215 total aligned reads displayed.

GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762

gCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCg                                                                                    >  1:732291/1‑93 (MQ=255)

 cATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGa                                                                                   >  2:924812/1‑93 (MQ=255)

   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGaca                                                                                 >  2:1412624/1‑93 (MQ=255)

   ttACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGaca                                                                                 >  1:4111954/1‑93 (MQ=255)

    tACGCAAATAAAACAAAATTGTTAATGATTATAATAGCCTGCGTTAATTT‑TTTTTAATAAAAAAGGCAGTGTCGGAGGCGGCGC‑ACCCGacac                                                                                <  2:856498/93‑1 (MQ=255)

     aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGAcaca                                                                               >  2:1276049/1‑93 (MQ=255)

     aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGAcaca                                                                               <  1:597748/93‑1 (MQ=255)

       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGt                                                                             <  2:2621732/93‑1 (MQ=255)

       gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGt                                                                             >  1:293632/1‑93 (MQ=255)

       gcgAGTACAACGAAAGGGTTGATGATTATAATAGCCTGCGTTAAACT‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGt                                                                             <  2:3211385/93‑1 (MQ=255)

         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                                           <  2:2093182/93‑1 (MQ=255)

         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                                           >  2:1830743/1‑93 (MQ=255)

         gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                                           >  1:4226566/1‑93 (MQ=255)

          aggaCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGCGTCGGATGCGGCG‑TCTCCGACACAGTTCa                                                                          >  2:4655234/4‑93 (MQ=255)

          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTcc                                                                          >  2:3130569/1‑92 (MQ=255)

          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCa                                                                          >  1:2146803/1‑93 (MQ=255)

          aGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGAAGTGTCGGATGCGGCGC‑ATCCGACACAGTTcc                                                                          >  1:3433200/1‑92 (MQ=255)

           gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAc                                                                         >  1:519490/1‑93 (MQ=255)

           gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAc                                                                         >  1:3163600/1‑93 (MQ=255)

            tACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAcc                                                                        <  1:221857/93‑1 (MQ=255)

             aCAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccc                                                                       >  2:431568/1‑93 (MQ=255)

              cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAcccc                                                                      >  2:1909710/1‑93 (MQ=255)

               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccc                                                                     >  1:3672460/1‑92 (MQ=255)

               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccc                                                                     >  1:2614074/1‑92 (MQ=255)

               aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccc                                                                     >  2:4406428/1‑92 (MQ=255)

                aCGAAATGTTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ACCCGACACAGTTCAccccct                                                                    <  2:2720158/93‑3 (MQ=255)

                aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccct                                                                    >  2:2425285/1‑91 (MQ=255)

                aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccct                                                                    >  1:2112083/1‑91 (MQ=255)

                  gAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCtgt                                                                  >  2:4413312/1‑93 (MQ=255)

                    aaTGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGa                                                                >  2:2497191/1‑93 (MQ=255)

                     aTGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGaa                                                               <  1:405776/93‑1 (MQ=255)

                     aTGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGaa                                                               <  1:3814494/93‑1 (MQ=255)

                      tGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAg                                                              <  1:3694819/93‑1 (MQ=255)

                       ggTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGt                                                             >  1:453854/1‑93 (MQ=255)

                          tgatgaCTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAAc                                                          >  1:4636185/1‑93 (MQ=255)

                          tgatgaCTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAAc                                                          >  1:1341988/1‑93 (MQ=255)

                             tgaCTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTcc                                                       >  2:1505633/1‑93 (MQ=255)

                              gaCTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTCACTCCa                                                      >  1:4005056/1‑93 (MQ=255)

                               aCTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                     <  1:1927065/93‑1 (MQ=255)

                                cTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    <  1:3786470/93‑1 (MQ=255)

                                cTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    >  2:4474936/1‑93 (MQ=255)

                                 tataATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaa                                                   >  1:2553543/1‑93 (MQ=255)

                                  aaaatgGCCTGC‑TTAAACGTTTTTTGATGAAAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  <  2:1720452/91‑1 (MQ=255)

                                   taataGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 <  2:2670433/93‑1 (MQ=255)

                                   taataGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 <  2:3620739/93‑1 (MQ=255)

                                   taataGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAAt                                                 >  1:2011068/1‑93 (MQ=255)

                                        gCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                            >  1:3903152/1‑93 (MQ=255)

                                        gCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTc                                            <  1:1192263/93‑1 (MQ=255)

                                         ccTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                           <  1:3471869/93‑1 (MQ=255)

                                           tGCGTTAACCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:4568581/1‑93 (MQ=255)

                                           tGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGCAGTAACTCCATAAATGCCTCata                                         >  2:3537526/1‑93 (MQ=255)

                                           tGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:3432876/1‑93 (MQ=255)

                                           tGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  1:2035415/1‑93 (MQ=255)

                                             cGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       <  1:3087489/93‑1 (MQ=255)

                                             cGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATc                                       <  2:3867406/93‑1 (MQ=255)

                                              gTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      >  1:3056393/1‑93 (MQ=255)

                                              gTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGAGGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATcc                                      <  1:1861167/93‑1 (MQ=255)

                                               ttAAACT‑TTTTGGAGAATAAAGGCGGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     <  2:2564357/93‑1 (MQ=255)

                                               ttAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     >  2:2229190/1‑93 (MQ=255)

                                               ttAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCt                                     <  1:2258375/93‑1 (MQ=255)

                                                 aaaCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGt                                   <  1:1777950/93‑1 (MQ=255)

                                                  aaCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  >  1:1999376/1‑93 (MQ=255)

                                                  aaCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  >  1:4413657/1‑93 (MQ=255)

                                                  aaCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  >  2:3339382/1‑93 (MQ=255)

                                                  aaCG‑TTTTTGATAATAAATGCAGTGTCGGAGGCGGCGC‑ACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGtt                                  <  2:811844/93‑1 (MQ=255)

                                                   aCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  1:3060180/1‑93 (MQ=255)

                                                   aCG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  1:3630753/1‑93 (MQ=255)

                                                    cG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                >  1:1231329/1‑93 (MQ=255)

                                                    cG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                >  2:3156936/1‑93 (MQ=255)

                                                     g‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTTACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               >  1:3253268/1‑93 (MQ=255)

                                                       tttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAcc                              <  1:382575/93‑1 (MQ=255)

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

                                                        ttttGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAccc                             >  2:3659381/1‑93 (MQ=255)

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

                                                             atgatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTACATAAATGCCTCATATCCTGTTAACCCGGCaa                        >  1:2342117/1‑93 (MQ=255)

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

                                                              tgatAAATGCAGTGTCGGATGCGGC‑CCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaat                       <  1:1625416/93‑1 (MQ=255)

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

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

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

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

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

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

                                                                atAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  1:312463/1‑93 (MQ=255)

                                                                 tAAATGCAGTGTCGGATGCGGCGCG‑TCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    <  2:4156785/93‑1 (MQ=255)

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

                                                                 tAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataat                    <  1:4150128/93‑1 (MQ=255)

                                                                   aaTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtct                  <  1:1030548/93‑1 (MQ=255)

                                                                    aTGCAGTGTCGGAGGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                 <  1:990520/93‑1 (MQ=255)

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

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

                                                                     tGCAGTGTCGGAGGCGCCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTcc                <  2:3941235/93‑1 (MQ=255)

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

                                                                              ggATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  2:3764903/93‑1 (MQ=255)

                                                                              ggATGCGGCGC‑ACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  1:2432943/93‑1 (MQ=255)

                                                                              ggAGGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGc       <  1:3184596/93‑1 (MQ=255)

                                                                               gATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      >  1:2238874/1‑93 (MQ=255)

                                                                               gATGCGGCGC‑ATCCGACACAGTTCACACCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGCTGcg      >  1:3416690/1‑92 (MQ=255)

                                                                               gATGCGGCGC‑ATCCGAAACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      <  2:2394829/93‑1 (MQ=255)

                                                                                 tGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGa    >  2:3969005/1‑93 (MQ=255)

GCATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACG‑TTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085587‑4085762

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