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%	350.4 / NA	114	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 (56/58); total (56/58)

Only 100 of 147 total aligned reads displayed.

GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225051‑4225226

gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGccc                                                                                     <  1:257898/93‑1 (MQ=255)

gcgcTGGCTGGCACCGAATCGGGGGTATGACGCGGACTGATTCACAAATCTGTCATTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGccc                                                                                     <  2:501722/93‑1 (MQ=255)

 cgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCCCTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCt                                                                                    <  2:4509422/93‑1 (MQ=255)

    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGtc                                                                                 <  1:2236429/93‑1 (MQ=255)

    tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGtc                                                                                 <  1:4561444/93‑1 (MQ=255)

       ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCCCAATGAGCGCCCTGTctct                                                                              <  2:4578607/93‑1 (MQ=255)

        tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                             >  2:3781631/1‑93 (MQ=255)

         ggcACGGAATCTGGGGTATGACGCGGACTGATCCACAAATCTGTTATTTTTCCTA‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCttt                                                                            <  1:1104435/93‑1 (MQ=255)

         gccgccGAATCGGGGTTATGCCGCGGACTGATTCACAAATCTGTCACTTTCCCTA‑ACAACAAA‑CAGGGCGCTTAATGAGCGCCCTGTCTCttt                                                                            <  2:4311963/89‑1 (MQ=255)

           cACCGAATGTGGGGTATGACGCGGACTGATTCACAAATCTGTCATTTTTCCTT‑ACAACAAA‑CAGGGCGCTAAATGAGCGCCCTGTCTCTttat                                                                          <  1:1050529/93‑1 (MQ=255)

               gAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGAGCGCTCAATGAGCGCCCTTTCTCTTTATTAat                                                                      >  2:4126306/1‑93 (MQ=255)

                attCTGGGTTATGACGCGGACTAATTCACAATTCTGTCACTTTTCATA‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAata                                                                     <  2:2084794/91‑1 (MQ=255)

                 atgtGGGGAATGCCGCGGATTAATTCAAAATTCTTTTTCTTTCCCTA‑AAAACAAC‑CGGGGCGTTAAAGGACCCCCCTGTCTCTTTATTAatat                                                                    <  2:3282509/90‑1 (MQ=255)

                  tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATg                                                                   >  1:3634977/1‑93 (MQ=255)

                   cTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGa                                                                  <  2:1446516/93‑1 (MQ=255)

                     ggggTATGCCGCGGACTATTTAACAATTTTGTC‑TTTTTCCTTAACAAAAAA‑CGGGGCCCTCAATGAGCCCCCTGTCTCTTTATTAATATGaaa                                                                <  2:122863/93‑1 (MQ=255)

                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaaa                                                                >  1:3101449/1‑93 (MQ=255)

                     ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAACGCCCTGTCTCTTTATTAATATGaaa                                                                >  2:3491946/1‑93 (MQ=255)

                      gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                               <  1:4535136/93‑1 (MQ=255)

                       ggTATGACGCGGACTGATTCACAATTCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGGGCGCCCTGTCTCTTTATTAATATGAAaca                                                              <  2:1426410/93‑1 (MQ=255)

                       ggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAaca                                                              >  2:1479946/1‑93 (MQ=255)

                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAacac                                                             <  1:716032/93‑1 (MQ=255)

                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAacac                                                             <  1:1311201/93‑1 (MQ=255)

                        gTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTAAATGGGCCCCCTGTCTCTTTATTAATATGAAacac                                                             <  1:1391324/93‑1 (MQ=255)

                         tATGCCGGGAACTTATTCAAAATTCTGTCATTTTTCCTA‑AAAACAAA‑CGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACt                                                            <  2:2544560/93‑1 (MQ=255)

                         tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACt                                                            >  2:3695053/1‑93 (MQ=255)

                         tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑AAGGGGGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACt                                                            >  1:2608682/1‑93 (MQ=255)

                          aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACtt                                                           >  2:2906634/1‑93 (MQ=255)

                           tgccgcgGACTGATTCACAAATCTTTCACTTTTCCTT‑ACAACAAG‑CGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTta                                                          <  1:2015448/90‑1 (MQ=255)

                           tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                          <  1:1428806/93‑1 (MQ=255)

                           tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta                                                          <  1:2162159/93‑1 (MQ=255)

                            gACGCGTACTAATTCACAAATCTGTCATTTTTCCTT‑ACAACAAC‑CAGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTtat                                                         <  2:2754163/93‑1 (MQ=255)

                              cgcgGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                       <  2:4529066/93‑1 (MQ=255)

                              cgcgGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATAc                                                       >  2:739355/1‑93 (MQ=255)

                               gcgGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                      >  1:1515067/1‑93 (MQ=255)

                                cgGACTTATCCACAAATCTGTCTTTTTTCTTT‑ACAACAAA‑CAGGGCGCGAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                     <  1:419720/93‑1 (MQ=255)

                                cgGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                     <  2:2789490/93‑1 (MQ=255)

                                  gACTGATTCACAAATCTGTCCCTTTTCCTT‑ACAAAAAA‑CGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                   <  1:3466601/93‑1 (MQ=255)

                                  gACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTTATATGAAACACTTATACTTGa                                                   >  2:322931/1‑93 (MQ=255)

                                  gACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                   >  2:2294864/1‑93 (MQ=255)

                                    cTGATCCACAAATCTGTCATTTTCCCTT‑ACAACAAG‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                 <  2:4462897/93‑1 (MQ=255)

                                    atGATTCACAAATCTGTCACTTTCCCTT‑ACAACAAG‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGaaa                                                 <  1:3405068/92‑1 (MQ=255)

                                      taTTCACAAGTCGTTCATTTTCCTTT‑ACAAAAAA‑CAGGCCGCTGAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                               <  1:2840827/92‑1 (MQ=255)

                                      gATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                               >  1:4448749/1‑93 (MQ=255)

                                      gATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTACACTGGAAACa                                               >  1:1863535/1‑93 (MQ=255)

                                      aaTTCACAAATCTGTCTCTTTTCCTT‑ACAACAAA‑CGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                               <  2:3463681/92‑1 (MQ=255)

                                       aTTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAg                                              <  2:3787309/93‑1 (MQ=255)

                                        ttCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                             >  1:2912904/1‑93 (MQ=255)

                                        ttCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCCGTCTCTTCTTTAATATGAAAAACATATAATTGAAACaag                                             >  1:1648429/1‑91 (MQ=255)

                                         tCACAAATTTGACACTTTTGCAT‑GCATCAAA‑CGGGGCCCGCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                            <  1:541848/93‑1 (MQ=255)

                                         tCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGGTCCATGAGCGCCCTGTCTCTTTATTTATATGAAACACTTATACTTGAAACAggg                                            >  2:2696683/1‑92 (MQ=255)

                                         tCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGAGCTCAATGAGCGCCCTGTCTCTTTATTAATATTAAACACTTATTATTGAAACAggg                                            >  1:3564047/1‑92 (MQ=255)

                                          cacaAATCTGTCTTTTTTCCTT‑ACAACAAA‑CAGGGCGCTAAAGGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCt                                           <  2:321518/93‑1 (MQ=255)

                                          cacaAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCt                                           <  1:2201062/93‑1 (MQ=255)

                                          cacaAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCt                                           >  2:4646877/1‑93 (MQ=255)

                                            caAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                         >  1:328498/1‑93 (MQ=255)

                                             aaaTCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGa                                        >  1:2837501/1‑93 (MQ=255)

                                                tCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                     >  1:4522990/1‑93 (MQ=255)

                                                  tGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                   >  2:1976313/1‑93 (MQ=255)

                                                  tGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCCGTCTCTTTTTTAATATTTAAAACTTATACTTGAAAAAAGCTGGgataat                                   >  2:4009818/1‑93 (MQ=255)

                                                   gTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGgataatt                                  >  2:2089079/1‑92 (MQ=255)

                                                   gTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGagtaaat                                  >  1:353640/1‑91 (MQ=255)

                                                    tCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATTAGCGCCCCGTCTCTTTTTTTATATTAAACACTTATACTCGAAACAAGCTGGAACAATCt                                 >  2:2637818/1‑93 (MQ=255)

                                                    tCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                 <  1:1462417/93‑1 (MQ=255)

                                                     cACTTTTCCTA‑ACAACAAA‑CGGGGCGCACAATGCCCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCtt                                <  1:1641258/93‑1 (MQ=255)

                                                      cattttCCTT‑ACAAAAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               <  1:714097/91‑1 (MQ=255)

                                                      aCTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               <  2:768079/93‑1 (MQ=255)

                                                       cTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                              <  1:2355863/93‑1 (MQ=255)

                                                       cTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                              <  1:1466787/93‑1 (MQ=255)

                                                       cTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                              >  1:3890942/1‑93 (MQ=255)

                                                       cTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                              <  1:405426/93‑1 (MQ=255)

                                                        ttttTCTT‑AAAAAAAACCGGGGCTCTAAATGCGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAg                              <  2:1642272/93‑1 (MQ=255)

                                                        ttttCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGCAACAGGCTGGAATAATCTTAGc                             >  2:2983110/1‑93 (MQ=255)

                                                        ttttCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTatc                             >  1:320961/1‑91 (MQ=255)

                                                        ttttCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             <  1:3778988/93‑1 (MQ=255)

                                                        tttcCCTA‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             <  1:1735376/93‑1 (MQ=255)

                                                          ttCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                           >  1:293779/1‑93 (MQ=255)

                                                          ttCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCg                           <  2:3034301/93‑1 (MQ=255)

                                                             cat‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAATCCTTATACTGGAAACAGGCTGGAAATATCTTAGGCggga                        >  2:4618455/3‑91 (MQ=255)

                                                              tt‑ACAACAAG‑CAGGGCGCTCAAGGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                       <  1:4247472/93‑1 (MQ=255)

                                                               t‑ACAACAAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTTGAAACAGGCTGGAATAATCTTAGCCGGaaaa                      >  2:3377284/1‑93 (MQ=255)

                                                                 acaacaAA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                     >  2:1376949/1‑93 (MQ=255)

                                                                 aaaacaaa‑caGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                     <  1:1043247/91‑1 (MQ=255)

                                                                  caacaaA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                    >  1:1359430/1‑93 (MQ=255)

                                                                  caacaaA‑CAGGGCGCTCAATGAGCGCCCTGTCTCTTTAATAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                    <  1:3085563/93‑1 (MQ=255)

                                                                  aaacaaa‑caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                    <  1:3761349/92‑1 (MQ=255)

                                                                   aacaaa‑caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAca                   >  1:1486118/1‑92 (MQ=255)

                                                                      aaa‑caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                >  1:547784/1‑93 (MQ=255)

                                                                       aa‑caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg               <  2:2317191/93‑1 (MQ=255)

                                                                       aa‑caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg               <  1:1793582/93‑1 (MQ=255)

                                                                        a‑caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt              >  1:126767/1‑93 (MQ=255)

                                                                          cggggCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc             <  2:427238/91‑1 (MQ=255)

                                                                             ggCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGCAATAATCTTAGCCGGAAAAACGAATGTCaaa          <  2:640595/93‑1 (MQ=255)

                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGGAAAAAGAATGTCaaaaag       >  2:2623390/1‑91 (MQ=255)

                                                                                 cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgga      >  1:1295502/1‑92 (MQ=255)

                                                                                  tCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGca     >  2:3698082/1‑93 (MQ=255)

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

                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGATAAACGAATGTCAAAAGGCACAg  >  1:605100/1‑93 (MQ=255)

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

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

GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTT‑ACAACAAA‑CAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225051‑4225226

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