BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	2,819,255	A→G	5.9%	T471T (ACT→ACC)	alaS ←	alanyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	2,819,255	0	A	G	5.9%	277.9 / 6.3	101	T471T (ACT→ACC)	alaS	alanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand): ref base A (52/43); new base G (3/3); total (55/46)
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.23e-01

Only 100 of 103 total aligned reads displayed.

CAGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGATCATTG  >  W3110S.gb/2819189‑2819325

cAGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTc                                                                     <  1:1342915/70‑1 (MQ=255)

cAGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCa                                                                    <  1:2784808/71‑1 (MQ=255)

cAGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCa                                                                    <  1:2851566/71‑1 (MQ=255)

 aGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATTAACAAACAGCGCAGTCAc                                                                   <  1:52107/71‑1 (MQ=255)

 aGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCAc                                                                   <  1:220713/71‑1 (MQ=255)

 aGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCAc                                                                   <  1:1960749/71‑1 (MQ=255)

    aCGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACttt                                                                <  1:2493087/71‑1 (MQ=255)

          acaGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGtt                                                          >  1:2131547/1‑71 (MQ=255)

          acaGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGtt                                                          >  1:2842646/1‑71 (MQ=255)

           caGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGtt                                                          >  1:559499/1‑70 (MQ=255)

                  cTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCa                                                  <  1:3695548/71‑1 (MQ=255)

                  cTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCa                                                  <  1:1756217/71‑1 (MQ=255)

                    ggCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCa                                                        >  1:2122508/1‑63 (MQ=255)

                    ggCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCa                                                        >  1:86526/1‑63 (MQ=255)

                      ccTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATg                                              >  1:1736049/1‑71 (MQ=255)

                       cTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATgg                                             <  1:860364/71‑1 (MQ=255)

                        tGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGtt                                                     <  1:3849706/62‑1 (MQ=255)

                        tGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGAt                                               <  1:3049631/68‑1 (MQ=255)

                        tGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGAt                                               <  1:2093338/68‑1 (MQ=255)

                         gCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTc                                                         <  1:1005091/57‑1 (MQ=255)

                         gCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTc                                                         <  1:800472/57‑1 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCa                                                                    >  1:1483092/1‑45 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:370551/1‑65 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:375102/1‑65 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:1099868/1‑65 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATgg                                             <  1:2570973/68‑1 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTc                                           >  1:2240521/1‑70 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCa                                          >  1:2422276/1‑71 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCa                                          >  1:1403665/1‑71 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCa                                          >  1:3738941/1‑71 (MQ=255)

                          cATTGATGGCATCAACCGCTTTACCATCAACAAAAAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:1708241/1‑65 (MQ=255)

                            ttGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACttt                                                                <  1:2552196/47‑1 (MQ=255)

                            ttGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACttt                                                                <  1:2150094/47‑1 (MQ=255)

                            ttGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGt                                            <  1:286015/67‑1 (MQ=255)

                            ttGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGt                                            <  1:2500353/67‑1 (MQ=255)

                              gATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTcc                                                   >  1:2142506/1‑58 (MQ=255)

                              gATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAg                                       >  1:2605207/1‑70 (MQ=255)

                              gATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGc                                      >  1:3445066/1‑71 (MQ=255)

                              gATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGc                                      >  1:1005331/1‑71 (MQ=255)

                                tGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCa                                                  <  1:3619920/57‑1 (MQ=255)

                                 ggCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATg                                              >  1:1008168/1‑60 (MQ=255)

                                 ggCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATg                                              >  1:3625621/1‑60 (MQ=255)

                                 ggCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCtt                                   >  1:642128/1‑71 (MQ=255)

                                 ggCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCtt                                   >  1:2487639/1‑71 (MQ=255)

                                 ggCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCtt                                   >  1:3593925/1‑71 (MQ=255)

                                   cATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTc                                                    <  1:1622909/52‑1 (MQ=255)

                                   cATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTc                                                    <  1:2740965/52‑1 (MQ=255)

                                   cATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTc                                                    <  1:1106326/52‑1 (MQ=255)

                                   cATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTT‑‑AGATGGTCATAGCCTTTaa                                >  1:1343856/1‑70 (MQ=255)

                                     tCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                <  1:2223508/54‑1 (MQ=255)

                                      cAACCGCTTTACCATCAACAAACAGCGCGGTCACTTTGCCGTTCa                                                        >  1:2302196/1‑45 (MQ=255)

                                      cAACCGCTTTACCATCAACAAACAGCGCGGTCACTTTGCCGTTCa                                                        >  1:2865154/1‑45 (MQ=255)

                                      cAACCGCTTTACCATCAACAAACAGCGCGGTCACTTTGCCGTTCa                                                        >  1:2292039/1‑45 (MQ=255)

                                      cAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAAt                              >  1:381533/1‑71 (MQ=255)

                                      cAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAAt                              >  1:111843/1‑71 (MQ=255)

                                      cAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAAt                              >  1:192372/1‑71 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCTGTTCAGTTCCAGa                                                >  1:2060670/1‑52 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTc                                                    >  1:1489940/1‑48 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:104451/1‑52 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:3260680/1‑52 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCa                                          >  1:3266236/1‑58 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCAt                                         >  1:3404340/1‑59 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAAtt                             >  1:1006415/1‑71 (MQ=255)

                                       aaCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTAAGTTCCAGa                                                >  1:486275/1‑52 (MQ=255)

                                       aaCAGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGa                                                >  1:774372/1‑52 (MQ=255)

                                        aCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGtt                                                          <  1:1228800/41‑1 (MQ=255)

                                        aCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGtt                                                          <  1:1233294/41‑1 (MQ=255)

                                        aCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGtt                                                          <  1:554804/41‑1 (MQ=255)

                                        aCCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAAt                              >  1:1419621/1‑69 (MQ=255)

                                         ccGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCa                           <  1:3557778/71‑1 (MQ=255)

                                          cGCTTTACCATCAACAAACAGCGCGGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAg                          <  1:374388/71‑1 (MQ=255)

                                          cGCTTTACCATCAACAAACAGCGCGGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAg                          <  1:3318000/71‑1 (MQ=255)

                                          cGCTTTACCATCAACAAACAGCGCGGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAg                          <  1:3633665/71‑1 (MQ=255)

                                          cGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATg                                              <  1:2634474/51‑1 (MQ=255)

                                             tttACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGc                                      <  1:2245436/56‑1 (MQ=255)

                                              ttACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAg                                                       >  1:1154974/1‑38 (MQ=255)

                                              ttACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAg                                                       >  1:3476955/1‑38 (MQ=255)

                                               tACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTaaa                               <  1:470419/61‑1 (MQ=255)

                                               tACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCa                     >  1:3646593/1‑71 (MQ=255)

                                               tACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCa                     >  1:558587/1‑71 (MQ=255)

                                                aCCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGt                                                      <  1:2805036/37‑1 (MQ=255)

                                                   aTCAACAAACAGCTCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATa                                        <  1:158889/48‑1 (MQ=255)

                                                    tCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGt                 >  1:2734110/1‑70 (MQ=255)

                                                      aacaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGc                                      >  1:2007911/1‑47 (MQ=255)

                                                      aacaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCaa              <  1:246704/71‑1 (MQ=255)

                                                      aacaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCaa              >  1:668812/1‑71 (MQ=255)

                                                      aacaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCaa              >  1:2376383/1‑71 (MQ=255)

                                                      aacaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCaa              >  1:2906839/1‑71 (MQ=255)

                                                      aacaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCaa              <  1:3082503/71‑1 (MQ=255)

                                                       acaaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAac             <  1:1885544/71‑1 (MQ=255)

                                                         aaacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAacac           >  1:2579327/1‑71 (MQ=255)

                                                          aacaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGt                                            <  1:1985159/37‑1 (MQ=255)

                                                           acaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGAt                        >  1:3266954/1‑56 (MQ=255)

                                                           acaGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGAt                        >  1:17624/1‑56 (MQ=255)

                                                            caGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGa        <  1:3459207/71‑1 (MQ=255)

                                                             aGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCaa              <  1:2702911/64‑1 (MQ=255)

                                                              gcgcAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGATc      <  1:241070/71‑1 (MQ=255)

                                                              gcgcAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGATc      <  1:3068408/71‑1 (MQ=255)

                                                                  aGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGATCATTg  <  1:3145553/71‑1 (MQ=255)

                                                                  aGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGATCATTg  <  1:226352/71‑1 (MQ=255)

CAGCACGACCACAGCTTCCTGGCCTGCATTGATGGCATCAACCGCTTTACCATCAACAAACAGCGCAGTCACTTTGCCGTTCAGTTCCAGATGGTCATAGCCTTTAAATTCAGATGCACTGTCAACACGGATCATTG  >  W3110S.gb/2819189‑2819325

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