BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,491,023	A→G	Y318H (TAT→CAT)	idnD ←	L‑idonate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,491,023	0	A	G	100.0%	757.8 / NA	235	Y318H (TAT→CAT)	idnD	L‑idonate dehydrogenase
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (125/110); total (125/110)

Only 100 of 245 total aligned reads displayed.

TATTTATCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGC  >  NC_000913/4490933‑4491112

tatttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGtt                                                                                         <  1:5118477/93‑1 (MQ=255)

 atttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTc                                                                                        >  1:250497/1‑93 (MQ=255)

  tttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCa                                                                                       <  2:2339620/93‑1 (MQ=255)

  tttatCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCa                                                                                       <  2:483619/93‑1 (MQ=255)

   ttatCTCTTAGAAAACAAGCTGGACTTTTACTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAg                                                                                      <  2:2422310/93‑1 (MQ=255)

     atCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCa                                                                                    <  1:5490337/93‑1 (MQ=255)

     atCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCa                                                                                    >  1:2486574/1‑93 (MQ=255)

       ctctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACt                                                                                  >  1:730343/1‑93 (MQ=255)

       ctctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACt                                                                                  >  1:901671/1‑93 (MQ=255)

        tctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTc                                                                                 <  1:4857362/93‑1 (MQ=255)

         ctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCa                                                                                >  1:1635162/1‑93 (MQ=255)

         ctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCa                                                                                <  1:2593198/93‑1 (MQ=255)

         ctTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCa                                                                                >  2:4270025/1‑93 (MQ=255)

          ttAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAg                                                                               >  2:1683297/1‑93 (MQ=255)

          ttAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAg                                                                               >  1:3222478/1‑93 (MQ=255)

           tAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGt                                                                              <  2:1139802/93‑1 (MQ=255)

           tAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGt                                                                              <  2:2350635/93‑1 (MQ=255)

            aGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTa                                                                             <  1:1094901/93‑1 (MQ=255)

            aGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTa                                                                             <  2:6085559/93‑1 (MQ=255)

              aaaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTaaa                                                                           >  2:1339719/1‑93 (MQ=255)

              aaaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTaaa                                                                           >  1:1550685/1‑93 (MQ=255)

               aaaCAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAg                                                                          >  2:4721272/1‑93 (MQ=255)

                  cAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCa                                                                       >  1:3623236/1‑93 (MQ=255)

                  cAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGCAAAGGCa                                                                       >  2:3875212/1‑93 (MQ=255)

                   aaGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAg                                                                      >  1:4969489/1‑93 (MQ=255)

                   aaGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAg                                                                      <  1:578396/93‑1 (MQ=255)

                   aaGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAg                                                                      <  2:1411774/93‑1 (MQ=255)

                   aaGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAg                                                                      <  2:2838110/93‑1 (MQ=255)

                   aaGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAg                                                                      <  1:4363114/93‑1 (MQ=255)

                    aGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGt                                                                     >  2:5849280/1‑93 (MQ=255)

                     gCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTg                                                                    >  2:603936/1‑93 (MQ=255)

                     gCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTg                                                                    <  2:4837573/93‑1 (MQ=255)

                     gCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTg                                                                    <  2:4514920/93‑1 (MQ=255)

                      cTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTgg                                                                   <  1:5829417/93‑1 (MQ=255)

                      cTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTgg                                                                   >  2:4764650/1‑93 (MQ=255)

                        ggACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGAt                                                                 <  1:3449384/93‑1 (MQ=255)

                        ggACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGAt                                                                 >  2:4059140/1‑93 (MQ=255)

                        ggACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGAt                                                                 <  1:1003045/93‑1 (MQ=255)

                             tttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATa                                                            <  1:6210792/93‑1 (MQ=255)

                             tttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATa                                                            >  1:2223987/1‑93 (MQ=255)

                              ttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATaa                                                           <  1:4728924/93‑1 (MQ=255)

                              ttGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATaa                                                           <  1:5623776/93‑1 (MQ=255)

                               tgctgcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAAc                                                          <  2:195431/93‑1 (MQ=255)

                                 ctgcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGc                                                        >  1:5370338/1‑93 (MQ=255)

                                  tgcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGcc                                                       >  1:4072720/1‑93 (MQ=255)

                                  tgcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGcc                                                       >  1:3790306/1‑93 (MQ=255)

                                   gcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCa                                                      >  1:5330509/1‑93 (MQ=255)

                                   gcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCa                                                      >  1:5930400/1‑93 (MQ=255)

                                   gcCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCa                                                      >  2:4502810/1‑93 (MQ=255)

                                     cTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCAtt                                                    >  2:651851/1‑93 (MQ=255)

                                      tGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTc                                                   >  1:621135/1‑93 (MQ=255)

                                      tGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTc                                                   >  2:3525090/1‑93 (MQ=255)

                                      tGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTc                                                   <  1:4584575/93‑1 (MQ=255)

                                      tGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTc                                                   >  1:3958264/1‑93 (MQ=255)

                                       gggTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCg                                                  >  1:42726/1‑93 (MQ=255)

                                            ttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccagc                                             <  2:5153727/93‑1 (MQ=255)

                                            ttATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccagc                                             >  2:678868/1‑93 (MQ=255)

                                             tATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCgccagcc                                            <  1:3208136/93‑1 (MQ=255)

                                                cACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATg                                         <  1:836272/93‑1 (MQ=255)

                                                 aCCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGa                                        >  2:499562/1‑93 (MQ=255)

                                                 aCCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGa                                        >  2:37075/1‑93 (MQ=255)

                                                   cGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGaca                                      >  2:2530448/1‑93 (MQ=255)

                                                   cGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGaca                                      >  2:3326230/1‑93 (MQ=255)

                                                   cGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGaca                                      <  1:4913475/93‑1 (MQ=255)

                                                    ggCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGacac                                     <  1:1957246/93‑1 (MQ=255)

                                                     gCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACt                                    <  2:763693/93‑1 (MQ=255)

                                                     gCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACt                                    <  1:3951564/93‑1 (MQ=255)

                                                       gAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAGGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGc                                  <  2:261105/93‑1 (MQ=255)

                                                        aaaCGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCg                                 >  2:2038252/1‑93 (MQ=255)

                                                         aaCGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCgg                                >  1:3893623/1‑93 (MQ=255)

                                                         aaCGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCgg                                >  1:4177997/1‑93 (MQ=255)

                                                          aCGTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGt                               <  2:883188/93‑1 (MQ=255)

                                                            gTAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTAt                             >  1:141316/1‑93 (MQ=255)

                                                             tAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTAtt                            <  1:5639361/93‑1 (MQ=255)

                                                             tAGCGCCTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTAtt                            >  1:4035005/1‑93 (MQ=255)

                                                               gcgcCTCTTCCAGGTCAGTGATGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTaa                          >  2:3599466/1‑93 (MQ=255)

                                                                 gcCTCTTCAAGGTCAGTGAAGGTAGCTTCAGCACCAAGTAAAGGCAGTGGATTAAAAACGCTATTCGCCAGCCATGACACTGCAGTATTAAAt                        <  1:2382180/93‑1 (MQ=255)

                                                                  ccTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       <  2:5787959/93‑1 (MQ=255)

                                                                  ccTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       >  1:6079557/1‑93 (MQ=255)

                                                                  ccTCTTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       >  1:3544798/1‑93 (MQ=255)

                                                                  ccTCTTCCAGGGCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAAtt                       <  2:5694027/93‑1 (MQ=255)

                                                                   ctctTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTc                      <  1:2899428/93‑1 (MQ=255)

                                                                    tctTCCAGGTCAGTGCAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCg                     >  2:3770865/1‑93 (MQ=255)

                                                                    tctTCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCg                     >  1:517100/1‑93 (MQ=255)

                                                                       tCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTg                  <  1:114534/93‑1 (MQ=255)

                                                                       tCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTg                  <  2:3372258/93‑1 (MQ=255)

                                                                       tCCAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTg                  >  2:4365482/1‑93 (MQ=255)

                                                                        ccAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTgg                 >  1:3118134/1‑93 (MQ=255)

                                                                         cAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGt                >  2:1848648/1‑93 (MQ=255)

                                                                         cAGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGt                >  1:2706365/1‑93 (MQ=255)

                                                                          aGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTa               >  2:5332833/1‑93 (MQ=255)

                                                                          aGGTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTa               >  2:3001235/1‑93 (MQ=255)

                                                                            gTCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTaaa             >  2:3986102/1‑93 (MQ=255)

                                                                             tCAGTGAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTaaaa            >  2:5214939/1‑93 (MQ=255)

                                                                                  gAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGaa       >  2:5671777/1‑93 (MQ=255)

                                                                                  gAAGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGaa       <  1:48369/93‑1 (MQ=255)

                                                                                   aaGGGATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGaaa      <  1:2032743/93‑1 (MQ=255)

                                                                                     gggATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAaga    >  2:4621537/1‑93 (MQ=255)

                                                                                      ggATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAagag   <  2:5587223/93‑1 (MQ=255)

                                                                                       gATGTTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGc  >  2:3175743/1‑93 (MQ=255)

TATTTATCTCTTAGAAAACAAGCTGGACTTTTGCTGCCTGGGTTTTATCACCGGCGAAACGTAGCGCCTCTTCCAGGTCAGTGAAGGGATATTCAGCACTCAGTAAAGGCAGTGGATTGATAACGCCATTCGCCAGCCATGACACTGCGGTATTAAATTCGCTGGTAAAACGGAAAGAGC  >  NC_000913/4490933‑4491112

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