BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,448,541	+C	intergenic (+46/‑92)	ytfQ → / → ytfR	putative LACI‑type transcriptional regulator/YtfR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,448,541	1	.	C	100.0%	617.8 / NA	208	intergenic (+46/‑92)	ytfQ/ytfR	putative LACI‑type transcriptional regulator/YtfR
Reads supporting (aligned to +/- strand): ref base . (0/0); new base C (135/73); total (135/73)

Only 100 of 230 total aligned reads displayed.

ATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑G‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448449‑4448623

aTACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCg                                                                                       >  1:957090/1‑93 (MQ=255)

 tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                      >  2:6166208/1‑92 (MQ=255)

 tACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgc                                                                                      >  1:4895766/1‑92 (MQ=255)

  aCCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcg                                                                                     <  2:3563671/93‑3 (MQ=255)

   ccGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCgcgc                                                                                    <  2:2756755/93‑4 (MQ=255)

             agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGc                                                                         >  1:3502656/1‑93 (MQ=37)

             agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGc                                                                         >  1:5121479/1‑93 (MQ=37)

             agaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGc                                                                         >  1:3623063/1‑93 (MQ=37)

              gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCg                                                                        >  1:4730919/1‑93 (MQ=37)

              gaaTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCg                                                                        >  2:5919096/1‑93 (MQ=37)

              gaaTTAGAGAAGAAGAAAAATATGGGGAATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCg                                                                        >  1:1407616/1‑93 (MQ=25)

                aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGAt                                                                      >  2:1532589/1‑93 (MQ=37)

                aTTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGAt                                                                      <  2:5306169/93‑1 (MQ=37)

                  tAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATgg                                                                    <  1:1855677/93‑1 (MQ=37)

                   agagAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGt                                                                   >  2:2805107/1‑93 (MQ=37)

                    gagaAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTg                                                                  >  2:4934388/1‑93 (MQ=37)

                     agaagaagaAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGc                                                                 >  1:5829941/1‑93 (MQ=37)

                     agaagaagaAAAATAGGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGc                                                                 <  2:2413753/93‑1 (MQ=25)

                      gaagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCa                                                                >  2:4089406/1‑93 (MQ=37)

                       aagaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCAc                                                               >  1:2750468/1‑93 (MQ=37)

                        agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACt                                                              >  2:5214200/1‑93 (MQ=37)

                        agaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACt                                                              >  1:3163196/1‑93 (MQ=37)

                         gaagaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTg                                                             >  1:345035/1‑93 (MQ=37)

                           agaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGaa                                                           >  2:3467055/1‑93 (MQ=37)

                            gaaAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAg                                                          >  1:370480/1‑93 (MQ=37)

                             aaaaaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGc                                                         >  2:5166983/1‑93 (MQ=37)

                              aaaaTATGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGcc                                                        <  2:4569262/93‑1 (MQ=25)

                                aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      >  2:609402/1‑93 (MQ=37)

                                aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      <  2:2077199/93‑1 (MQ=37)

                                aaTATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTg                                                      >  2:5739251/1‑93 (MQ=37)

                                 atatGGGGTATGGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAGGGCGCTCGC‑GCCGCATCAGGCGATGGTGCACTGAAGCCTGa                                                     <  1:5015409/93‑1 (MQ=25)

                                  tatGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGAt                                                    <  2:2542589/93‑1 (MQ=37)

                                   atGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATg                                                   <  1:205331/93‑1 (MQ=37)

                                    tGGTGTATTGAGGGTTGCTCTGCCTGATGCCGATCCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGc                                                  <  2:1432397/93‑1 (MQ=21)

                                        gTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACg                                              >  1:3026522/1‑93 (MQ=25)

                                         tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             <  2:1889237/93‑1 (MQ=25)

                                         tATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGc                                             >  1:6131709/1‑93 (MQ=25)

                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            >  1:172718/1‑93 (MQ=25)

                                          aTTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCt                                            <  2:1754872/93‑1 (MQ=25)

                                           ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           <  1:5642075/93‑1 (MQ=25)

                                           ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           >  2:3131807/1‑93 (MQ=25)

                                           ttGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCtt                                           >  2:4184617/1‑93 (MQ=25)

                                           tcgAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCCCTGAAGCCTGATGCGACGCtt                                           <  1:3257650/91‑1 (MQ=21)

                                            tGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTa                                          <  1:4174649/93‑1 (MQ=25)

                                             gAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAc                                         >  2:4162660/1‑93 (MQ=25)

                                              aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGTATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        <  2:3752201/93‑1 (MQ=25)

                                              aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGGGACGCGTagc                                        >  1:5099783/1‑91 (MQ=21)

                                              aGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTAcc                                        <  1:4644018/93‑1 (MQ=25)

                                               gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  2:2447480/1‑93 (MQ=25)

                                               gggTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcg                                       >  1:2701183/1‑93 (MQ=25)

                                                ggTTGCTATGCCTGATGCCGATTAGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgc                                      >  2:2825740/1‑93 (MQ=25)

                                                 gTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACcgcg                                     >  1:2513496/1‑93 (MQ=25)

                                                  ttGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGt                                    <  2:774639/93‑1 (MQ=25)

                                                    gCTATGCCTGATGCCGATTCGTAGGCTGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCt                                  >  2:6394915/1‑93 (MQ=21)

                                                     cTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCtt                                 >  1:5922490/1‑93 (MQ=25)

                                                      tATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                >  2:5028851/1‑93 (MQ=25)

                                                      tATCCCTGATGCCCATTCGTAGGCCGGATAAGGCCCTCGCCC‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTa                                <  2:3594948/93‑1 (MQ=21)

                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                               >  2:6313554/1‑93 (MQ=21)

                                                       aTGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTAt                               >  2:2005146/1‑93 (MQ=21)

                                                        tGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATc                              >  1:222789/1‑93 (MQ=21)

                                                         gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  1:4325677/1‑93 (MQ=21)

                                                         gCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCa                             >  2:3315186/1‑93 (MQ=21)

                                                          ccTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCAt                            <  2:3684017/93‑1 (MQ=21)

                                                           cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:1318856/1‑93 (MQ=21)

                                                           cTGATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATg                           >  1:5474899/1‑93 (MQ=21)

                                                            tGATGCCGATTCGTAGGCGGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGc                          <  1:83956/93‑1 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:1022002/93‑1 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         >  2:3022217/1‑93 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:6010751/93‑1 (MQ=21)

                                                             gATGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGcc                         <  1:4357177/93‑1 (MQ=21)

                                                              aTGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCt                        <  2:4664065/93‑1 (MQ=21)

                                                               tGCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTa                       >  2:1281497/1‑93 (MQ=21)

                                                                gCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      <  2:487827/93‑1 (MQ=21)

                                                                gCCGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTAc                      <  1:909479/93‑1 (MQ=21)

                                                                 ccGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACt                     <  1:3276984/93‑1 (MQ=21)

                                                                  cGATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTg                    <  1:6213706/93‑1 (MQ=21)

                                                                   gATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   >  1:4547515/1‑93 (MQ=21)

                                                                   gATTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTgg                   >  2:5143517/1‑93 (MQ=21)

                                                                    aTTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  <  2:1699976/93‑1 (MQ=18)

                                                                    aTTCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTggg                  >  1:4146983/1‑93 (MQ=18)

                                                                     ttCGTGGGCATGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  2:3750036/93‑1 (MQ=18)

                                                                     ttCGTAGGCGGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  1:5573361/93‑1 (MQ=21)

                                                                     ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 <  2:218188/93‑1 (MQ=21)

                                                                     ttCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGa                 >  1:513088/1‑93 (MQ=21)

                                                                      tCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                <  2:1497504/93‑1 (MQ=21)

                                                                      tCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                <  1:6054489/93‑1 (MQ=21)

                                                                      tCGTAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAg                <  2:911382/93‑1 (MQ=21)

                                                                        gTAGCCCGGATAAGGCCCCCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCa              <  2:4807328/93‑1 (MQ=21)

                                                                         tAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             >  1:1236580/1‑93 (MQ=21)

                                                                         tAGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCAc             >  2:5401021/1‑93 (MQ=21)

                                                                          aGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACg            >  2:4768402/1‑93 (MQ=21)

                                                                          aGGCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACg            >  1:19335/1‑93 (MQ=21)

                                                                           ggCCGGATAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGc           >  2:1147612/1‑93 (MQ=21)

                                                                                 aTAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     >  1:6108957/1‑93 (MQ=25)

                                                                                 aTAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTaca     >  1:2972150/1‑93 (MQ=25)

                                                                                  tAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    >  2:3197327/1‑93 (MQ=25)

                                                                                  tAAGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    <  1:1266728/93‑1 (MQ=25)

                                                                                  aaaGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTacac    <  1:2420245/92‑1 (MQ=25)

                                                                                   aaGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACAcc   >  1:3668747/1‑93 (MQ=25)

                                                                                    aGGCGCTCGC‑GCCGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCg  >  1:402262/1‑93 (MQ=25)

ATACCGCAAAAGAAGAATTAGAGAAGAAGAAAAATATGGGGTATTGAGGGTTGCTATGCCTGATGCCGATTCGTAGGCCGGATAAGGCGCTCG‑‑G‑CGCATCCGGCGATGGTGCACTGAAGCCTGATGCGACGCTTACCGCGTCTTATCATGCCTACTGGGAGCACGCTTTACACCG  >  NC_000913/4448449‑4448623

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