BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	AP009048	2,929,500	A→C	8.7%	K423N (AAA→AAC)	sdaB →	L‑serine deaminase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AP009048	2,929,500	0	A	C	8.7%	316.3 / 27.0	127	K423N (AAA→AAC)	sdaB	L‑serine deaminase II
Reads supporting (aligned to +/- strand): ref base A (66/50); new base C (3/8); total (69/58)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 132 total aligned reads displayed.

TAAAAGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGG  >  AP009048/2929433‑2929568

tAAAAGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGtt                                                                   <  1:4857065/71‑1 (MQ=255)

  aaaGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTa                                                                  <  1:718699/70‑1 (MQ=255)

  aaaGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTa                                                                  >  1:69276/1‑70 (MQ=255)

  aaaGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTAt                                                                 >  1:5027836/1‑71 (MQ=255)

   aaGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                <  1:7027248/71‑1 (MQ=255)

    aGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                >  1:5839882/1‑70 (MQ=255)

     gCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGa                                                              <  1:1203676/71‑1 (MQ=255)

      cGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATc                                                                >  1:586335/1‑68 (MQ=255)

      cGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCATCGATAAAGTTATCGaa                                                             <  1:4587840/71‑1 (MQ=255)

       ggTGAACGCCGCACGTATGGCGCTTCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGaaa                                                            >  1:3007212/1‑71 (MQ=255)

           aaCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCa                                                         >  1:5034058/1‑70 (MQ=255)

            aCGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATg                                                       >  1:46447/1‑71 (MQ=255)

             cgccgcACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATg                                                       >  1:1738498/1‑70 (MQ=255)

               ccgcACGTATGGCGCTGCGCCGTACCAGCGCGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTAc                                                    >  1:6694455/1‑71 (MQ=255)

               ccgcACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCAt                                                        <  1:197884/67‑1 (MQ=255)

                cgcACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGtt                                                                   >  1:6525127/1‑55 (MQ=255)

                   aCGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                <  1:4951919/71‑1 (MQ=255)

                    cGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:5446248/1‑70 (MQ=255)

                     gTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACa                                              <  1:712133/71‑1 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTTCGAAACAg                                             >  1:433729/1‑71 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAg                                             <  1:3538357/71‑1 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAg                                             <  1:4618812/71‑1 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAg                                             <  1:5963014/71‑1 (MQ=255)

                      tATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAg                                             <  1:835677/71‑1 (MQ=255)

                         ggCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGt                                                                    <  1:6773974/45‑1 (MQ=255)

                          gcgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTa                                          >  1:4266160/1‑70 (MQ=255)

                           cgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaa                                        <  1:2727544/71‑1 (MQ=255)

                           cgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTaaa                                        <  1:2718470/71‑1 (MQ=255)

                           cgcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTa                                          >  1:940751/1‑69 (MQ=255)

                            gcTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       <  1:1392129/71‑1 (MQ=255)

                               gcgcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGtt                                                                   <  1:3470631/40‑1 (MQ=255)

                               gcgcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGtt                                                                   <  1:5296372/40‑1 (MQ=255)

                                cgcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:3833386/1‑71 (MQ=255)

                                 gcCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAAc                                                           <  1:6614253/46‑1 (MQ=255)

                                  ccGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:1233290/71‑1 (MQ=255)

                                   cGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCa                                                         >  1:3277740/1‑46 (MQ=255)

                                   cGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:5545298/70‑1 (MQ=255)

                                    gTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATg                                                       >  1:5842242/1‑47 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:724734/1‑53 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:1532428/1‑53 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGaaa                                                >  1:1541778/1‑53 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGa                                      <  1:767731/63‑1 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGa                                      <  1:6028287/63‑1 (MQ=255)

                                     tACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACg                              >  1:1761350/1‑71 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGATACAGGTAAAGATATGAACGcc                            >  1:3390138/1‑68 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       >  1:6468716/1‑57 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAg                         >  1:2913585/1‑71 (MQ=255)

                                          gCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATAAACGCCAAg                         >  1:679844/1‑71 (MQ=255)

                                           cGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAAc                                               <  1:3830034/48‑1 (MQ=255)

                                           cGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAg                         >  1:5299627/1‑70 (MQ=255)

                                              gCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGc                             <  1:1068568/63‑1 (MQ=255)

                                              gCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGc                             <  1:871746/63‑1 (MQ=255)

                                               ccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatg                                    >  1:6996491/1‑54 (MQ=255)

                                               ccGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcg                    >  1:4747291/1‑71 (MQ=255)

                                                cgcgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATg                                  <  1:2528305/56‑1 (MQ=255)

                                                cgcgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   <  1:2721011/71‑1 (MQ=255)

                                                 gcgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:4570855/1‑70 (MQ=255)

                                                  cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGcc                            <  1:4002609/60‑1 (MQ=255)

                                                  cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgcg                  >  1:2886278/1‑70 (MQ=255)

                                                  cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgcg                  >  1:3682488/1‑70 (MQ=255)

                                                  cgcgTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGa                 <  1:7384204/71‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGCAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAAc              >  1:3255967/1‑71 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGat                                     <  1:3477806/48‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGat                                     <  1:5728343/48‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGat                                     <  1:620418/48‑1 (MQ=255)

                                                     gTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAAc              >  1:4961947/1‑71 (MQ=255)

                                                       cTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACct            >  1:7637019/1‑71 (MQ=255)

                                                        tGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACctc           >  1:5655541/1‑71 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAcc                     >  1:4811241/1‑60 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:747903/1‑62 (MQ=255)

                                                         gCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:2982016/1‑62 (MQ=255)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:1069425/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:6524508/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:7398388/65‑2 (MQ=25)

                                                          tttCGATAACTTTGTCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGTGAAAcg             <  1:2200057/65‑2 (MQ=25)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       >  1:7033860/1‑41 (MQ=255)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGa                                 <  1:2006492/47‑1 (MQ=255)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACctc           <  1:4130382/69‑1 (MQ=255)

                                                          ccTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACct            <  1:2203227/68‑1 (MQ=255)

                                                           ttCGATAACTTTATCGAGACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcg        >  1:7591951/2‑71 (MQ=21)

                                                           cTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcg        <  1:4448795/71‑1 (MQ=255)

                                                            tCGATAACTTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACcgc                   >  1:4956487/1‑59 (MQ=21)

                                                            tCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       >  1:4734270/1‑71 (MQ=255)

                                                            tCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       >  1:1726656/1‑71 (MQ=255)

                                                             cGATAACTTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:2014821/1‑42 (MQ=38)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGatat                                   >  1:7041894/1‑42 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                    >  1:3450165/1‑41 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGata                                    >  1:7862862/1‑41 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgcg      >  1:7886075/1‑71 (MQ=255)

                                                             cGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgc       >  1:6215451/1‑70 (MQ=255)

                                                              gATAACTTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTcgcg      <  1:5247900/70‑1 (MQ=255)

                                                              gATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACct            <  1:5956545/64‑1 (MQ=255)

                                                              gATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcgg     <  1:3310168/71‑1 (MQ=255)

                                                               aTAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGcc                            <  1:1392488/47‑1 (MQ=255)

                                                               aTAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcggc    <  1:6721552/71‑1 (MQ=255)

                                                                tAAAGTTATCGAAACCATGTACGAAACAGGTAAAg                                       <  1:1398492/35‑1 (MQ=255)

                                                                tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcggcg   <  1:1631415/71‑1 (MQ=255)

                                                                tAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCgcggcg   >  1:6041998/1‑71 (MQ=255)

                                                                 aaaGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTAcc                     >  1:7816489/1‑52 (MQ=255)

                                                                 aaaGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAAc              >  1:3373483/1‑59 (MQ=255)

TAAAAGCGGTGAACGCCGCACGTATGGCGCTGCGCCGTACCAGCGAGCCGCGCGTCTGCCTCGATAAAGTTATCGAAACCATGTACGAAACAGGTAAAGATATGAACGCCAAGTACCGCGAAACCTCTCGCGGCGG  >  AP009048/2929433‑2929568

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: `‑`