BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,861,751	T→C	78.9%	N423S (AAC→AGC)	lldP ←	L‑lactate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,861,751	0	T	C	78.9%	155.6 / 35.6	71	N423S (AAC→AGC)	lldP	L‑lactate permease
Reads supporting (aligned to +/- strand): ref base T (5/10); new base C (29/27); total (34/37)
Fisher's exact test for biased strand distribution p-value = 2.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTCAGC  >  W3110S.gb/3861687‑3861816

gaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCg                                                                >  1:3065575/1‑68 (MQ=255)

gaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGa                                                               >  1:1554083/1‑69 (MQ=255)

gaagGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGa                                                               >  1:1700126/1‑69 (MQ=255)

   ggTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGATGATGACAGTCCGGAATAGCTCGaaa                                                             >  1:1904176/1‑68 (MQ=255)

         tGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGc                                                      <  1:343010/69‑1 (MQ=255)

         tGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGc                                                      <  1:1332272/69‑1 (MQ=255)

          gCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                     <  1:3573458/69‑1 (MQ=255)

          gCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCg                                                     <  1:2527509/69‑1 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAgg                                                       <  1:2924588/66‑1 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAgg                                                       <  1:2340699/66‑1 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAgg                                                       <  1:3079096/66‑1 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCg                                                     <  1:884635/68‑1 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGa                                                    >  1:2858565/1‑69 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGa                                                    >  1:1308405/1‑69 (MQ=255)

           cATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGa                                                    >  1:2121686/1‑69 (MQ=255)

             tGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAAt                                                  <  1:3696706/69‑1 (MQ=255)

               aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAAt                                                  >  1:3608829/1‑67 (MQ=255)

               aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAAt                                                  >  1:1984792/1‑67 (MQ=255)

               aCCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAAt                                                  >  1:1231013/1‑67 (MQ=255)

                  ggTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCCAATGCCAg                                             <  1:2634946/69‑1 (MQ=255)

                     gtgCGCCAGTGCCAGCGCCAGTGTTGATGAAAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCAc                                          <  1:2870946/69‑1 (MQ=255)

                       gcgcCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATaa                                                          <  1:1465375/51‑1 (MQ=255)

                            aGTGCCAGCGCCAGTGTTGATGACAGTCTGGAATAGTTCGAAATaa                                                          <  1:1382530/46‑1 (MQ=255)

                               gCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAAt                                                            <  1:1946869/41‑1 (MQ=255)

                               gCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCg                                                     <  1:558822/48‑1 (MQ=255)

                               gCCAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCg                                                     <  1:2594497/48‑1 (MQ=255)

                                ccAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               <  1:2633021/69‑1 (MQ=255)

                                ccAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               <  1:2442971/69‑1 (MQ=255)

                                ccAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               >  1:3123975/1‑69 (MQ=255)

                                ccAGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               >  1:1210162/1‑69 (MQ=255)

                                  aGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATaaa                                                         >  1:3714322/1‑41 (MQ=255)

                                  aGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               >  1:1918636/1‑67 (MQ=255)

                                  aGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               >  1:3435016/1‑67 (MQ=255)

                                  aGCGCCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATgg                               >  1:766373/1‑67 (MQ=255)

                                    cgcCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTa                           <  1:348296/69‑1 (MQ=255)

                                    cgcCAGTGTTGATGACAGTCCGGAATAGCTCGAAAT‑AAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:2260954/69‑1 (MQ=255)

                                     gcCAGTGTTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:2516065/69‑1 (MQ=255)

                                           gtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                                >  1:2994268/1‑57 (MQ=255)

                                           gtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                                >  1:1035841/1‑57 (MQ=255)

                                           gtTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATCCCAGCACCAt                                       >  1:2104285/1‑50 (MQ=255)

                                           gtTGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAg                                             <  1:2218259/44‑1 (MQ=255)

                                           gtTGATGACAGTCCGGAATAGCTCCAAATAAAGGCGAATGCCAg                                             <  1:3616266/44‑1 (MQ=38)

                                            ttGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAg                                             <  1:2213623/43‑1 (MQ=255)

                                            ttGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:2319478/62‑1 (MQ=255)

                                            ttGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:2291147/62‑1 (MQ=255)

                                            ttGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:3040138/62‑1 (MQ=255)

                                            ttGATGACAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:742096/62‑1 (MQ=255)

                                              gatgaCAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                                >  1:1811347/1‑54 (MQ=255)

                                                 gaCAGTCCGGAATAGCTCGAAATAAAGGCGAATGCCAGCAc                                          >  1:3170318/1‑41 (MQ=255)

                                                    aGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGt           >  1:683755/1‑69 (MQ=255)

                                                       ccGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAg                          <  1:170870/51‑1 (MQ=255)

                                                       ccGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGtt          >  1:2324924/1‑67 (MQ=255)

                                                       ccGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGtt          >  1:2320589/1‑67 (MQ=255)

                                                       ccGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGtt          >  1:537668/1‑67 (MQ=255)

                                                       ccGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGtt          >  1:864118/1‑67 (MQ=255)

                                                       ccGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGt           >  1:179348/1‑66 (MQ=255)

                                                        cGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCtt       <  1:1722987/69‑1 (MQ=255)

                                                        cGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                                <  1:89378/44‑1 (MQ=255)

                                                        cGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                                <  1:3652133/44‑1 (MQ=255)

                                                        cGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATg                                <  1:1230764/44‑1 (MQ=255)

                                                        cGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTa                           <  1:2061889/49‑1 (MQ=255)

                                                        cGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTa                           <  1:2057610/49‑1 (MQ=255)

                                                        cGGAATAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTa                           <  1:173927/49‑1 (MQ=255)

                                                           aaTAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGcc              >  1:2399606/1‑59 (MQ=255)

                                                           aaTAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGcc              >  1:1644165/1‑59 (MQ=255)

                                                           aaTAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGcc              >  1:1281663/1‑59 (MQ=255)

                                                           aaTAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGc               >  1:2671740/1‑58 (MQ=255)

                                                             tAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTGAGc  >  1:246947/1‑69 (MQ=255)

                                                             tAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTCAGc  >  1:3098598/1‑69 (MQ=255)

                                                             tAGCTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTCAGc  >  1:2491973/1‑69 (MQ=255)

                                                               gTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAg                             <  1:618741/40‑1 (MQ=255)

GAAGGTGAATGCATGACCGGTGTGCGCCAGTGCCAGCGCCAGTGTTGATGACAGTCCGGAATAGTTCGAAATAAAGGCGAATGCCAGCACCATACCGATGGAGTAGATGGGCAGAGCCAGTTCTTTCAGC  >  W3110S.gb/3861687‑3861816

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