BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,225,139	T→C	intergenic (+49/‑171)	metH → / → yjbB	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,225,139	0	T	C	100.0%	393.9 / NA	126	intergenic (+49/‑171)	metH/yjbB	B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein
Reads supporting (aligned to +/- strand): ref base T (0/0); new base C (67/59); total (67/59)

Only 100 of 162 total aligned reads displayed.

AGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4225050‑4225228

ggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGcc                                                                                        <  1:1081569/92‑1 (MQ=255)

ggcgcTGGCTGGCACCGAATCTGGGGCATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGcc                                                                                        >  2:2129138/2‑93 (MQ=255)

 gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGccc                                                                                       <  1:5713551/93‑1 (MQ=255)

 gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGccc                                                                                       >  2:3153860/1‑93 (MQ=255)

 gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTCCCTT‑ACAACAAACAGG‑GCGCTAAATGAGCGccc                                                                                       <  1:558871/93‑1 (MQ=255)

  ctctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCt                                                                                      <  1:5988469/91‑1 (MQ=255)

  cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCt                                                                                      <  1:3987316/93‑1 (MQ=255)

   gctggctgGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTg                                                                                     <  1:4155285/93‑1 (MQ=255)

    ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGt                                                                                    <  2:5715395/93‑1 (MQ=255)

    ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCCGt                                                                                    >  2:2028037/1‑93 (MQ=255)

    agagctggcACCGAATCGGGGTTATGACGGGGATTGATCCAAAAATCTGTCATTTTCCCTT‑ACAACAAACAGG‑GCGCAAAATGAGCGCCCTGt                                                                                    <  2:3545193/90‑1 (MQ=255)

     tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGGCGCTCTCAATGGGCCCCCTgc                                                                                    >  1:1476652/1‑92 (MQ=255)

     tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGtc                                                                                   >  2:3106270/1‑93 (MQ=255)

       gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGtctc                                                                                 >  2:2972464/1‑93 (MQ=255)

         tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCtt                                                                               >  2:5022499/1‑93 (MQ=255)

          ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCttt                                                                              <  1:4989006/93‑1 (MQ=255)

          ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCttt                                                                              <  2:1032179/93‑1 (MQ=255)

          ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGACCTGTCTCttt                                                                              >  1:3843229/1‑93 (MQ=255)

            cACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTttat                                                                            >  2:1073679/1‑93 (MQ=255)

             aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACGGG‑GCGCTCAATGAGCGCCCTGTCTCTttatt                                                                           <  1:3163346/93‑1 (MQ=255)

             aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCGCTGTCTCTttatt                                                                           >  2:2022166/1‑93 (MQ=255)

                gAATCTGGGGTATGACGCGGACTGATTCACAAATCGGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAat                                                                        <  1:4632607/93‑1 (MQ=255)

                  aTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAatat                                                                      <  1:575522/93‑1 (MQ=255)

                   tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATg                                                                     <  1:2928461/93‑1 (MQ=255)

                    cTGGGGTATGACGCGGACTAATTCAAAAATCTGTCACTTTCCCTT‑ACAACAAACAGG‑GCCCTCAATGAGCGCCCTGTCTCTTTATTAATATGa                                                                    <  1:3672141/93‑1 (MQ=255)

                       gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                 >  1:2794605/1‑93 (MQ=255)

                        ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATTAGCGCCCTGTCTCTTTATTAATATTAAacc                                                                >  1:3591349/1‑92 (MQ=255)

                          tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACt                                                              <  1:5578086/93‑1 (MQ=255)

                            tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                            <  1:2043298/93‑1 (MQ=255)

                            tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                            <  2:4271257/93‑1 (MQ=255)

                             gaagcgGAGTGATTCAAAAATCTGTCACTTTCCCTT‑ACAACAAACAGG‑GCGCTAAATGAACGCCCTGTCTCTTTATTAATATGAAACACTtat                                                           <  1:4946697/90‑1 (MQ=255)

                             gACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtat                                                           >  2:5816754/1‑93 (MQ=255)

                             gACGCGGACTGATCCACAAATCTGTCATTTTTCCTT‑ACAACAAACAGG‑GCCCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtat                                                           <  1:4824334/93‑1 (MQ=255)

                              gcgcgGCCTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACGGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                          <  1:2613227/92‑1 (MQ=255)

                              aCGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                          <  2:789800/93‑1 (MQ=255)

                              aCGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                          <  2:3005962/93‑1 (MQ=255)

                                gcgTACTGATTCAAAAATCTGTC‑TTTTTCCTTAAAAACAAACGGG‑GCGCTAAAGGCGCGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                        <  1:5378603/93‑1 (MQ=255)

                                   gACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGGCCTGTCTCTTTATTAATATGAAACACTTATACTTGa                                                     >  2:519372/1‑93 (MQ=255)

                                   gACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GAGCGCCATGTGAGCCCCCGCCCTTTTTTTATATTAAAAACACATACTTGa                                                     >  1:2240054/1‑93 (MQ=255)

                                     cTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGAGCACTGAGGGCCCTGCGCCTTTGTTTATATTAAACACTTTTACTCGata                                                   >  1:1206780/1‑91 (MQ=255)

                                      tGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCCATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaaa                                                  >  2:3603649/1‑92 (MQ=255)

                                      tGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaaa                                                  >  2:1378426/1‑92 (MQ=255)

                                      tGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                  >  2:615457/1‑93 (MQ=255)

                                      tGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                  >  2:1047282/1‑93 (MQ=255)

                                       tatccacaAATCTTTCTCTTTTCCTT‑ACAACAAACAGG‑GCGCTCAAGGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                 <  1:212401/89‑1 (MQ=255)

                                       gATTCACAACTCTTTCACTTTTCCTT‑CCAACAACCAGG‑CCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                 <  1:1676983/93‑1 (MQ=255)

                                       gATTCACAAATCTGTCATTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                 <  1:443376/93‑1 (MQ=255)

                                        aTTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                                <  2:878133/93‑1 (MQ=255)

                                         ttCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTTATATTAAACACTTATACTTGAAACaag                                               >  2:2433329/1‑91 (MQ=255)

                                         ttCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                               <  2:2781442/93‑1 (MQ=255)

                                          tCACAAATCTGTCATTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                              <  1:1413047/93‑1 (MQ=255)

                                           cacaAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCt                                             <  2:3097269/93‑1 (MQ=255)

                                             caAATCTGTCACTTTTCCTT‑ACAACAAACGGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           <  2:5199090/93‑1 (MQ=255)

                                             caAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           >  2:4239143/1‑93 (MQ=255)

                                             caAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           >  2:3861163/1‑93 (MQ=255)

                                             caAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                           <  2:579793/93‑1 (MQ=255)

                                               aaTCTGTCTTTTTTCCTA‑ACAACAAACAGG‑GCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                         <  2:1357640/93‑1 (MQ=255)

                                                aTCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                        >  1:3133800/1‑93 (MQ=255)

                                                 tCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                       >  1:3918980/1‑93 (MQ=255)

                                                   tGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                     >  1:3167922/1‑93 (MQ=255)

                                                    ttCATTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                    <  2:2632986/92‑1 (MQ=255)

                                                    gTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                    >  2:3200161/1‑93 (MQ=255)

                                                       ctttttCCTT‑ACAACAAACAGG‑CCCCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                 <  2:1857935/91‑1 (MQ=255)

                                                       aCTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                 >  1:6239348/1‑93 (MQ=255)

                                                       aCTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                                 >  1:1888010/1‑93 (MQ=255)

                                                        cTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                                <  1:2921026/93‑1 (MQ=255)

                                                         ttttCTTA‑ACAACAAGCGGG‑CCGTTAAAGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                               <  2:5602901/93‑1 (MQ=255)

                                                         ttttCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGGAACAAGCTGGGATAATCTTAGc                               >  2:1664051/1‑93 (MQ=255)

                                                         ttttCCTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                               >  2:635840/1‑93 (MQ=255)

                                                           ttCCTA‑AAAAAAAACAGG‑GCGCTAAATGAGCCCCCTGTCTCTTTATTAATTTGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                             <  1:3372855/93‑1 (MQ=255)

                                                             ccTT‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGa                           >  1:3644474/1‑93 (MQ=255)

                                                               tt‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                         >  2:828444/1‑93 (MQ=255)

                                                               tt‑ACAACAAACAGG‑GCGCTCAATGAGCGCCCCGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAAACTTAGCCCGaaa                         >  2:5366816/1‑93 (MQ=255)

                                                                t‑ACAAAAAACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaa                        >  1:2770985/1‑93 (MQ=255)

                                                                   caacaaACAGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                      >  2:2840474/1‑93 (MQ=255)

                                                                    aacaaacgGG‑GCGTTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                     <  1:3544171/93‑1 (MQ=255)

                                                                    aacaaacaGG‑GCGCTCAATGAGCGCCCTTTCTCTTTTTTAATATTAAACACTTATACTGGAAACAGGCTTGAATAATCTTAGCCGGaaaaaaa                     >  2:4852915/1‑91 (MQ=255)

                                                                     aaaaacaGG‑GCGCTAAATGACCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                    <  2:371625/91‑1 (MQ=255)

                                                                        aacaGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 >  1:2120543/1‑93 (MQ=255)

                                                                        aacaGG‑GCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 <  1:4484353/93‑1 (MQ=255)

                                                                         agaGG‑CCGATCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                <  2:1241903/91‑1 (MQ=255)

                                                                         acggg‑gCGCAAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                <  2:3452107/90‑1 (MQ=255)

                                                                         acaGG‑GCGCTAAATGGGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt                <  1:3760737/93‑1 (MQ=255)

                                                                          caGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               <  2:6240782/93‑1 (MQ=255)

                                                                          caGG‑CCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc               <  1:5039759/93‑1 (MQ=255)

                                                                           aGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCGGAAAAAAGAAAGTCa              >  2:5417253/1‑93 (MQ=255)

                                                                           aGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa              >  2:5002671/1‑93 (MQ=255)

                                                                           aGG‑GCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa              <  2:1551250/93‑1 (MQ=255)

                                                                            gg‑gCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa             >  2:3690091/1‑93 (MQ=255)

                                                                               gcgcTCAATGAGCGCTCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa           <  1:2999108/93‑1 (MQ=255)

                                                                               gcgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa           <  1:3003903/93‑1 (MQ=255)

                                                                                  cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:5085616/1‑93 (MQ=255)

                                                                                  cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:5514547/1‑93 (MQ=255)

                                                                                  cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:4779674/1‑93 (MQ=255)

                                                                                     aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     <  2:2156281/93‑1 (MQ=255)

                                                                                     aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     >  1:6261598/1‑93 (MQ=255)

                                                                                      aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg    >  2:3331860/1‑93 (MQ=255)

                                                                                      aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg    <  2:5981192/93‑1 (MQ=255)

                                                                                      aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg    <  1:408254/93‑1 (MQ=255)

                                                                                       tGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt   <  2:479584/93‑1 (MQ=255)

AGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAACAGG‑GCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4225050‑4225228

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