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%	680.2 / NA	216	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 (117/99); total (117/99)

Only 100 of 322 total aligned reads displayed.

GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4225049‑4225228

gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgcgc                                                                                          <  2:6035865/91‑1 (MQ=255)

gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgcgc                                                                                          >  2:2775574/3‑93 (MQ=255)

     ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCATTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGt                                                                                     <  2:9401987/93‑1 (MQ=255)

     ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGt                                                                                     >  2:1124376/1‑93 (MQ=255)

     ctggctggCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGt                                                                                     >  1:2791295/1‑93 (MQ=255)

      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                    <  1:8394042/93‑1 (MQ=255)

      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                    >  2:9514735/1‑93 (MQ=255)

       ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGAGCGCTCAATGAGCGGCCTGtct                                                                                   >  2:808432/1‑93 (MQ=255)

          tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                                >  2:3159234/1‑93 (MQ=255)

           ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCttt                                                                               >  1:1312110/1‑93 (MQ=255)

            gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCCCTGTCTCtttt                                                                              >  2:5507809/1‑92 (MQ=255)

               ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGGGCTCAATCAGCGCCCTGTCTCttttttt                                                                           >  2:7471375/1‑92 (MQ=255)

               ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGGGCCCTGTCTCTttatta                                                                           >  2:2621339/1‑93 (MQ=255)

               ccGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatta                                                                           >  2:1660123/1‑93 (MQ=255)

                  aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAata                                                                        <  1:6701496/93‑1 (MQ=255)

                    tcgggggTATGACGCGGACTGATTCAAAAATCTGTCTCTTTTCCTTACAAAAAACGGGGCGCTCAAGAAGCTCCCTGTCTCTTTATTAATATg                                                                      <  2:922650/90‑1 (MQ=255)

                    tCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGGGGTCAATGAACCGCCTGTCTCTTTTTTTatata                                                                      >  2:239350/1‑92 (MQ=255)

                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                    <  2:5268866/93‑1 (MQ=255)

                       ggtgTATGACGCGGTCTGTTTCACAATTCTGTCTCTTTCCCTCACAACAAACAGGGCGCTCAAGGAGCCCCCTGTCTCTTTATTAATATGaaa                                                                   <  1:5659370/90‑1 (MQ=255)

                            aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACtt                                                              >  1:1679310/1‑93 (MQ=255)

                            aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACtt                                                              <  2:3442779/93‑1 (MQ=255)

                            aTGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACtt                                                              >  1:5352939/1‑93 (MQ=255)

                             tGACTCGGACTGATTCAAAAATCTGACGCTTTTCCTTACAAGAAACAGGGCTCTAAATAATGGCCCTGTCTCTTTATTAATATGAAACACTta                                                             <  2:2910904/93‑1 (MQ=255)

                             tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCGCAATGAGCGCCCTGTTTTTTTATTAATATGAAACACTta                                                             >  1:478684/1‑93 (MQ=255)

                              gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta‑t                                                           <  1:1056678/93‑1 (MQ=255)

                              gACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta‑t                                                           <  1:972364/93‑1 (MQ=255)

                               gcgcgGTGTGATTCACAAACCGATCACTTTCCCTGCCAACAAACGGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta‑ta                                                          <  2:9766341/92‑1 (MQ=255)

                               aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTta‑ta                                                          <  1:4635922/93‑1 (MQ=255)

                                 gcgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACt                                                        <  1:6275271/93‑1 (MQ=255)

                                  cgGACTGATTCACAAATTTGTCTTTTTTCTTAAAAAAAAACGGGGCGCTAAATGGGCCCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTg                                                       <  2:196527/93‑1 (MQ=255)

                                   ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGGGCTCAATGAGCGCCCCGTCCCTTTATTTATATTAAACACTTA‑TACTgg                                                      >  1:7401432/1‑93 (MQ=255)

                                   ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTgg                                                      >  1:5037276/1‑93 (MQ=255)

                                    gACTGATTCCCAAATCTGTCCCTTTCCCTTACAACAACCAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGa                                                     <  2:8311446/93‑1 (MQ=255)

                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGa                                                     >  1:9335055/1‑93 (MQ=255)

                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCCCAATGAGCGCCCCGTCTCTTTTTTTATTTTAAACACTTA‑TTCTTGa                                                     >  2:3895001/1‑93 (MQ=255)

                                      cTGATTCACAATCCTCTCCCTTTCCCTAACAACAAACAGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGaaa                                                   <  2:4303493/93‑1 (MQ=255)

                                       tGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAAc                                                  >  2:8210198/1‑93 (MQ=255)

                                       tGATTCACAAATCTGTCACTTTTCCTTAAAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAAc                                                  <  1:2647702/93‑1 (MQ=255)

                                        gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACa                                                 <  1:8517750/93‑1 (MQ=255)

                                        gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACa                                                 >  2:8769894/1‑93 (MQ=255)

                                        gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACa                                                 <  2:7014834/93‑1 (MQ=255)

                                         aTTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGGGCTCAATTAGTGTCCTGTCTCTTTTTTTATGTTAAAAACCTATTCCTGGaaaag                                                 >  2:5539177/1‑91 (MQ=255)

                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCAAGAGCGCCCTGTCTCTTTATTTATATGAAACAATTA‑TACTTGaaaaaag                                               >  1:1920566/1‑89 (MQ=255)

                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAgg                                               >  2:4423265/1‑93 (MQ=255)

                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAgg                                               >  1:2997445/1‑93 (MQ=255)

                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAgg                                               >  2:6430678/1‑93 (MQ=255)

                                            tgaaaaTCTGGCAGTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCt                                             <  2:760514/90‑1 (MQ=255)

                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATTAGCGCCCCGTCTCTTTTTTTATATTAAAAACTTA‑TACTTGAAACAAGCTgg                                           >  2:3805496/1‑93 (MQ=255)

                                                aaTCTGTCATTTTTCCTTACAACAAACAGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGaa                                         <  2:4037608/93‑1 (MQ=255)

                                                aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGaa                                         <  2:7920271/93‑1 (MQ=255)

                                                 agcTGTCACTTTTCCTTACAACAAACAGGGGCCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTCA‑TACTGGAAACAGGCTGGaat                                        <  1:1819844/91‑1 (MQ=255)

                                                  tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGaata                                       <  1:817536/93‑1 (MQ=255)

                                                    tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGaataat                                     >  1:758534/1‑93 (MQ=255)

                                                    tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGaataat                                     <  1:8261956/93‑1 (MQ=255)

                                                      ttatttttCGTTACAATAAAAAGGCCGGTAAATGAGCCCACTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCt                                   <  1:4107185/89‑1 (MQ=255)

                                                       cACTTTTCATTCCAACAAACAGGCCGCGCAGTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCtt                                  <  1:3896053/93‑1 (MQ=255)

                                                        ctttttCCTCACAACAACCAGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTa                                 <  1:4231736/91‑1 (MQ=255)

                                                        cttttcCCTAAAAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTa                                 <  2:473466/91‑1 (MQ=255)

                                                         cTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAg                                <  2:6840303/93‑1 (MQ=255)

                                                          ttttCCTAAAAACAAACGGGCCGCTCAAGGGCCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGc                               <  1:7214827/93‑1 (MQ=255)

                                                           tttCCTTAAAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGcc                              <  1:2695547/93‑1 (MQ=255)

                                                           ttcccTTACAACAACCGGGCCGCTCATTGAGCCCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGcc                              <  1:5912687/90‑1 (MQ=255)

                                                             tctttAAAACAGACAGGGCGCTCAATGAGCGCTCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCgg                            <  1:4110198/90‑1 (MQ=255)

                                                               cTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                          <  1:6513102/93‑1 (MQ=255)

                                                               cTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                          <  2:6109903/93‑1 (MQ=255)

                                                               cTTACAACAAACAGGGCGCTAAAGGGGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                          <  2:7748144/93‑1 (MQ=255)

                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGGATAATCTTAGCCAGaaa                         >  2:9624426/1‑93 (MQ=255)

                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                         <  1:3539665/93‑1 (MQ=255)

                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                         >  1:3109185/1‑93 (MQ=255)

                                                                 tACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAAAAGGCTGGAATAATATTAGCCGGaaaa                        >  1:2501900/1‑93 (MQ=255)

                                                                  tcaacaaACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACCCTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                       <  2:8711904/92‑1 (MQ=255)

                                                                  acaacaAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                       >  1:3282257/1‑93 (MQ=255)

                                                                  acaacaAACAGGGCGCTAAATGAGCTCACTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                       <  2:3547043/93‑1 (MQ=255)

                                                                   cagaaaacaGGGCACAAAATGAGCGCGCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                      <  1:465545/89‑1 (MQ=255)

                                                                   aaaaaaacgGGCCGCTAAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAc                      <  2:2239325/89‑1 (MQ=255)

                                                                    aacaaacaGGGCCCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                     <  2:9819593/93‑1 (MQ=255)

                                                                      caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   <  1:9501069/93‑1 (MQ=255)

                                                                      aaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                   <  1:9201915/92‑1 (MQ=255)

                                                                       aaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                  >  2:1734777/1‑93 (MQ=255)

                                                                        accggggCGCTAAAGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 <  1:9138682/89‑1 (MQ=255)

                                                                        aacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 <  1:3619839/93‑1 (MQ=255)

                                                                        aacaGGGCGCTAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg                 <  1:8951935/93‑1 (MQ=255)

                                                                            gggCGTTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa             <  2:1245914/93‑1 (MQ=255)

                                                                             ggCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCCGGAAAACGAATGTCaaa            >  2:902363/1‑93 (MQ=255)

                                                                             ggCGCTAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaa            <  2:5791045/93‑1 (MQ=255)

                                                                               cgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAg          >  1:5181847/1‑93 (MQ=255)

                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg         >  2:5917641/1‑93 (MQ=255)

                                                                                 cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:1397268/1‑93 (MQ=255)

                                                                                 cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:5384890/1‑93 (MQ=255)

                                                                                 cTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGc        >  2:9303107/1‑93 (MQ=255)

                                                                                    aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     <  2:656275/93‑1 (MQ=255)

                                                                                    aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     <  1:6879845/93‑1 (MQ=255)

                                                                                    aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     >  1:2392214/1‑93 (MQ=255)

                                                                                    aaTGAGCGCCCTGTCTCTTTATTAATATGAAACACCTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     >  2:8475043/1‑93 (MQ=255)

                                                                                    aaTGAGCCCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGcaca     <  2:5186114/93‑1 (MQ=255)

                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg    <  2:6706356/93‑1 (MQ=255)

                                                                                      ggAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGTAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGt   <  1:8662380/92‑1 (MQ=255)

                                                                                       gAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa  >  1:9898891/1‑93 (MQ=255)

                                                                                       gAGCGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa  >  1:4352915/1‑93 (MQ=255)

                                                                                       gAGCGCCCTGCCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTa  <  1:3403024/93‑1 (MQ=255)

GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTA‑TACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTA  >  NC_000913/4225049‑4225228

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