BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	3,993,460	A→T	H16Q (CAT→CAA)	grxC ←	glutaredoxin 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	3,993,460	0	A	T	100.0%	158.4 / NA	42	H16Q (CAT→CAA)	grxC	glutaredoxin 3
Reads supporting (aligned to +/- strand): ref base A (0/0); new base T (20/22); total (20/22)

ACTCCTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGATGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACATTGGCC  >  NC_017635/3993414‑3993505

aCTCCTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTgg                                            <  1:747984/50‑1 (MQ=255)

aCTCCTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTgg                                            <  1:1920037/50‑1 (MQ=255)

 cTCCTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTgg                                            >  1:590383/1‑49 (MQ=255)

  tCCTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGc                                           >  1:380366/1‑49 (MQ=255)

   ccTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCaa                                         >  2:15560/1‑50 (MQ=255)

    cTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCa                                          <  1:274259/48‑1 (MQ=255)

      gAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATAc                                      >  2:1616566/1‑50 (MQ=255)

       aaaaCTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATa                                       <  1:79337/48‑1 (MQ=255)

       aaaaCTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACg                                     <  2:1054435/50‑1 (MQ=255)

       aaaaCTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACg                                     <  1:229289/50‑1 (MQ=255)

        aaaCTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACgg                                    >  1:655844/1‑50 (MQ=255)

        aaaCTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACgg                                    <  2:470855/50‑1 (MQ=255)

          aCTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGc                                  <  1:293390/50‑1 (MQ=255)

           cTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACggg                                   >  2:158203/1‑48 (MQ=255)

           cTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCa                                 <  1:1605848/50‑1 (MQ=255)

              aCGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAgg                               <  1:45234/49‑1 (MQ=255)

              aCGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAgg                               <  2:607671/49‑1 (MQ=255)

              aCGCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGt                              >  1:585842/1‑50 (MQ=255)

                gCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGttt                            <  2:1204838/50‑1 (MQ=255)

                gCCCTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGttt                            <  1:1354777/50‑1 (MQ=255)

                 cccTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGTTTc                           >  2:1056696/1‑50 (MQ=255)

                  ccTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCt                          >  1:1879604/1‑50 (MQ=255)

                  ccTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCt                          >  1:2042570/1‑50 (MQ=255)

                   cTTGCTGCTCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCtt                         <  2:118873/50‑1 (MQ=255)

                       ctgctCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTgg                      <  2:1738348/49‑1 (MQ=255)

                        tgctCAGCAGTGCTTTTGCACGTTGGCAATACGCGCAGGTTTCTTTGGta                    >  2:231694/1‑50 (MQ=255)

                           tCAGCAGTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGta                    >  2:900470/1‑47 (MQ=255)

                             agcagTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGa                <  2:1817905/49‑1 (MQ=255)

                             agcagTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGa                <  2:1325035/49‑1 (MQ=255)

                             agcagTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGAt               >  1:150174/1‑50 (MQ=255)

                              gcagTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGAtt              >  2:567481/1‑50 (MQ=255)

                              gcagTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGAtt              >  1:946774/1‑50 (MQ=255)

                               cagTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGAttt             >  1:157936/1‑50 (MQ=255)

                                agTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGAttt             >  1:5157/1‑49 (MQ=255)

                                agTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTc            <  2:1574500/50‑1 (MQ=255)

                                 gTGCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTc            >  1:753342/1‑49 (MQ=255)

                                   gCTTTTGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCaa          >  1:1659456/1‑49 (MQ=255)

                                     ttttGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACAt       >  2:1088626/1‑50 (MQ=255)

                                        tGCACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACATTg     <  1:1356187/49‑1 (MQ=255)

                                          cACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACATTGGcc  <  2:1734718/50‑1 (MQ=255)

                                          cACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACATTGGcc  <  2:1567785/50‑1 (MQ=255)

                                          cACGTTGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACATTGGcc  <  2:2065696/50‑1 (MQ=255)

ACTCCTGAAAACTCACGCCCTTGCTGCTCAGCAGTGCTTTTGCACGATGGCAATACGGGCAGGTTTCTTTGGTATAGATTTCAACATTGGCC  >  NC_017635/3993414‑3993505

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