BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,829,233	C→A	L97M (CTG→ATG)	xanP →	xanthine:H(+) symporter XanP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,829,233	0	C	A	100.0%	172.3 / NA	42	L97M (CTG→ATG)	xanP	xanthine:H(+) symporter XanP
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (27/15); total (27/15)

TGGCATCGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGCTGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGCCCCGCTG  >  NC_000913/3829186‑3829277

tgGCATCGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGAtg                                            <  1:2226740/50‑1 (MQ=255)

  gCATCGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGAtgtt                                          <  1:1960842/50‑1 (MQ=255)

   cATCGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGAtgttg                                         >  1:1197626/1‑50 (MQ=255)

    aTCGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGAtgttgt                                        >  1:567538/1‑50 (MQ=255)

      cGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTc                                       >  1:3676091/1‑49 (MQ=255)

      cGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTc                                       >  1:4424689/1‑49 (MQ=255)

       gATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCt                                      <  1:4207247/49‑1 (MQ=255)

       gATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTa                                     >  1:3703182/1‑50 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    <  1:1157261/50‑1 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    <  1:336428/50‑1 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    >  1:3858987/1‑50 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    >  1:2280540/1‑50 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    >  1:1495414/1‑50 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    >  1:771856/1‑50 (MQ=255)

        attattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAt                                    <  1:1234517/50‑1 (MQ=255)

         ttattCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTAtt                                   >  1:3143915/1‑50 (MQ=255)

              cAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAggg                              >  1:1837333/1‑50 (MQ=255)

              cAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAgg                               <  1:430962/49‑1 (MQ=255)

               aaaTTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGc                             <  1:4234023/50‑1 (MQ=255)

               aaaTTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGc                             <  1:3781435/50‑1 (MQ=255)

                 aTTAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCAc                           >  1:737813/1‑50 (MQ=255)

                  ttAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCAc                           >  1:4264773/1‑49 (MQ=255)

                   tAAGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCAc                           >  1:4012938/1‑48 (MQ=255)

                    aaGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAg                        >  1:1193856/1‑50 (MQ=255)

                     aGGCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGc                       >  1:2918066/1‑50 (MQ=255)

                      ggCCTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCt                      >  1:3949896/1‑50 (MQ=255)

                        ccTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCtt                     >  1:4071232/1‑49 (MQ=255)

                         cTGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTc                    >  1:879823/1‑49 (MQ=255)

                          tGGGGTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCaa                  >  1:3310689/1‑50 (MQ=255)

                           ggggTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAAc                 >  1:4388072/1‑50 (MQ=255)

                            gggTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCaa                  <  1:1415741/48‑1 (MQ=255)

                            gggTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAAc                 <  1:1912066/49‑1 (MQ=255)

                             ggTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACtt               <  1:4197716/50‑1 (MQ=255)

                             ggTCCGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACtt               >  1:3489349/1‑50 (MQ=255)

                                ccGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTttgt            >  1:18871/1‑50 (MQ=255)

                                 cGGTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTttgt            <  1:2415476/49‑1 (MQ=255)

                                  ggTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTttgttg          >  1:2648020/1‑50 (MQ=255)

                                   gTTGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGc         >  1:1416830/1‑50 (MQ=255)

                                    ttGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGcc        <  1:4375904/50‑1 (MQ=255)

                                     tGGCTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGccc       >  1:2026136/1‑50 (MQ=255)

                                        cTCCGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGCCCCGc    >  1:2259886/1‑50 (MQ=255)

                                           cGGGATGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGCCCCGCTg  <  1:192677/49‑1 (MQ=255)

TGGCATCGATTATTCAAATTAAGGCCTGGGGTCCGGTTGGCTCCGGGCTGTTGTCTATTCAGGGCACCAGCTTCAACTTTGTTGCCCCGCTG  >  NC_000913/3829186‑3829277

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