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%	1070.3 / NA	233	intergenic (+38/+13)	yiiG/frvR	hypothetical protein/putative frv operon regulatory protein
Reads supporting (aligned to +/- strand): ref base G (0/0); new base . (143/90); total (143/90)

Only 100 of 247 total aligned reads displayed.

CATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085588‑4085762

   tACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAGTGCAGTGTCGGATGCGGCGC‑ATCCGacac                                                                                <  1:2018759/93‑1 (MQ=255)

    aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGAcaca                                                                               <  1:669172/93‑1 (MQ=255)

    aCGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGAcaca                                                                               <  2:3915206/93‑1 (MQ=255)

      gcgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGt                                                                             >  2:5553339/1‑93 (MQ=255)

       cgAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGtt                                                                            >  2:2783465/1‑93 (MQ=255)

        gAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTc                                                                           >  2:4886435/1‑93 (MQ=255)

          gTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAc                                                                         >  2:2957061/1‑93 (MQ=255)

             cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAcccc                                                                      >  1:6182619/1‑93 (MQ=255)

             cAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCCGCGC‑ATCCGACACAGTAGAcccc                                                                      >  2:4839604/1‑93 (MQ=255)

              aaCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccc                                                                     >  2:5808884/1‑92 (MQ=255)

               aCGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAccccct                                                                    >  1:2972927/1‑91 (MQ=255)

                  aaaTGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCtgtg                                                                 >  1:2608808/1‑93 (MQ=255)

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

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

                     tGGTTGATGACTATAATAGCCTGCGTTAAACTTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAg                                                              >  1:3629898/1‑93 (MQ=255)

                     tGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAg                                                              >  2:1840989/1‑93 (MQ=255)

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

                      ggTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGt                                                             >  2:712616/1‑93 (MQ=255)

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

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

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

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

                              aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                     >  2:2315354/1‑93 (MQ=255)

                              aCTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCAt                                                     >  1:6033210/1‑93 (MQ=255)

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

                               cTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTATCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATa                                                    <  1:1967225/93‑1 (MQ=255)

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

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

                                 ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  <  1:3803985/93‑1 (MQ=255)

                                 ataataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATaaa                                                  <  2:138217/93‑1 (MQ=255)

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

                                   aataGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATg                                                <  1:5097905/93‑1 (MQ=255)

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

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

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

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

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

                                        ccTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCa                                           >  1:35925/1‑93 (MQ=255)

                                         cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCat                                          >  2:314796/1‑93 (MQ=255)

                                         cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACACCCTGTGAAGTAACTCCATAAATGCCTCat                                          >  2:6185806/1‑93 (MQ=255)

                                         cTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCAACCCCTGTGAAGTAACTCCATAAATGCCTCat                                          >  2:1056045/1‑93 (MQ=255)

                                          tGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCata                                         >  2:3826029/1‑93 (MQ=255)

                                           gCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCatat                                        >  2:891385/1‑93 (MQ=255)

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

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

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

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

                                               tAAACGTTTTTGATGATAAATCCAGTGTCGGATGCGGC‑CCACCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTg                                    <  1:4302518/93‑1 (MQ=255)

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

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

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

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

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

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

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

                                                  aCGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTa                                 >  1:4062939/1‑93 (MQ=255)

                                                   cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                <  1:454271/93‑1 (MQ=255)

                                                   cGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTaa                                <  2:1590583/93‑1 (MQ=255)

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

                                                    gTTTTTGATGATAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAAc                               >  2:4311031/1‑93 (MQ=255)

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

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

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

                                                         tgatgaTAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGc                          <  2:2343433/93‑1 (MQ=255)

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

                                                           atgatAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaa                        >  1:831896/1‑93 (MQ=255)

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

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

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

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

                                                              atAAATGCAGTGTCGGATGCGGCGC‑ATCCGACACAGGTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCaataa                     >  2:5997634/1‑93 (MQ=255)

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

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

                                                                  aTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                 <  1:1180015/93‑1 (MQ=255)

                                                                  aTGCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAAtctc                 <  1:5826463/93‑1 (MQ=255)

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

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

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

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

                                                                    gCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               <  1:5631597/93‑1 (MQ=255)

                                                                    gCAGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccg               >  2:608225/1‑93 (MQ=255)

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

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

                                                                      aGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             >  1:132467/1‑93 (MQ=255)

                                                                      aGTGTCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTccgcc             >  2:590693/1‑93 (MQ=255)

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

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

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

                                                                       gtgtCGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACACGGCAATAATCTccgccg            >  2:247988/1‑93 (MQ=255)

                                                                        tgtCGGATGCGGCGC‑ATCCGAAACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCgg           >  1:1621020/1‑93 (MQ=255)

                                                                         ggcGGATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGa          <  1:4697761/91‑1 (MQ=255)

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

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

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

                                                                             gATGCGGCGC‑ATCCGACACAGTTCACCCCCTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCt      <  2:1927780/93‑1 (MQ=255)

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

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

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

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

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

CATTACGCGAGTACAACGAAATGGTTGATGACTATAATAGCCTGCGTTAAACGTTTTTGATGATAAATGCAGTGTCGGATGCGGCGCGATCCGACACAGTTCACCCCTTGTGAAGTAACTCCATAAATGCCTCATATCCTGTTAACCCGGCAATAATCTCCGCCGGATGCTGACA  >  NC_000913/4085588‑4085762

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