BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	3,804,808	C→A	6.7%	G107G (GGC→GGA)	yicJ →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	3,804,808	0	C	A	6.7%	428.6 / 8.1	133	G107G (GGC→GGA)	yicJ	predicted transporter
Reads supporting (aligned to +/- strand): ref base C (79/45); new base A (2/7); total (81/52)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 133 total aligned reads displayed.

TTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATCCCTTAC  >  W3110S.gb/3804740‑3804877

tttGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCa                                                                      >  1:2527925/1‑70 (MQ=255)

  tGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCaaaa                                                                   >  1:654152/1‑71 (MQ=255)

         cTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                            <  1:763437/71‑1 (MQ=255)

              aTTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcag                                                       <  1:1158094/71‑1 (MQ=255)

              aTTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcag                                                       >  1:2359110/1‑71 (MQ=255)

              aTTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcag                                                       >  1:1833858/1‑71 (MQ=255)

                tCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCt                                                             >  1:3719361/1‑63 (MQ=255)

                tCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagc                                                      >  1:1089962/1‑70 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagca                                                     >  1:233946/1‑70 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGaaa                                                    >  1:3596523/1‑71 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCaa                                                    >  1:1084039/1‑71 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCaa                                                    >  1:2067882/1‑71 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCaa                                                    >  1:2396968/1‑71 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCaa                                                    >  1:1656742/1‑71 (MQ=255)

                 cGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCaa                                                    >  1:3071221/1‑71 (MQ=255)

                  gggATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAt                                                   <  1:1823604/71‑1 (MQ=255)

                      tCGTCTGTGTACTGGCCTATAGCACGTCAGATCTCAGTATGAACGGCAAAATg                                                                 <  1:851979/53‑1 (MQ=255)

                      tCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttact                                               >  1:1185445/1‑71 (MQ=255)

                       cGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttact                                               >  1:2747207/1‑70 (MQ=255)

                           tgtgTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTacac                                          >  1:1873778/1‑71 (MQ=255)

                           tgtgTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTacac                                          >  1:48042/1‑71 (MQ=255)

                           tgtgTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGAAAAATGATTTATgcagc                                                      >  1:400059/1‑59 (MQ=255)

                             tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTacac                                          <  1:360422/69‑1 (MQ=255)

                             tgtACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTacac                                          <  1:3132552/69‑1 (MQ=255)

                                 cTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                             <  1:3532623/62‑1 (MQ=255)

                                    gCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtac                                 <  1:650151/71‑1 (MQ=255)

                                    gCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtac                                 >  1:1176894/1‑71 (MQ=255)

                                    gCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtac                                 >  1:3538276/1‑71 (MQ=255)

                                    gCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActt                                   <  1:3850815/69‑1 (MQ=255)

                                    gCCTATAGCACGCCAGATCTCAGTATGAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTActt                                   <  1:1760471/69‑1 (MQ=255)

                                      cTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                               <  1:2934468/71‑1 (MQ=255)

                                       tataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                               >  1:963489/1‑70 (MQ=255)

                                       tataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                               >  1:25382/1‑70 (MQ=255)

                                       tataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                               >  1:2840534/1‑70 (MQ=255)

                                       tataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                               >  1:151222/1‑70 (MQ=255)

                                       tataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacct                               >  1:1770008/1‑70 (MQ=255)

                                        ataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgctgca                                                     >  1:1368391/1‑47 (MQ=255)

                                        ataGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtacctta                             >  1:366383/1‑71 (MQ=255)

                                         taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTa                                                            <  1:1478601/39‑1 (MQ=255)

                                         taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagca                                                     <  1:358802/46‑1 (MQ=255)

                                         taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATgcagca                                                     <  1:3768353/46‑1 (MQ=255)

                                         taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactt                                              >  1:1106296/1‑53 (MQ=255)

                                         taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACAcc                                         >  1:1343964/1‑58 (MQ=255)

                                         taGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                            >  1:194311/1‑71 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACCCCCTa                                      >  1:1837825/1‑59 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTa                                      >  1:3271342/1‑59 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACtt                          >  1:2465710/1‑71 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACtt                          <  1:3383599/71‑1 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACtt                          <  1:2236532/71‑1 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACt                           >  1:541202/1‑70 (MQ=255)

                                           gCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACt                           >  1:2996658/1‑70 (MQ=255)

                                           gCACGCCAGATCTCAGTATCAACGGCAAAATGATCTATGCAGCAAttactta                                             >  1:3334697/1‑52 (MQ=255)

                                             aCGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        <  1:479364/71‑1 (MQ=255)

                                             aCGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        <  1:2446726/71‑1 (MQ=255)

                                              cGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        >  1:3368545/1‑70 (MQ=255)

                                              cGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        >  1:2681544/1‑70 (MQ=255)

                                                  aGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTtata                      >  1:692324/1‑68 (MQ=255)

                                                   gATCTCGGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACcgt                  >  1:3702066/1‑71 (MQ=255)

                                                   gATCTCAGTATGAACGGCAAAATGATCTATGCAGCAAt                                                   <  1:3314170/38‑1 (MQ=255)

                                                   gATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACcgt                  >  1:3714408/1‑71 (MQ=255)

                                                    aTCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                             <  1:1308581/43‑1 (MQ=255)

                                                    aTCTCAGTATGAACGGCAAAATGATCTATGCAGCAAttactta                                             <  1:2538141/43‑1 (MQ=255)

                                                    aTCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACcgtc                 >  1:2279824/1‑71 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAAttac                                                >  1:3145950/1‑38 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                            >  1:1437085/1‑58 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                            >  1:1084351/1‑58 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                            >  1:600278/1‑58 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                            >  1:1180884/1‑58 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTtat                       >  1:1230166/1‑63 (MQ=255)

                                                      ctcAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACcgtcgt               >  1:1896710/1‑71 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCgtcgtc              <  1:101789/69‑1 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCgtcgtc              <  1:3012643/69‑1 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCgtcgtc              <  1:264013/69‑1 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa            <  1:2185796/71‑1 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa            >  1:1775/1‑71 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa            <  1:1238926/71‑1 (MQ=255)

                                                         aGTATGAACGGCAAAATGATCTATGCAGCAATTACTTAAACCCTACTTACCTTACTTTATACCGTCGTCaa            <  1:1416690/71‑1 (MQ=255)

                                                           tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTActtaccttac                            <  1:1550873/53‑1 (MQ=255)

                                                           tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        <  1:2944864/57‑1 (MQ=255)

                                                           tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACcgtcgt               <  1:1159708/66‑1 (MQ=255)

                                                           tATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACCTTATACcgtcgt               <  1:77061/66‑1 (MQ=255)

                                                            aTGAACGGCAAAATGATCTATGCAGCAATTACTTACAACCTActtaccttac                            >  1:417891/1‑52 (MQ=255)

                                                             tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTa                                      <  1:121801/41‑1 (MQ=255)

                                                             tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa            >  1:2184002/1‑67 (MQ=255)

                                                             tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCaa            >  1:2889074/1‑67 (MQ=255)

                                                             tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAatat         <  1:3778497/70‑1 (MQ=255)

                                                             tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAatat         <  1:1662562/70‑1 (MQ=255)

                                                             tGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAatat         <  1:3218239/70‑1 (MQ=255)

                                                              gAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCATCGTCa             <  1:2957680/65‑1 (MQ=255)

                                                              gAACGGAAAAATGATTTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCATCGTCa             <  1:689017/65‑1 (MQ=255)

                                                               aaCGGCAAAATGATCTATGCAGCAATTACTTacac                                          >  1:2239415/1‑35 (MQ=255)

                                                               aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        >  1:2356868/1‑53 (MQ=255)

                                                               aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTta                        >  1:1410118/1‑53 (MQ=255)

                                                               aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATcc       >  1:643781/1‑70 (MQ=255)

                                                               aaCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACAGTCGTCAATATccc      >  1:755976/1‑71 (MQ=255)

                                                                aCGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCa             >  1:1648583/1‑63 (MQ=255)

                                                                  ggCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAata          >  1:2776034/1‑64 (MQ=255)

                                                                  ggCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAata          >  1:913335/1‑64 (MQ=255)

                                                                  ggCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATCCCTTa   <  1:2614584/71‑1 (MQ=255)

                                                                   gCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATCCCTTAc  >  1:3560575/1‑71 (MQ=255)

TTTGGCGCACTGCCATTCGGGATCGTCTGTGTACTGGCCTATAGCACGCCAGATCTCAGTATGAACGGCAAAATGATCTATGCAGCAATTACTTACACCCTACTTACCTTACTTTATACCGTCGTCAATATCCCTTAC  >  W3110S.gb/3804740‑3804877

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