BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	753,654	G→A	A272V (GCG→GTG)	gltA ←	citrate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	753,654	0	G	A	95.9%	154.7 / ‑2.5	49	A272V (GCG→GTG)	gltA	citrate synthase
Reads supporting (aligned to +/- strand): ref base G (2/0); new base A (22/25); total (24/25)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

AGCCCATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGC  >  NC_000913/753559‑753745

aGCCCATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGcaccgc                                                                                         <  2:1160172/100‑1 (MQ=255)

  cccATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCtt                                                                                       <  2:1224481/100‑1 (MQ=255)

   ccATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTc                                                                                      >  1:535890/1‑100 (MQ=255)

      tCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGtt                                                                                   <  2:1264970/100‑1 (MQ=255)

       cAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTa                                                                                  <  1:365204/100‑1 (MQ=255)

        aGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAg                                                                                 <  1:1109527/100‑1 (MQ=255)

         ggCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGc                                                                                <  1:1114369/100‑1 (MQ=255)

         ggCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGc                                                                                <  1:18608/100‑1 (MQ=255)

         ggCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCAGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCACTTCGTTAGc                                                                                >  2:612167/1‑100 (MQ=255)

                aGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCAccgccg                                                                         <  1:1150413/100‑1 (MQ=255)

                 gAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGt                                                                        >  1:286014/1‑100 (MQ=255)

                    tCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCACCGTgcg                                                                     >  1:227145/1‑100 (MQ=255)

                    tCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTgcg                                                                     <  1:629571/100‑1 (MQ=255)

                     cATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTgcgc                                                                    <  1:591888/100‑1 (MQ=255)

                        ttttGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAgg                                                                 <  2:884090/100‑1 (MQ=255)

                        ttttGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAgg                                                                 <  1:677600/100‑1 (MQ=255)

                               tttCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCcac                                                          >  2:218675/1‑100 (MQ=255)

                                 tCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAg                                                        <  2:658746/100‑1 (MQ=255)

                                  ngCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGt                                                       >  1:244908/2‑100 (MQ=255)

                                    cACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGa                                                     >  1:1159213/1‑100 (MQ=255)

                                      cgacgaACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAg                                                   >  1:1303199/1‑100 (MQ=255)

                                       gacgaACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGc                                                  <  2:580500/100‑1 (MQ=255)

                                         cgaACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCaa                                                >  1:1089758/1‑100 (MQ=255)

                                           aaCAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATa                                              >  2:720649/1‑100 (MQ=255)

                                           aaCAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGAAACGCTTCGTTAGCACAGGCGTGCGCAGGTCCCCACAGTGAAGCAATa                                              >  2:1037369/1‑100 (MQ=255)

                                             cAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATAcc                                            <  1:803583/100‑1 (MQ=255)

                                                aTTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACctgc                                         <  2:84267/100‑1 (MQ=255)

                                                   ccGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCAGCGATTCGTTAGCAGCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCtgctgc                                      >  2:196388/1‑100 (MQ=255)

                                                   ccGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCtgctgc                                      <  1:849281/100‑1 (MQ=255)

                                                    cGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCg                                     <  1:607506/100‑1 (MQ=255)

                                                    cGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCg                                     >  1:395710/1‑100 (MQ=255)

                                                      gAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGAt                                   <  2:535886/100‑1 (MQ=255)

                                                          gtgtTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAg                               <  2:5229/100‑1 (MQ=255)

                                                            gtTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGc                             >  2:852919/1‑100 (MQ=255)

                                                            gtTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGc                             >  1:295076/1‑100 (MQ=255)

                                                               tAACGGAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCaaa                          >  1:83086/1‑100 (MQ=255)

                                                                aaCGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCacac                         >  1:519482/1‑100 (MQ=255)

                                                                   ggAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACggg                      <  2:286011/100‑1 (MQ=255)

                                                                   ggAGCTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACggg                      >  2:1072567/1‑100 (MQ=255)

                                                                       cTGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCg                  <  1:858203/100‑1 (MQ=255)

                                                                        tGATTTCTTCCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGc                 >  1:568302/1‑100 (MQ=255)

                                                                            ttcttcCAGCATTTTCGGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCAccc             >  1:442637/1‑100 (MQ=255)

                                                                              cttcCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGa           >  2:866745/1‑100 (MQ=255)

                                                                               ttcCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGaa          <  1:1078542/100‑1 (MQ=255)

                                                                                tcCAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAag         <  2:912851/100‑1 (MQ=255)

                                                                                 ccAGCA‑TTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAagag       >  1:218187/1‑100 (MQ=255)

                                                                                  cggCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAagag       >  1:884383/3‑100 (MQ=255)

                                                                                  cAGCATTTTCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAagag       >  2:359644/1‑100 (MQ=255)

                                                                                    gCAGTTGCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGcc     <  2:1122749/100‑1 (MQ=255)

                                                                                       tgttCAGCACCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGc  <  1:407175/98‑1 (MQ=255)

AGCCCATCAGGCGGAAAGAATCATTTTTGTCTTTCGCACGACGAACAAATTCCGGAATGTGTTTAACGGAGCTGATTTCTTCCAGCATTTTCAGCGCCGCTTCGTTAGCACCGCCGTGCGCAGGTCCCCACAGTGAAGCAATACCTGCTGCGATACAGGCAAACGGGTTCGCACCCGAAGAGCCAGC  >  NC_000913/753559‑753745

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