BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,481,820	G→A	13.2%	A673V (GCG→GTG)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,481,820	0	G	A	13.2%	68.4 / 3.0	37	A673V (GCG→GTG)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand): ref base G (18/14); new base A (3/2); total (21/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TCAGCTCGAGATACCAGTCACAGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGA  >  NC_000913/4481680‑4481955

gcaggccGAGGTACCAGTCACAGAACGGGTTCCAGGTGAACTCAGACAGAATGCCGGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGcc                                                                                                                                       <  1:572236/137‑1 (MQ=255)

  aGCTCGAGATACCAGTCACAGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCGGCGGTCCGCCAg                                                                                                                                     >  2:402212/1‑143 (MQ=255)

   gCTCGAGATACCAGTCACAGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAGGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGATGAGCCCCGCCAGAATCCAGCGGTCCGCCAGc                                                                                                                                    >  1:112546/1‑143 (MQ=255)

    cTCGAGATACCAGTCACAGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTCCGCTTTGATGGTCGCGTTGACCTCCTCCAGAATCCCGCCGTccgcccgcc                                                                                                                                   >  2:849392/1‑142 (MQ=255)

    cTCGAGATACCAGTCACAGAAATGGTTCCAGGTGAACTCATACAGAATGACTGCAGCGATAGCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCGGACACAATACAGCGGTCCGcccact                                                                                                                                   >  2:868379/1‑139 (MQ=255)

         gagaCCAGGCACAGAACCGGTGCCAGGTGAACTCAGACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGc                                                                                                                              <  1:506058/140‑1 (MQ=255)

            ggcAGTCACAGAACGGGTTCCAGGGGAGCTGTGGCAGGAGGCCTGCGGCGATATGGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCACCAGCGACAGCGTc                                                                                                                           <  2:356808/141‑1 (MQ=255)

            aCCAGTCACAGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCACCAGCGACAGCGTc                                                                                                                           >  2:259291/1‑143 (MQ=255)

                   acaGTACTGGTTCCGGGTGATCTCATACAGGATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCCGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTcgc                                                                                                                    <  2:76721/143‑1 (MQ=255)

                     aGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTcgccg                                                                                                                  >  1:179401/1‑143 (MQ=255)

                     aGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTcgccg                                                                                                                  >  2:406350/1‑143 (MQ=255)

                          tGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTACAGCGCTTCGCGGTACGATTTGATGGACTGGTTGAACTCCGCCAGAATCCAGCGGTCCACCCGCGACAGCGTCACTTCCCCGACGtt                                                                                                             >  2:281598/1‑143 (MQ=255)

                                   gTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCTGCCGTTGAGGCCGCa                                                                                                    >  2:535140/1‑143 (MQ=255)

                                       aCTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCCTGTTGAACTCCGCCAGACTCCCGCGGTCCGCCAGCGACACCGTCCTTTCGCCGCCCTTGAATCCGccaacc                                                                                                >  2:63019/1‑138 (MQ=255)

                                         tCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTg                                                                                              <  1:566484/143‑1 (MQ=255)

                                                 aaTGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGGATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACctgctg                                                                                      >  1:734420/1‑143 (MQ=255)

                                                        gcggcgTGATCGGAGCGGAAGCTGGCCGTCGCGTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTtcat                                                                               <  2:274367/143‑1 (MQ=255)

                                                        gcggcgATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGGCATTTCGCCGGCGTTGAAGCCGCAATGCTGACCTTCTGCGTtcct                                                                               >  2:369557/1‑141 (MQ=255)

                                                          tgcgAGATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGCCTGGTTGGACTCCGCCACAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTtcatca                                                                             <  2:706942/142‑1 (MQ=255)

                                                           gcgATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAg                                                                            <  1:544915/143‑1 (MQ=255)

                                                                   gAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGACCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGCCCTCCTGTGTCCCTCAGcacacccc                                                                    >  1:407923/1‑139 (MQ=255)

                                                                   gAAGCGGAAGCTGTCCAGCGCTTCACAGTACGCTAAGATAAACAAGTTGAACTCCGCCAAAATCCAGCCGACAACCAGAGACACCATCAACTAGCCGACGCTGAAGCAGCAATACCCACCCTCTGTAAGCATGAGCCCAAAGc                                                                    >  2:620017/1‑143 (MQ=255)

                                                                    aagcggaagcggTGCAGTGTTTCGCGGTACCGTATGATGGTGTGGTTGAACTCCGCCAGAAGCCAGCGGTCCGCCAGCGACAGCGTCATGTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAGGCg                                                                   <  1:779713/143‑1 (MQ=255)

                                                                    aaGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCg                                                                   <  1:305470/143‑1 (MQ=255)

                                                                             cTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAACGACAGCGTCCTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTCCTGCGTTCATCAGCCCAAAGCGGCTGGCGtt                                                          >  2:468802/1‑143 (MQ=255)

                                                                                  cAGCGCTTCGCGGTACGCTTTGATGGTCAGGTTGAACTCCGCCAAAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCATCTGAGTTCAACAGCACAAAGCGGCCAGCGATCcaca                                                     >  2:369228/1‑143 (MQ=255)

                                                                                   aGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGCTCCGCCAGCGACAGCGTCCTCTCGCCGCCGTTCAAGCCGCAATCCCGACCTTCTGCGTTCCTCAGCACACAGCCGCCGGCGTTCcacgg                                                    >  2:617452/1‑141 (MQ=255)

                                                                                        ttCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGt                                               <  1:40938/143‑1 (MQ=255)

                                                                                          cgcgGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGACTGGCGTTCCACAGCTTGTTa                                             <  1:406350/143‑1 (MQ=255)

                                                                                               tACGCTTTGATGGTCTTGTTGTACTCCGCGAGAATCGAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCTTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGaa                                        <  2:697547/143‑1 (MQ=255)

                                                                                                       gATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCACCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                                <  1:259291/143‑1 (MQ=255)

                                                                                                        aTGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCCATCCTGACCTTCTGTGTTCATCCGCACAACGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTa                               >  1:527037/1‑143 (MQ=255)

                                                                                                                gTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                       >  1:700223/1‑143 (MQ=255)

                                                                                                                          gCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCat             >  1:792645/1‑143 (MQ=255)

                                                                                                                                aTCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAACTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCa       <  1:791890/143‑1 (MQ=255)

                                                                                                                                 tCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAg      >  1:258479/1‑143 (MQ=255)

                                                                                                                                   cAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCAGCAATCCTGACATTCTGTGTTCATCAGCACGAGGCGGCTGGCGTTCCACAGCGTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATAGCCCGGtg    >  1:817903/1‑142 (MQ=255)

                                                                                                                                     gCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCGGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa  <  1:627518/143‑1 (MQ=255)

TCAGCTCGAGATACCAGTCACAGAACTGGTTCCAGGTGAACTCATACAGAATGCCTGCGGCGATATCGAAGCGGAAGCTGTCCAGCGCTTCGCGGTACGCTTTGATGGTCTGGTTGAACTCCGCCAGAATCCAGCGGTCCGCCAGCGACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGA  >  NC_000913/4481680‑4481955

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